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FOREWORD 


The  Florida  Agricultural  and  Mechanical  University  (FAMU)  in 
Tcdlahassee  was  the  site  of  the  16th  Annual  Meeting  and  the  20th  Day 
of  Scientific  Lectures  of  the  National  Society  of  Black  Physicists 
(NSBP)  on  April  21  to  24.  1993.  The  theme  of  the  meeting  was 
"Physics,  Technology  and  Business",  The  first  day  of  the  meeting  was 
held  in  the  Grand  Ballroom  on  FAMU's  beautiful  campus  and  the 
second  day  and  third  half-day  were  held  at  the  Radisson  Hotel  in 
downtown  Tallahassee.  Over  20  scientific  talks  were  given  by 
professional  and  student  physicists  during  the  meeting.  In  addition  a 
poster  session  was  held  on  Thursday  night.  On  Friday  a  panel 
discussion  on  the  subject  of  the  meeting's  theme  was  held  and  had  as 
its  title  "Minority  Owned  High  Technology  Businesses:  Techniques  for 
Growing  a  Successful  One"  and  was  led  by  Dr.  Donald  Butler  of  BLES 
Scientific.  The  banquet  was  held  on  Friday  night  at  the  Governor's 
Club  and  the  kejmote  speaker  was  Dr.  James  Stith  of  the  United  States 
Military  Academy,  West  Point.  The  Saturday  morning  session  featured 
talks  by  undergraduate  students.  On  Saturday  afternoon  a  tour  was 
held  of  Wakulla  Springs,  a  nearby  wildlife  preserve  park.  Dr.  Earl 
Shaw  was  elected  Fellow  of  the  National  Society  of  Black  Physicists. 
Dr,  Cynthia  McIntyre  received  the  NSBP  Outstanding  Dissertation 
Award  for  1993, 
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Turbulent  Enhancement  of  Second  Viscosity  in  Nonequilibrium  Flow 

Jean  Chabi  Orou*  and  Joseph  A.  Johnson  HI''' 

CeNNAs,  Florida  A&M  University,  Tallahassee,  FL  32310 

Abstract 

A  physical  model  for  entropy  production  associated  with  reactive  flow  is 

f 

presented.  The  general  form  j(f) «  (1-e  ^has  the  correct  asymptotic  behav¬ 
iors  for  t->0,  t->'C  and  t->oo  for  finite  t  where  t  is  the  relaxation  time.  From  this,  a 
relationship  between  the  second  viscosity  and  the  relaxation  time  for 
A+B->P  is  modeled  from  the  macroscopic  entropy  rate  equation.  Finally,  using 
a  first  Order  approximation  for  reaction  rate  distortion  from  reduced  molecular 
chaos,  enhanced  second  viscosities  can  be  predicted  for  strongly  turbulent  fluid 
systems  with  long  relaxation  times  and  V  •  V  0. 


Nonmenclature 

S  =entropy 

s  =specific  entropy 

Q  =quantity  of  heat 

T  =absolute  temperature 

e  =specific  internal  energy 

A  =area 

p  =density  of  the  fluid 

V  =velocity 

dv  =volume  element 

<1)  =dissipation  function 

a  =entropy  generated  per  unit  volume  per  unit  time 
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X  =relaxation  time 

X  =entropy  generated  by  chemical  reaction 

=derivative  of  the  component  of  the  velocity  vector  in  x  direction 
with  respect  to  x 
Re  =Reynolds  number 

t  =time  (laboratory  time) 

(I  =viscous  stress  tensor 

T|  ^dynamic  viscosity 

^  =second  viscosity 


Introduction 

It  is  usually  assumed  that  the  second  (or  bulk)  viscosity  and  the  chemical 
relaxation  time  are  connected.  However,  the  second  viscosity  is  only  relevant  in 
processes  where  there  is  a  compression  or  an  expansion  and  it  is  usually  re¬ 
garded  as  small  even  when  it  is  relevant.^  Nonetheless,  it  now  seems  appropri¬ 
ate  to  determine  in  general  whether  or  not  there  can  be  circ\amstances  where  the 
influence  of  second  viscosity  might  be  important.^  This  seems  especially  interest¬ 
ing  since,  in  some  situations,  effects  due  to  turbulence  are  often  treated  as  effects 
due  to  anomalous  viscosity.  The  present  paper  uses  the  macroscopic  entropy 
production  rate  equation  in  a  reactive  flow  in  order  to  approach  the  possibility  of 
a  connection  between  the  second  viscosity  and  the  relaxation  time.  In  addition, 
by  using  the  qualitative  implications  of  recent  experimental  results,  we  will 
speculate  on  a  role  for  turbulence  in  these  phenomena.  Generally  stated,  this 
model  is  proposed  on  intuitive  arguments. 

Entropy  and  a  Relaxation  Process 
The  total  change  of  entropy  in  a  system  is 
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o  c  —  f  dQrev 

^2-^1  -  J — T 

1 

where  we  have  assumed  the  system  is  brought  from  state  1  to  state  2  so  that 

S2  —  Si  =  ASo  +  ASi, 

The  first  term  in  the  right  is  the  entropy  carried  into  the  system  over  the  bound¬ 
aries  from  outside  and  the  second  term  is  the  entropy  produced  in  the  system 
during  the  process. 

In  what  follows,  we  are  going  to  deal  with  a  system  in  which  the  entropy 
carried  into  the  system  is  zero  and  the  last  term  in  the  above  equation  is  gener¬ 
ated  by  viscosity,  thermal  conduction  and  chemical  reaction.  Specifically,  we 
define  a  to  be  the  entropy  per  unit  time  and  to  be  the  dissipation  function. 
Specifically,  O  can  be  interpreted  as  the  irreversible  dissipation  of  the  mechanical 
energy  into  heat  caused  by  the  viscosity  per  unit  time  per  unit  volume. 

With  these  definitions,  the  energy  balance  can  be  expressed  as  follows:^ 

in p[^  -  = UK*  -  Mm'’ 


where 


O  =  div{V}i)-VdiviJL 


^  du  dv  dw 


+M, 


f  dv 


dz 

^dw  du\ 
V 


f  du 
dy  dx 


V 


'\dz  dy)  "^\dx  dz. 

in  which  the  velocity  vector  V  has  components  (w, v,  w),  h  is  the  viscous  stress 


tensor  and  q  is  the  energy  flow  vector.  The  previous  equation  becomes: 

pf  De 


Dt  n2  Dt 
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Since 


TdS  =  de- 


f  ^ 


we  now  have: 


where  ^  = 


When  the  entropy  generated  by  a  chemical  reaction  is  taken  into  account, 
an  additional  term  is  required  in  the  expression  of  the  entropy  source.  The  new 
expression  for  the  entropy  becomes 


Tds  =  de- 


(z'' 


dp  +  Tdx 


where  F 


=  vanishes  if  the  gas  attains  an  unconstrained  equilibrium  at 


all  time. 

In  order  to  derive  a  relationship  between  the  additional  term  and  the  en¬ 
tropy  source,  one  notices  that  the  expression  for  the  entropy  source  becomes: 

°-r  V 

We  will  limit  ourselves  to  the  case  of  a  chemical  reaction  of  the  type  A+B — >P 
where  A,  B  and  P  are  assumed  to  be  nonmonoatomic  perfect  gases^  and  %  is  the 
entropy  generated  by  the  chemical  reaction  during  the  process.  Because  the  reac¬ 
tion  stops  after  a  finite  time,  %  must  reach  a  limit  value  which  is  its  value  at  equi¬ 
librium  and,  at  that  time,  the  derivative  of  x  with  respect  to  time  becomes  zero  . 

On  the  other  hand,  the  process  being  irreversible,  the  entropy  must  in¬ 
crease  from  zero  to  the  limit  value  in  a  way  which  agrees  with  the  active  relax¬ 
ation  process.  Assuming  that  %  is  only  a  function  of  time  and  that  it  satisfies  the 
above  statements,  the  profile  of  %  could  be  a  function  of  the  form: 
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The  general  behavior  of  x=x(t)  during  a  relaxation  process  would  be  as  is  shown 
in  Fig.  1  where  x  is  the  relaxation  time.  We  also  assume  that  Boltzmann’s  equa¬ 
tion  holds.5 


Second  Viscosity  and  Relaxation  Time 

Let's  consider  a  chemical  reaction  which  starts  at  t=0  with  %  =0. 

/ 

Elaborating  the  approximation  above,  we  assume  that  the  value  of  x  becomes 

t 

proportional  to  (l-(l/e))  at  t=T  and  stops  when  t  is  large  enough  so  that^  r  is 
almost  zero.  Whatever  value  x  has,  provided  it  is  a  finite  one,  equilibrium  will 
eventually  be  reached.  If  the  entropy  at  equilibrium  is  independent  of  x,  then  Xeq 
is  independent  of  x  and  the  rate  of  change  of  x  is  fixed.  The  expression  of  the  en¬ 
tropy  source  becomes: 

O  q.VT  pT  j. 

a  — - - - +  e  r 

T  Tt 

The  last  term  is  the  entropy  generated  per  unit  time  per  unit  volume  by  the 
chemical  reaction . 

For  t  fixed  and  a  small  relaxation  time,  the  equilibrium  is,  of  course, 
reached  more  quickly  than  when  the  relaxation  time  is  large.  For  one-dimen¬ 
sional  flow,  one  now  finds: 

Next  we  restrict  the  treatment  to  an  isolated  system  at  a  constant  temperature; 
this  restriction  ignores  heat  released  (or  absorbed)  by  the  process  A  +  B  ->P .  The 
entropy  source  then  becomes: 

:  pT  .L 
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Thus  we  achieve  the  previously  derived  energy  balance  equation  takes  the 


form; 


///P 


De  P  Dp 
Dt  Dt 


dv  =  JJJ 


2  pr 


dv 


V  J  . 

By  restricting  our  treatment  to  a  perfect  and  nonmonoatomic  gas  at  a  constant 
temperature  in  an  isolated  system,  the  integrands  in  the  energy  balance  equation 
cannot  depend  on  time  macroscopically.  This  means  that  the  derivative  of  a  with 
respect  to  time  must  be  zero  and,  consequently,  the  entropy  generated  by  the 
molecular  process  must  be  balanced  by  a  change  in  the  viscosity.  Since  the  sec¬ 
ond  viscosity  alone  is  free  to  show  a  change  with  time,  it  is  constrained  by  the 
energy  balance  equation  above  (since  do/dt  =  0)  as  follows: 

^_pr  1 

dt 


2  2 

U  T 

X 


These  results  are  summarized  in  Figs.  2  and  3.  Setting  b^=(prxeq/ux^),  we 

t 

show  ^  ^  range  in  t  and  x  indicated,  a  substantial  sensi- 

dt  X 

tivity  in  d^/dt  is  observed;  this  is  particularly  true,  for  example,  at  low  velocities 
Ux  and  high  densities  p.  The  relative  magnitudes  are  indicated  in  Fig.  3  which  is 
a  slice  through  Fig.  2  at  x=30  |isec. 

Notice  that  as  the  relaxation  time  decreases,  the  overall  entropy  change  due 
to  the  chemical  reaction  increases.  For  decreasing  relaxation  times,  there  corre¬ 
spondingly  (from  the  equations  above)  is  a  decreasing  role  for  the  evolution  of 
second  viscosity  to  play  in  the  entropy  balance. 

Second  Viscosity,  Relaxation  Time,  and  Turbulence 
It  has  been  observed  that  fully  developed,  as  well  as  transitional,  turbu¬ 
lence  can  be  described  by  qualitatively  predictable  Orr-Sommerfeld  like  behav- 
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iors  in  the  fonn  I=I(Re,  Repeak)  where  I  is  the  hirbulent  intensity  and  Repeak  is  the 
Reynolds  number  at  maximum  turbulent  intensity.^  It  is  also  observed  that  the 
reaction  rate  decreases  with  increasing  turbulent  intensity 7  For  x  —  x(Re)  and  I  - 


I(Re),we  get  x  =  x(Re(I))  implicitly  where  I  is  the  turbulent  intensity.  Notice 
T  =  cnst .  Then  the  form 


is  a  suitable  intuitive  first  approximation. 

This  expression  is  now  combined  with  the  expression  above  for  d^/dt  to 
achieve  the  behaviors  shown  in  Figs.  4  and  5.  In  Fig.  4  one  notices  that  the  range 
of  non-zero  turbulent  intensities  increases  as  the  value  of  Rcpeak  increases.  A 
typical  value  for  Repgak  would  be  roughly  10^;  using  this  value  and  a  value  of  x 
=30  iisec,  one  obtains  the  range  of  values  for  d%/dt  shown  in  Fig.  5.  Here  as  in 
Fig.  2  the  most  dramatic  changes  are  seen  at  very  low  values  of  laboratory  time. 
Nonetheless,  values  for  t  (roughly  1  jisec  <  t  <  40  |isec)  are  found  during  which  a 
significant  influence  of  relaxation  time  and  turbulent  intensity  on  the  evolution 
of  second  viscosity  is  possible. 

Conclusions 

By  using  a  physical  model  for  entropy  production  associated  with  reactive 
flow,  a  relationship  between  the  second  viscosity  and  the  relaxation  time  is  de¬ 
termined  from  the  macroscopic  entropy  production  rate  equation.  This  relation¬ 
ship  enables  us  to  predict  the  behavior  of  the  second  viscosity  when  turbulence 
evolves  in  a  nonequilibrium  flow.  Taking  advantage  of  previous  studies  which 
have  concluded  that  the  reaction  rate  decreases  when  hirbulence  increases,  one 
can  argue  that  the  second  coefficient  of  viscosity  increases  when  turbulence  in¬ 


creases. 
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It  is  difficult  to  normalize  X=X(>);  therefore  the  results  obtamed  here 
second  Viscosity  as  a  function  of  the  relaxahon  time  and  the  Reynolds  n^ber 
give  only  qualitative  behaviors.  Nonetheless,  a  relationship  has  been  estabbshed 
between  the  second  viscosity  and  the  relaxation  dme  for  a  nonequilibrium  pro- 
cess  We  have  shown  how  this  relationship  can  lead  to  a  dependence  o  e  sec 
ond  viscosity  on  turbulence.  For  a  flow  where  V  •  V  ^  0,  anomalous  v.cous  e  - 
facts  can  therefore  be  expected,  under  some  circumstances,  for  turbuen 

nonequilibrium  systems. 
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Figure  2.  The  Serisitivity  of  the  Rate  of  Change  in  Second  Viscosity  (d^  dt)  to 
Changes  in  Laboratory  Time  (t)  and  Chemical  Relaxation  Time  (t). 


15 


arb.  units 


t  (sec) 


Figvire  3.  A  Slice  through  the  Coi\tour  Plot  in  Fig.  2  at  t  -  30  ^sec. 
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(arb.  units) 


Figure  4.  The  Sensitivity  of  Turbulent  Intensity  to  Changes  in  Reynolds 
Number.  Turbulent  Intensity  is  indicated  by  Int  and  Repeak  is  the  value  of 
the  Reynolds  number  at  peak  txirbxilent  intensity  as  discussed  in  Ref.  6. 
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Re 


Contour  Plot: 


1 

dt 


vs,  (Re  ,  t ) 


Figure  5.  The  Sensitivity  of  the  Rate  of  Change  in  Second  Viscosity  (d^/dt)  to 
Changes  in  Laboratory  Time  (tiab)  and  Reynolds  Number.  A  value  of 
Rcpeak  =1.0xl06  is  used  for  this  contour  plot. 
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ABSTRACT 


Addition  theorems  are  derived  for  Coulomb  Sturmian  fimctions  and  reduced  Coloumb 
Sturmian  functions  in  coordinate  and  momentum  space.  The  Coulomb  Sturmians  and  their 
Fourier  Transforms,  are  complete  in  the  sense  of  distributions  in  Sobolev  space  W2^^(R^), 
which  is  a  proper  subset  of  Hilbert  space  L^(R^).  The  derivations  involve  the  use  of 
the  relationship  between  Coulomb  Sturmians  and  four  dimensional  spherical  harmonics 
(Fock  hyperspherical  projection)  and  proceeds  by  representing  the  translation  operator 
in  coordinate  and  momentum  space,  on  a  basis  of  Coulomb  Sturmians  and  their  Fourier 
Transforms  respectively.  The  resulting  matrix  representations  of  the  translation  operator 
in  coordinate  and  momentum  space  axe  respectively,  finite  linear  combinations  of  Coulomb 
Sturmians  and  their  Fourier  Transforms,  and  are  both  imitary.  The  addition  theorems  are 
useful  for  application  to  problems  which  involve  the  evaluation  of  multicenter  integrals 
over  exponentially  decaying  orbitals. 
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I.  INTRODUCTION 


The  problem  of  performing  multicenter  integrals  is  endemic  to  the  solution  of  many 
problems  of  functional  analysis.  The  solution  may  be  approached  by  the  method  of  numer¬ 
ical  quadrature  or  by  analytical  analysis.  The  present  paper  is  concerned  with  the  method 
of  anzdytical  analysis.  Given  this  delimitation,  two  further  possible  approaches  are  gen¬ 
erally  available:  either  a  transformation  of  variables  is  sought  or  an  addition  theorem  is 
used.  However,  it  is  fequently  the  case  that  numericzil  quadratures  axe  combined  with  vari¬ 
able  transformations  and/or  addition  theorems.  The  literature  on  these  techniques  is  large 
and  diverse  and  no  attempt  will  be  made  here  to  review  or  contrast  these  methods.  The 
proceedings  of  a  conference  on  multicenter  integrals^  and  a  recent  review  paper^  provide 
a  prospective  on  the  literature. 

The  present  work  is  concerned  with  the  derivation  of  addition  theorems.  Such  addition 
theorems  may  be  generally  defined  by 

nfur^)  =  E  E  (1) 

»  i 

If  no  integration  over  one  of  the  vectors  (say  ^2)  is  to  be  performed,  Eq.  (1)  is  sometimes 
rewritten  as 


(2) 

I 

where 

Bi(f2)  =  ECyJ>(r-,)  (3) 

i 

Thus,  an  addition  theorem  explicitly  separates  the  variables  n  and  r2  so  that  integrations 
may  be  performed  over  one  or  both  of  them. 

The  addition  theorems  indicated  schematicaiUy  in  Eqs.  (1)  and  (2)  are  to  be  distin¬ 
guished  from  the  Lowdin  a-function  (or  a-fimction)  technique,^  which  is  defined  by 

(4) 

•  i 

The  characteristic  which  distinguishes  the  a-function  method  (Eq.  (4))  from  the  present 
method  (Eq.  (1))  is  that  the  addition  theorem  in  the  a-function  method  has  an  assumed 
anguleir  behavior  (spherical  harmonics  Yi^{f))  and  the  radial  behavior  is  then  worked  out 
exactly  (a-function).  The  a-function  has  assumed  many  forms  which  vary  in  detail  (see 
for  example  Refs.  3  -  6),  but  which  are  very  similar.  Some  of  these  methods  are  not 
called  a-function  methods,  but  they  share  the  fundamental  form  of  Eq.  (4).  One  key 
characteristic  is  the  two-range  nature  of  the  a-function  and  the  slowly  convergent  infinite 
series  which  obtain.^"® 
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The  problem  of  performing  multicenter  integrals  is  especially  important  in  quantum 
chemistry  when  the  LCAO  (linear  combination  of  atomic  orbitals)  method  is  used.^** 
In  the  LCAO  method,  wavefimctions  for  molecules  are  constructed  as  a  solution  of  the 
Schrodinger  equation  by  taking  linear  combinations  of  atomic  orbitals  with  coordinate 
centers  at  various  locations  in  space.  In  pau’ticulax,  the  atomic  orbitals  are  frequently 
placed  at  the  origin  of  the  atomic  nuclei  that  compose  the  molecule.  It  is  known  that  the 
convergence  of  basis  set  expansions  is  critically  controlled  by  how  well  the  basis  functions 
satisfy  the  boundary  conditions  which  the  wavefunction  is  required  to  satisfy.  There  are 
fundamentally  three  such  requirements:  (1)  long-range  exponential  decay;^^  (2)  electron- 
nucleeir  cusp;^^  and  electron-electron  cusp.^*  The  basis  functions  used  by  most  quantum 
chemists  (Gaussians)  satisfy  none  of  the  three  requirements.  They  are  chosen  because  of 
the  relative  simplicity  of  the  multicenter  integrals  which  result. 

The  proper  atomic  orbitals  $  are  well  known  to  be  composed  of  linear  combinations 
of  Slater  Type  Orbitals  (STOs) 

^  (5) 


where 


and 


(6) 


Nn,{a)  = 


|(2n,)!li 


(7) 


The  function  represented  by  5  is  a  normalized  STO.  The  STOs  satisfy  requirements  one 
and  two. 

The  diffictilty  of  performing  multicenter  integrals  when  the  STOs  are  used  as  the 
atomic  orbitals  is  well  known®  and  led  to  the  use  of  Gaussian  Type  Orbitals  (GTOs).^°  The 
GTOs  have  a  dependence  on  r  that  involves  a  form  exp(— ar^),  as  opposed  to  the  STO  form 
exp{—ar).  Although  the  multicenter  integrals  are  easier  to  perform  than  the  corresponding 
integrals  are  for  STOs,  they  suffer  from  the  fact  that  the  cusp  at  the  nuclei  and  the 
proper  exponential  long-range  behavior,  are  difficxilt  to  represent  when  the  wavefunction 
is  calculated  from  linear  combinations  of  GTOs.^^ 

It  is  the  intent  of  the  present  work,  to  suggest  another  set  of  exponentially  decaying 
functions  (Coulomb  Stxirmians^^),  related  to  the  STOs,  which  have  a  number  of  properties 
which  aid  in  the  evaluation  of  multicenter  integrals,  when  the  integrand  is  composed  of 
them.  Further,  an  addition  theorem  is  derived  which  facilitates  the  derivation  of  explicit 
analytic  formulas  for  the  evaluation  of  the  multicenter  integrals.  In  section  II,  some  prop¬ 
erties  of  the  Coloumb  Sturmians  (CSs),  which  are  useful  for  the  derivation  of  the  addition 
theorem,  for  evaluating  special  cases  of  the  addition  theorem,  or  for  evaluating  multicen¬ 
ter  integrals  over  CSs,  are  presented.  In  section  III,  the  addition  theorem  is  derived  for 
the  CSs.  In  section  V,  procedures  useful  for  evaluating  certain  multicenter  over  CSs  are 
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discussed.  In  section  V,  certain  special  cases  of  the  addition  theorem  and  multicenter 
integrals  are  discussed.  Finally,  section  V  gives  the  summary  and  conclusions. 


II.  SOME  PROPERTIES  OF  COULOMB  STURMIANS 


An  important  set  of  exponentially  decaying  functions  was  first  formally  introduced 
by  Rotenberg,^^  in  1962,  which  he  called  Sturmians.  These  fimctions  are  more  commonly 
called  Coulomb  Sturmians  now^^  to  distinguish  them  from  other  types  of  Sturmian  func¬ 
tions  with  different  potentials.  An  important  feature  of  the  CSs  is  that  they  axe  a  complete 
set  in  Sobolev  space^^>^®  which  is  a  proper  subset  of  Hilbert  space 

Thus  completeness  in  W2^^^(R®)  implies  completeness  in  as  explicited  pointed  out 

by  Weniger.^^  The  CSs  axe  defined  by^^ 

=  Nn,(a)e-"L^;!*i\(2ar)yr(2af )  (8) 

where 


Nni(a)  =  (2a)3/2 


(n  —  /  —  1)! 
2n(n  -1-  /)! 


ii/2 


(9) 


is  the  normalization  constant,  L  is  a  Laguerre  polynomial  (same  definition  as  given  in 
Arfken^^)  given  by 


+  (10) 

where  iFi  is  a  confluent  hypergeometric  function^®  and  is  a  solid  spherical  harmonic 
given  by 


yr(rj  =  r'Yr(f)  (n) 

and  K  is  a  regular  spherical  harmonic.*®  The  series  for  L  terminates  in  the  present  case 
and  gives 


r(2/+l) 


(2ar) 


n—/— 1 

Y,  Brni{2arY 

r=0 


(12) 


(-l)^(n+i)! 

(n  —  I  —  I  —  r)!(2/  -|-  1  -f  r)!r! 


(13) 


In  particular,  the  Coulomb  Sturmians  are  complete  with  a  purely  discrete  spectrum. 
It  is  easy  to  show  that  a  normalized  STO  (S)  can  be  expanded  in  a  terminating  series  of 
CSs  (^)'® 
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(14) 


n/ =/,*+! 


with 


‘2anf{U  +  —  U  - 

(2ni)! 


1 

2 


X 


-If- 

E 


<r=:0 


(-l)‘^(n,-  +  /,•  +  <7)! _ 

{nf-li-l-a)\(2li  +  l  +  <r)\(7\ 


(15) 


and  vice  versa 


ni 


^{nl) 


.{a,r)=  ' 

n/=/,-hl 


("').(a)S”^,{a.f) 


(16) 


4";>‘(«)  = 


(«j  -  ^  -  l)!(ni  +  li)!'  ^ 

[(n/)!l^ 


arii 


X  (-I)"/-'*-! 


(ni-n/)!(n/  +  /i)!(n/-/i-l)! 


(17) 


For  reasons  that  will  be  discussed  in  section  IV  below,  it  is  useful  to  “equalize”  the 
screening  parameters  in  one  electron,  single  center  densities,  composed  of  CSs  basis  func¬ 
tions.  To  accomplish  this,  it  is  first  useful  to  equalize  the  screening  parameters  for  densities 
composed  of  STOs.  Some  trivial  algebra  produces 


■Vni(Q^)-Vn2(/^)  >ymim2 


(7,7,^ 


where  7  =  and  T  is  the  STO  one  electron  density  defined  by 


(18) 


^mi  m2 
^(n0i(n/)2 


iaj,r)  ^  SlZl 


Then,  defining  the  one  electron,  single  center,  CS  density  by 


(19) 


(20) 


the  required  expression  (which  equalizes  the  CS  screening  parameters  for  the  density  Q), 
is  obtained  by  first  using  Eq.  (16)  for  each  CS  in  then  using  Eq.  (18)  to  equalize 
the  screening  parameters  in  T,  and  then  using  Eq.  (14)  to  recover  the  required  Q  with 
equalized  screening  parameters.  The  result  is 
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(21) 


Til  1\T  f  \  ^ 

n3=:/i4“l  ns— /i+l 

-  t  t  ‘S'- (7) 


n4=/s+I 

It  is  important  to  note  that  all  of  the  sums  in  Eq.  (21)  represent  terminating  series.  The 
b  coefficients  are  given  by  Eq.  (15)  and  the  c  coefficients  are  given  by  Eq.  (17). 

The  CSs  satisfy  a  partial  differential  equation^"^  given  by 


2an 


This  can  be  contrasted  with  the  partial  differentied  equation  satisfied  by  the  STOs 

[Vl  -  a2]5^(a,  f)  =  Vni{a,  r)5r,(a,  r) 


where 


(22) 


(23) 


(24) 


It  can  be  easily  demonstrated  that  the  CSs  satisfy  a  homogeneous  Fredholm  integral  equa¬ 
tion  of  the  second  kind^^ 


diG(a,|r  -  i|)-^^;[)(or,x)  (25) 

where  G  is  the  Green’s  function  (same  as  for  the  modified  Helmholtz  equation)  and  is 
given  by 


G{a,  |r-i|) 


^1 

J\ 


The  Green’s  function  satisfies  the  partial  differential  equation 


(26) 


[Vl  -  a‘^]G{a,  If  -  x|)  =  -S{f-  x)  (27) 

The  trcinslation  operator  (note  that  hcis  to  be  expressed  in  Cartesian 

coordinates  for  this  expression  to  be  valid)  is  defined  by 

fjFir)  =  F{f+  A)  =  F{fA)  (28) 

where  =  r  +  A.  Now  an  important  characteristic  of  the  Laplacian,  which  is  well  known, 
is  its  invariance  with  respect  to  translations 

-  =  Vl  (29) 

?+A  r  ^ 
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Applying  the  translation  operator  to  both  sides  of  Eq.  (22)  and  using  Eq.  (29)  yields 


2an 


[V^-  -  a2]^m(a,r-‘+  A)  = 


(30) 


It  may  be  noted,  in  passing,  that  Eq.  (30)  can  be  used  to  generate  an  algebraic  eigenvalue 
equation  by  expanding  and  letting  V^operate  Then 

premultiplication  by  used  and  finally,  integration  over  dr.  The  resuh  is  an 

eigenvalue  equation  in  which  the  exact  eigenvalues  are  known.  However,  the  coefficients 
(eigenvectors)  in  the  aforementioned  expansion,  which  result  from  the  diagonalization  of 
the  coefficient  matrix,  are  the  desired  quantities  and  unfortunately,  these  components  have 
first  order  errors  in  them,  as  opposed  to  second  order  errors  in  the  approximate  eigenvalues. 
Thus  this  procedure  proves  too  slowly  convergent^®  to  be  useful,  without  some  systematic 
way  of  correcting  the  eigenvectors. 

The  integral  equation  solution  to  Eq.  (30)  can,  nevertheless,  be  written 


* J,(a,r-+  A)  =  jdx  G(a,  |f-  -4) 


(31) 


This  can  be  contrasted  with  the  result  of  applying  the  translation  operator  directly  to  both 
sides  of  Eq.  (25)  to  give 


dx  G{a,  \f+A-  xD— 

Finally,  the  CSs  satisfy  the  orthonormality  condition^® 


/ 


and  the  closure  condition^^ 


nlm 


(32) 


(33) 


(34) 


where  a;n(«^,^)  =  ^  overlap  of  two  CSs  is  given  by^®’^® 


J  df  ^*;;;)^(a,r  )'i’p)^(a,r  )  =  5(/,„),(/mh[7;l,/,^n,,n,+l  +^nt,n,  +7„,l/ni,n2-l]  (35) 


where 


±  _  p(n_T0(n±J_±_l)ni/2 

*  ““  I  7Z7znrT\  j 


(36) 


4n(n  ±  1) 

Note  that  the  integrands  of  Eqs.  (33)  and  (35)  are  composed  of  CSs  with  the  same  screening 
parameters.  The  expression  given  in  Eq.  (35)  is  derived  from  Eq.  (33)  using  the  following 
recursion  relation; 
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(37) 


ar  = 


7;^/ 


(n  +  l) 


X  «  +  (n-^ 

These  equations  will  prove  useful  in  the  discussions  given  below. 


III.  COULOMB  STURMIAN  ADDITION  THEOREM 


Filter  and  Steinborn^®  have  given  an  addition  theorem  for  orthonormalized  Slaters 
(ONSs,  orthonormedized  functions  built  up  from  the  STOs  by  analytical  Grzim-Schmidt 
orthonormalization)  by  using  the  translation  operator  explicitly  and  deriving  a  formula 
for  the  matrix  elements  of  the  translation  operator.  These  results  have  been  applied 
to  the  calculation  of  the  groimd  state  wavefunction  of  with  somewhat  disappointing 
results.^^  The  results  were  disappointing  because  of  the  slowness  of  the  convergence  of  the 
series  expressions  that  resulted  from  the  use  of  the  addition  theorem. 

The  present  work  derives  an  addition  theorem  for  the  CSs,  similar  in  structure,  to  the 
addition  theorem  of  Filter  and  Steinbom  for  the  ONSs.  The  present  derivation,  however, 
is  quite  different  in  approach.  As  is  presented  below,  the  present  work  exploits  the  Fourier 
transform  of  a  and  the  relationship  of  the  Fourier  transform  to  four-dimensional 

sphericcil  harmonics^®  using  the  Fock  hyperspherical  projection.^®  The  method  is  similar 
in  many  respects  to  techiniques  of  Shibuya  and  Wulfinaan^^  and  to  the  work  of  Alper^®  and 
Novosadov.^®  The  differences  between  these  techniques  eind  the  present  technique  will  be 
implicit  in  the  derivation  presented,  eind  in  the  approach  by  which  the  addition  theorem 
is  used  to  evaluate  the  multicenter  integrals.  The  work  of  Alper^®  is  perhaps  most  closely 
related  to  the  present  methods.  The  differences  and  improvements,  compared  to  the  work 
of  Alper,  will  be  explicitly  discussed  below. 

The  symmetric  definition  of  the  three-dimensional  Fourier  transform  is  first  appHed 
to  a  CS  so  that 


f  dp  e‘^'^  ^'„7(«.^ 

(38) 

[  dfe~'P'^  ^;0(a,r) 

(39) 

where  ^*^7  the  Fourier  transform  of  the  CS  (CSFT).  Weniger^®  demonstrated  the  re¬ 
lationship  between  a  CSFT  and  a  four-dimensional  spherical  harmonic  (the  FDSH  is  the 
one  used  by  Shibuya  and  Wulfman,  Ref.  27  and  Alper,  Ref.  28)  given  by 


X-,(r,^,0)  =  (-l) 


\-i)‘Trnt2ir^ 


+  p2)2 


(2r)5/2 


(40) 
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where  the  relationship  between  r  (hyperangle)  and  a  (screening  parameter)  is  defined  by 


i  =  sin{T)sin{d)cos{<l>) 

+p^ 

rj  =  ■  =  sin{T)sin{9)sin{<f>) 

+  p 

^  =  sin{T)cos{9) 

+p^ 


(41a) 

(416) 

(41c) 

(41d) 


Thus  a  is  a  scaling  parameter  and 

«^  +  ^"  +  C"+x’  =  i  (^2) 

It  is  clear  then  that  a  point  p  of  three-dimensional  space  is  mapped  (one-to-one)  onto  a 
surface  of  a  four-dimensional  unit  sphere  described  by  the  angular  vanables  The 

following  relationships  are  then  determined. 


sin 


X  r  2ap  ii 

(^)  =  iTJx;;?  J 


2a 


dO,  =  sin^{T)sin{9)dTd6d<j>  =  [^2  +p2] 


(43) 


(44) 


Do  not  confuse  the  four-dimensional  spheric^  differential  surface  element  dO.  with  the 
single  center  charge  density  of  Eq.  (20).  The  quantity  7r„,  has  been  defined  in  various  ways 
(as  pointed  out  by  Weniger^®).  For  example,  Biedenham  uses  7r„/  =  (-i)  and  Shibuya 
and  Wulfman^^  and  Alper^*  use  7r„,  =  1.  In  any  event,  |7r„,|  =  1.  The  orthonormality 
conditions  for  the  CSFTs  and  FDSHs  are  given  by 


1‘ir 

dr  sin} {r)  I  d9  sin{9)  /  d<i>  X*jJ^^iT,9,(f>)Xj^]^^{T,9,<f>) 
Jo  Jo 

=  j  dp%Z'  (“■  ?  ?  ) 

—  27r  6(nJm)i(nlm)2 


(45) 


The  convention  of  ir„/  =  1  will  be  used  from  this  point  on.  The  Shibuya  and  Wulfman 
expansion  of  a  plane  wave  in  terms  of  CSs  is  used  and  is  given  by 


=  4(7ra)^/^ 


n/m 


(a^  +  p2 ) 


(46) 


Eq.  (44)  may  be  compared  with  the  expansion  of  a  plane  wave  in  terms  of  CSs  and  CSFTs 
given  by  Weniger^® 
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(47) 


n/m 

For  later  use,  the  Fourier  transform  representation  of  a  three-dimensional  Dirac  delta 
function  and  the  Coulomb  potential  may  be  written  as 


-Pi)  =  (2^ 

1  1  /“ 

|ri-r2|  27r2  y  P  p2 


(48) 

(49) 


The  derivation  of  the  addition  theorem  begins  by  applying  the  translation  operator 
defined  by  Eq.  (28)  to  both  sides  of  Eq.  (38).  The  result  is 


$;;*,(a,r -b  A  )  =  (27r)-=*/2  j  dp  ^„7(a,^  (50) 

The  Shibuya  and  Wulfman  plane  wave  exp^lnsion  of  Eq.  (46),  is  used  for  both  plane  wave 
terms  (separately  expanded)  in  the  integrand  of  Eq.  (50).  Then  Eq.  (44)  is  used  to 
transform  the  volume  integral  over  p  to  an  integral  over  the  surface  of  a  four-dimensional 
unit  sphere  with  the  result 

^;;)(a,r-b44)=  E  [2(-l)’’‘+'**-"-'‘-'*-'-'i''+'’"'7r"®/"] 

{nlm)i  (n/m)2 

At  the  stage  of  Eq.  (51),  there  are  six  infinite  sioms.  Unless  some  of  the  sums  can  be 
terminated,  the  results  that  derive  from  Eq.  (51)  would  seem  to  of  limited  use.  This 
problem  will  be  alleviated  below.  Then,  a  critical  step  is  to  use  the  result^”^'^* 


(2Xf„(n)  +  Xj„(n)]  (52) 

where  Xi%(n)  =  1  and  X%{n)  =  cos(r).  Thus,  using  Eq.  (50)  in  Eq.  (49)  produces 


’^n/(«7^+  >1  )  =  Q! 


=  ^5/2 


(n/m)i  {nlm)2 


where 

^(nl)i(nl)2  ^  ^ 


(54) 
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(55) 


where  Xiq{CI)  =  1  h2is  explicitly  been  used.  It  should  be  noticed  that  C  and  J  in  Eq.  (53) 
are  explicitly  independent  of  the  screening  constant  a  and  the  translation  vector  A.  The 
relation 


(56) 

is  also  of  use.  The  property  of  the  FDSHs  which  allows  the  product  of  two  FDSHs  on  the 
same  hypersurface,  to  be  written  as  a  finite  linear  combination  of  a  single  FDSH  on  the 
same  surface,  is  centred  to  the  present  strategy.  This  expression  was  given  by  Shibuya  and 
Wulfman^’^  and  tised  to  great  advantage  by  Alper.^®  It  may  be  expressed  as 


— 1 

n3=|Tii  — n2|+l,2  ?3=Ui“^2|t2 

The  \2’  on  the  sum  indicates  every  other  term  is  used.  The  C  is  defined  by 

C(ni/imi,n2/2^T^2;^3^3^3)  =  +  l)(2/2  + 

(n2-l)  (ni  -1)  (na-l)  ^ 

X  {I2lim2mi\l2lihm3)  { 

h  h 


(57) 


2 

h 


(58) 


The  (•  •  •)  term  is  a  3j  symbol^”  and 


is  a  Wigner  9j  coefiicient.^°  The  efficient  evaluation  of  the  9j  symbol  is  critical  to  the 
ultimate  computational  success  of  the  present  program.  However,  present  and  projected 
computer  technology  offers  great  encouragement  in  this  regard.  Even  now,  textbooks  have 
reasonably  adequate  computer  codes  included  in  the  appendices.®^ 

The  evaluation  of  the  J  coefificient  proceeds  by  using  Eqs.  (56)  and  (57)  in  Eq.  (55) 
and  then  using  the  orthogonality  condition  expressed  by  Eq.  (45),  to  give 


(n/m) 


Hi +712— 1 


7l3  =  |”l'"^2|  +  lj2  ^3— |h“^!)2 


X  Su 


TI3+I 


2'5n2n3  +  ^  C(200, 713/3, m  -  mi;  714/3, m  -  mi)^„ 

n4=|n3-2|  +  l,2 


2  714 


(59) 
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Then,  analyzing  the  hmits  of  the  integraJs  in  Eq.  (53)  resulting  from  Eq.  (59),  produces  a 
result  which  is  a  specfic  form  of  Eq.  (1)  of  the  present  paper  (and  analogous  to  Eq.  (5.11a) 
of  Ref.  (23)) 


a 


5/2 


OO  OO 

E  E  E 


Ji=0  =  +  1 


_  max 
nj 


^(nl)i(n/)2 


X 


E 


mi=— /i 


j(nlm) 

(712/2 


r)  KJC'(o:,A  ) 


where  the  limits  on  the  722  sum  axe  defined  by 


(60) 


=  max(l2  +  1,  \ni  —  n|  —  1);  =  ni  +  n  +  1  (61) 

Note  that  C  is  defined  by  Eq.  (58)  and  J  is  defined  by  Eq.  (59).  The  expression  in  Eq. 
(60)  is  very  similau’  to  the  expression  given  in  Eq.  (5.11a)  of  Ref.  (23).  The  summation 
structure  is  exactly  the  same  and  the  limits  of  the  summations  are  exactly  the  same.  Most 
of  the  comments  concerning  the  character  of  the  convergence  follow  also  with  the  exception 
that  Eq.  (5.11a)  of  Ref.  (23)  is  an  expansion  which  is  valid  in  Hilbert  space  L^(R^),  while 
Eq.  (61)  above,  is  an  expsmsion  which  is  valid  in  Sobolev  space  W2^^(R^),  a  proper  subset 
of  Hilbert  space.  Thus,  the  right  hand  side  of  Eq.  (61)  converges  pointwise  for  all  f  and 
R.  Also,  as  pointed  in  Ref.  (23),  “The  number  of  terms  in  the  summations  over  I2  and 
n2  is  completely  determined  by  the  fixed  indices  n  and  /.  It  is  given  by  (2/  +  l)(2n  +  2) 
for  arbitrary  order  /j  and  ni,  and,  therefore,  the  number  of  terms  in  the  series  does  not 
incre2ise  with  the  order  of  terms.”  Note  that  I  have  substituted  I  for  L  and  n  for  N  in  the 
quotation. 

Actually,  an  addition  theorem  for  CSs  can  be  obtained  from  the  addition  theorem 
for  the  orthonormalized  Slaters  (ONSs)  of  Ref.  (23)  by  using  the  terminating  expansion 
of  CSs  in  ONSs  which  can  be  developed  easily  from  the  expressions  given  in  ref.  (19). 
However,  that  expression  would  not  be  as  compact  as  Eq.  (60).  One  of  the  central  points 
of  the  present  derivation  is  to  emphasize  the  use  of  the  FDSHs  and  the  integration  over  the 
surface  of  the  four-dimensional  Fock  hypersphere.  The  utility  of  this  technique,  using  the 
work  of  Alper^®  as  a  guide  with  the  present  appropriate  extensions,  will  be  demonstrated 
in  the  next  section. 


IV.  MULTICENTER  INTEGRALS  OVER  COULOMB  STURMIANS 


The  work  of  Alper,^*  which  applied  the  Fock  hyperspherical  projection^®,  with  the  de¬ 
velopments  of  Shibuya  and  Wulfman,^’^  to  the  problem  of  performing  multicenter  integrals 
over  exponential  type  orbitals,  was  very  nearly  successful  in  terminating  in  finite  series  ex¬ 
pressions,  all  of  the  integrals  which  axe  required  in  the  LCAO  approach  to  abinitio  quantum 
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chemistry.  However,  two  problems  plagued  his  developments.  The  first  problem  was  that 
the  basis  set  for  which  he  developed  his  expressions,  namely  hydrogenic  energy  eigenstates, 
are  not  complete  without  inclusion  of  the  continuum.^^  Thus,  his  formulas  were  of  limited 
utility  because  they  only  applied  for  discrete  eigenstates.  The  second  problem  was  that  the 
allowed  screening  paxzuneters  were  restricted  so  as  to  allow  the  Fock  projection  only  onto 
the  same  hypersphere.  That  is,  the  screening  parameter  determines  the  hypersphere  onto 
the  surface  of  which  the  projection  is  being  made.  Thus  all  hydrogenic  orbitals  occuring 
in  a  particular  integral  had  to  have  the  same  screening  parameter.  In  order  to  accomplish 
this  with  an  expansion  of  finite  length,  the  nuclear  charge  of  the  hydrogenic  orbitals  had  to 
be  related  to  the  principal  quantum  number  by  Z  =  nko  where  n  is  the  principal  quantum 
number  and  fco  is  a  fixed  constant,  common  to  all  of  the  hydrogenic  orbitals  used  in  the 
expansion. 

The  present  analysis  relaxes  both  of  these  retrictions.  Firstly,  the  CSs  are  complete  in 
Sobolev  space  with  a  purely  discrete  spectrum.  Secondly,  as  will  be  demonstrated,  since 
the  translations  are  made  in  R^,  using  Eq.  (59),  to  the  final  coordinate  centers  before 
projecting  onto  the  Fock  hypersphere,  Eq.  (21)  may  be  used  to  transform  the  integrals, 
with  a  finite  length  expansion,  to  integrals  over  CSs  with  the  same  screening  parameters 
before  projecting  onto  the  same  Fock  hypersphere.  It  should  be  noted  that  the  initial 
projection  onto  the  Fock  hypersphere  to  determine  J  of  Eqs.  (56)  and  (59)  is  distinct  from 
the  second  Fock  projection  made  to  perform  the  final  multicenter  integrals. 
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A  Calculation  of  the  Magnetic  Moment  of  the  A^'*' 

Milton  Dean  Slaughter* 

Department  of  Physics,  University  of  New  Orleans,  New  Orleans,  LA  70148 


A  fully  relativistic,  gauge  invariant,  and  nonperturbative  calculation  of  the  A"*"*"  magnetic  mo¬ 
ment,  fi^++  ,  is  made  using  equal-time  commutation  relations  (ETCRs)  and  the  dynamical  concepts 
of  asymptotic  SUf{2)  flavor  symmetry  and  realization.  Physical  masses  of  the  A  and  nucleon  are 
used  in  this  broken  symmetry  calculation.  It  is  found  that  =  2.04/ip,  where  Pp  is  the  proton 

magnetic  moment.  This  result  is  very  similar  to  that  obtained  by  using  static  SUf{Q)  symmetry 
or  the  broken  SUf{6)  quark  model. 


I.  INTRODUCTION 

Why  should  one  study  the  A++  magnetic  moment? 

The  answer  is  that  such  a  study  can  yield  information  on  hadronic  structure.  This 
information  is  important  to  the  development  of  a  correct  quantitative  description  of  hadrons. 
By  comparing  experimental  and  theoretical  determinations  of  the  A++  magnetic  moment, 
one  can  obtain  constraints  on  QCD,  quark  models,  Bag  models,  etc.. 

In  this  talk  results  are  presented  which  are  applicable  to  pion-proton  bremsstralung 
experiments  at  energies  near  the  A  resonance. 


II.  ASYMPTOTIC  LEVEL  REALIZATION 

We  first  consider  the  well-known  chiral  517(2)  ®  517(2)  charge  algebra 

[A„^,A„-]=2V^.  (1) 

Sandwich  the  above  commutator  between  the  identical  baryon  SUf{N)  states  ^B(a,  A| 

oo),  helicity  A,  and  physical 


and 

SUf 


N)  index  a.  Then 


_ 

/ 

p  -*  oo, 

p 
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^  {(B(a,  A)|  A^+  Iiil)  (ulI  A^r-  |B(a,  A))  -  (B(a,  A)1  A,-  Iiil)  (hlI  Afl.+  lB(a,  A))}  -  C(a,  A), 

(2) 

where  ^(a,  A)  =  2  (B(a,  A)|  V3  |B(a,  A))  is  a  pure  number.  Among  the  sum  over  single 
particle  hadronic  intermediate  states,  we  distinguish  (for  given  a  and  A)  the  fractional  con¬ 
tribution  /^(a,  A)  coming  from  all  the  states  ni  belonging  to  a  level  L{L  =0,1,2,...),  where 
L  connotes  angular  momentum.  Note  that  upon  varying  the  SUf{N)  index  a  which  pro¬ 
duces  non- vanishing  values  of  C(a,  A),  the  intermediate  states  ni  also  undergo  corresponding 
SUf{N)  rotations. 

The  algebraic  level  realization  hypothesis  then  states  that  /^(a.  A)  for  a  given  level  L 
will  depend  on  A  but  not  on  a. 

Thus,  the  fractional  contribution  from  each  level  L  is  universal  and  invariant  under 
SUf{N)  rotation. 

We  next  define  relevant  axial-vector  matrix  elements  as  follows: 

(p,l/2lA,+  ln,l/2)  =  /  , 

(A+M/2I/!,.  IAM/2)  s  -yij  ,  (3) 

(A++  l/2M,+  |y,l/2)  =  -VS  h  . 

Using  the  commutator  Eq.  (1)  and  asymptotic  level  realization  for  the  L  =  0  ground 
state  baryons  =  i''’,  §'''),  it  can  be  shown  easily  that 

III.  COMMUTATORS  USED  TO  CALCULATE  THE  A++  MAG¬ 
NETIC  MOMENT 

[[JLiO),A..],A,.]  =  2J^{0)  (5) 

[[73^(0),  A„+],A„-]  =  2  J3^(0)  (6) 

Here,  electromagnetic  current,  and  A„+  and  A„-  are  axial-vector  charge 

operators  (generators).  can  be  written  as  +  Uj*,  where  is  the 

isovector  part  of  the  electromagnetic  current  and  V5  is  the  isoscalar  part. 


35 


Insert  the  algebra,  Eq.  (6),  between  the  ground  states,  A  =  — 1/2|  and 

B'(/3,7,A  =  +1/2^  with  's*'  oo  and  T  -*•  oo.  Then  the  following  10  equations  are 
generated: 

(A++I  \JS,  /I..]  A,-  |A++)  -  (A++I  A,-  \JS,  A,*)  1a++)  =  2  (A++|  J?  |a++)  (7) 

(a+I  [js,  a,*]  a,-  |a+)  -  (a+I  a,-  [js,  a,*  ]  |a*)  =  2  (a+|  7?  |a+)  (8) 

(a»|  [js,  a,*]  a,-  Ia”)  -  (a“|  a.-  \js,  a,*  ]  |a“)  =  2  (a”|  |a")  (9) 

(A’l  173“,  A,*l  A,-  |a-)  -  (A-j  A,-  [73»,  A.*  ]  jA')  =  2  (a-|  73'  |a-)  (10) 

(p|  173»,  a,*|  a,-  Ip)  -  (pI  a,-  [73",  A,*  1  Ip)  =  2  (p|  j;  |p)  (ii) 

{n|  [73-,  A„]  A,-  |n)  -  {n|  A,-  (73',  A,+  ]  |n>  =  2  (n|  JS  In)  (12) 

(a+I  [73»,  A„|  A,-  Ip)  -  (a+I  a,-  [73',  A.*  1  Ip)  =  2  (A+j  73“  |p)  (is) 

(pi  (73“.  A,*l  A,-  1a+)  -  (p|  A,-  [73“,  A,*  1  |A+)  =  2  (p|  73'  |A^-)  (14) 

(A"!  [73“,  A,+l  A,-  |n)  -  (A“|  a,-  [73',  A,*  ]  |n)  =  2  (a»|  JS  |n)  (15) 

(n|  |73“,  A„1  a,-  1a“)  -  (n|  A,-  [JS,  A.,  ]  jA”)  =  2  (n|  73“  |a")  (16) 

Define: 

<p,-l/2,5'|  J3lp,l/2,f>  =  b  =  -  <  n,-l/2,s|J3|n,l/2,r>, 

<  n,-l/2,i'|  J3I  A°,l/2,r>  =  c  =  +  <p,-1/2,s|J3|  A+,l/2,f>, 

<  A°,-1/2,J'|  J3ln,l/2,r>  =  d  =  +<  A+, -1/2, s\Jz\p, 1/2, t>  (17) 

<  A+, -1/2,5*1731  A+,l/2,r>  =  a  =  -  <  A°,-l/2,s  |  J3I  A°,l/2,r>  , 

<  A++, -1/2,5  IJ3I  A++,  1/2,  r>  =  3a  =  -  <  A", -1/2, 5  1 J3I  A",  1/2,  f  >  . 

Thus,  with  the  definitions  given  in  Eqs.  (17),  Eq.(7),  for  instance,  becomes  explicitly: 
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+  ISah^  +  Zdgh^  —  Zcgh  —  6bh^  =  6a  (18) 

One  can  solve  this  system  of  10  simultaneous  non-linear  algebraic  equations  and  finds 
that  (physical  solution): 

a:6:c:d=l:-2:  ±2v^/5  :qF2V2/5,  (19) 

corresponding  to  A  =  ±|  /. 

(Note  that  the  Lorentz  index  n  ha^  been  dropped  for  convenience). 

Clearly  then,  we  have  the  result: 


<  A+,-l/2,J'|7|‘|  A+,l/2,f>  =  -(l/2)<p,-l/2,J‘|J3np,l/2,r>,  sj^oo.  (20) 

What  about  the  matrix  element  <  A''',  — 1/2,  s|T^|  A''',l/2,f  >  ? 

To  obtain  it,  we  first  write  all  matrix  elements  in  terms  of  their  isovector  and  isoscalar 
parts,  i.e. 


{p\  Ip) 

i  -{■  S, 

(n|  Jtm  |n) 

= 

—i  -h  6, 

(A++|J„„|A++) 

= 

3/  +  0, 

(A+|J.,„|A+) 

= 

i  + 

(A°|J„„|A0) 

= 

-j  + 

(AiJ.,„|A-) 

= 

—Zj 

(n|Jem|A°) 

= 

(Pl  1 A+) 

(A°|Je,,|n) 

= 

(A+|Jem  Ip) 

In  Eq.(21)  i,  j,  k  and  /  denote  the  isovector  matrix  elements  whereas  /3  and  S  denote  the 
isoscalar  matrix  elements. 

Next  insert  the  algebra,  Eq.  (6),  between  the  ground  states,  =  —1/2  and 

B'(^,  t'  ,X  =  -1-1/2^  with  s  -♦  oo  and  t 
We  find  that  /3  =  a  =  0,  i.e. 


oo. 


<  A,-1/2,J*|F/|  A,l/2,f>=0  (22) 

for  all  charge  states  of  the  A.  Therefore,  we  have  as  our  main  result: 

<  A+,-l/2,i'|J,';„|  A+,l/2,r>=-(l/2)<p,-l/2,s|J3np,l/2,r>,  s,f^oo.  (23) 

We  explore  the  consequences  of  Eq.  (22): 
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First: 


(p(i’,A)|J»„(0)|p(r,A') 

FM 

Fi(0)  =  1 


F^{q^)r  +  ,  (24) 

^2(0)  =  Kp  (proton  anomalous  magnetic  moment), 

(fi(0)  +  ■^2(0))(2^)  =  proton  magnetic  moment 


(A(s,A)|J.^,„(0)lA(t:A')) 

r-, 

F{{0)  =  Qc 

fiA 


Here,  q  =  yf  —  p. 


y^5i(j-,A)  (Ol  “i(O'), 

F2*(0)  =  ka  (a  anomalous  magnetic  moment), 

(F’i-(0)  +  f’2*(0))(^)(2^)  =  ^  magnetic  moment  . 

(25) 


IV.  Kinematics 

We  take  t  to  lie  along  the  z-axis  and  in  the  x  —  z  plane  such  that 


ji]  =  00  ,  =  rt^  00  ,  0  <  r  <  1  ,  r  held  constant. 

We  find  that 


21  (1“^^)  2  j  21 

q^\p^p+-,  =  — - — -m^  -  ,  and  r|A-A+7 


(1 

r 


m 


*2 


Note  that  in  the  limit  where  i — >  1,  that  -»  0  independently  of  m  or  m*  as  Sx  0. 
Now  choose  p  =  0,  evaluate  Eq.  (23),  and  use  the  fact  that  (0)  =  \  {kp-  K)  ■  We 
then  find  that: 


From  Eqs.  (21)  and  (22),  we  also  find  that: 

F2^+-^(0)  =  3^2^^  (0) 

=  {2  +  F,‘^“(o)}{5,}[,^],  (27) 

=  S^iK-  K)  ■ 
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Numerically,  we  obtain  (when  m*  =  1.232  GeV/c): 


/^P 


Theory 


—  2.04  ,  iTheory  “  2m ' 


(28) 


OTHER  THEORETICAL  CALCULATIONS 


/iA++  [e/2m] 

SC'(6)[ll 

5.58 

SU{6)  Mass  Corrected[2] 

4.25 

Bag  Model  Corrections  to  Quark  Model[3] 

4.41  -  4.89 

Lattice  Gauge[4] 

4.91  ±0.61 

QCD  Sum  Rule [5] 

2.7  ±0.4 

“EXPERIMENTAL”  VALUES 


[e/27Tij 


UCLA[6] 

3.7 -4.2 

SIN[7] 

4.6  —  4.9 

Musakhanov[8] 

3.6  ±2.0 

Pascual  and  Tarrach[9] 

5.6  ±2.1 

Heller  et  a/.[10] 

7.0  -  9.8 

Wittman[ll] 

5.58  -  7.53 

Bossard  et  a/. [12] 

4.58  ±0.33 

V.  CONCLUSIONS 

The  A++  magnetic  moment,  )Ua++)  can  be  calculated  in  a  fully  relativistic,  gauge  invariant 
manner  without  assuming  517(6)  symmetry. 

Estimated  theoretical  error  is  as  7%  since  radial  configuration  mixing  of  excited  A  states 
is  neglected. 

Results  are  in  good  agreement  with  experiment. 
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Abstract 


Spectra  of  He-like  Krypton  from  TFTR  Plasmas;  a  potential  Tj  diagnostic  for 
ITER.*  A.  J.  Smith+,  M.  Bitter,  H.  Hsuan,  P.  Beiersdorfer+‘'',  B  Fraenkel+‘*"^,  K. 
W.  Hill,  K.  M.  Young,  and  M.  Zamstorff,  Princeton  University,  Plasma  Physics 
Laboratory,  Princeton,  NJ  08543.  We  have  obtained  spectra  of  heliumlike  Krypton 
from  TFTR  plasmas.  Since  the  wavelength  of  the  Ka  line  of  He-like  krypton  is  .95 
A,  the  spectra  were  observed  in  second  order  Bragg  reflection  so  as  to  utilize  Bragg 
angles  0>4OO  and  thus  having  sufficient  resolution.  It  was  also  necessary  to  use  Xe 
as  a  detector  gas  to  enhance  the  efficiency  of  the  position  sensitive  detectors  for 
detection  of  the  13  keV  photons.  These  measurements  are  important  because 
Doppler  broadening  measurements  on  the  KrXXXV  spectrum  have  been  proposed 
for  Tj  diagnostic  on  the  International  Thermonuclear  Experimental  Reactor,  In  this 
study  we  also  look  at  the  effect  of  highly  ionized  krypton  on  plasma  performance. 


*Supported  by  DoE  Contract  No.  DE-AC02-76-CHO-3073 
■^Lock  Haven  University,  Lock  Haven  PA  17745; 
■'■■'■Lawrence  Livermore  National  Laboratory,  Ca  94550; 
■^■'■■'■Hebrew  University  Jerusalem,  Israel. 
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INTRODUCTION 


The  conceptual  design  of  the  International  Thermonuclear  Experimental 
Reactor  (ITER)l  indicates  that  this  next  generation  tokamak  will  have  a  stored 
energy  of  600  MJ,  a  plasma  volume  of  550  m^,  and  pulse  lengths  of  400-2500  sec. 
In  comparison  present  generation  tokamaks  such  as  TFTR  (Tokamak  Fusion  Test 
Reactor)  have  the  stored  energy  of  1-5  MJ,  a  volume  of  50  m^  and  pulse  lengths  of 
5  sec.  Two  of  the  main  issues  of  concern  with  ITER  operation  will  be,  a)  controlled 
release  of  stored  energy,  and  b)  measurement  of  the  central  ion  temperature  Ti(0). 

It  has  been  suggested  that  bpth  of  these  problems  be  solved  by  the  injection  of  high- 
Z  elements  into  the  plasma.  In  this  regard  krypton  (Z=36)  appears  to  be  the  natural 
choice;  it  is  an  inert  gas  and  therefore  easy  to  inject  into  the  plasma,  and  heliumlike 
krypton,  KrXXXV  is  the  dominant  ionization  state  at  ITER  operating  temperatures 
~7-30  keV.  Doppler  broadening  measurements  of  the  resonance  line  of  He-^like  Kr 
could  be  used  as  a  Ti(0)  diagnostic,  while  the  lower  charge  states  could  release 
through  the  emission  of  line  radiation  from  the  plasma  edges. 

It  is  therefore  important  to  study  the  effects  of  Kr  on  plasmas  and  to 
investigate  the  conditions  under  wliich  heliumlike  Kr  is  produced.  In  this  paper  we 
present  results  for  the  injection  of  krypton  into  ohmically  heated  TFTR  plasmas,  and 
also  spectra  of  heliumlike  krypton  obtained  from  tokamak  discharges. 

EXPERIMENTAL 

An  important  aspect  of  this  study  has  been  the  construction  of  a  high 
resolution  crystal  spectrometer  suitable  for  the  observation  of  the  satellite  spectmm 
of  KrXXXV.  The  wavelength  of  the  resonance  Ka  line  of 

heliumlike  Kr  is  0.95  A,  which  is  less  than  the  2d  spacing  of  natural  crystals,  the  Kr 
spectrum  must  therefore  be  observed  in  second-order  Bragg  reflection  at  the 
expense  of  a  significant  reduction  in  crystal  reflectivity.  The  13  keV  x-ray  photons 
of  He-like  krypton  were  recorded  with  one  arm  of  the  TFTR  Vertical  (high 
resolution)  spectrometer^  using  a  2023  -  quartz  crystal  ( 2d  =  2.7497  A;  with  a 
radius  of  curvature  1 1 .43  m)  in  second-order  Bragg  reflection.  This  arrangement 
leads  to  a  large  Bragg  angle  of  43.7 o  and  a  therefore  a  high  resolving  power  X/AX  = 
12,000.  The  spectrometer  also  reflects  6.5  keV  photons  in  first-order  Bragg 
reflection,  i.e.  with  much  higher  efficiency.  It  was  thus  necessary  to  include  an  8 
mil  A1  foil  absorber  in  front  of  the  detector;  this  has  a  1.5%  transmission  at  6.5  keV 
and  a  48%  transmission  at  13  keV  and  hence  a  significant  improvement  in  the 
signal-to-noise  ratio.  Furthermore,  xenon  was  used  as  detector  gas  in  the  position 


sensitive  multi-wire  proportional  counter,  this  increases  the  detection  efficiency  for 
13  keV  photons  by  40%. 

RESULTS  AND  DISCUSSION 

Krypton  was  injected  at  a  rate  of  0.2-0.35  Torr-liter/sec  for  0.2-0. 3  sec  into 
ohmically  heated  helium  discharges  at  three  different  plasma  currents  of  1.4,  1.6  and 
1 .8  MA.  These  discharges  had  electron  temperatures  in  the  range  5-6  keV  and 
electron  densities  1.5  -  2.5xl0l9  ni-3.  Fig.  1  shows  the  radial  electron  temperature 
profiles  for  these  discharges;  peaked  for  1 .4  and  1 .6  MA  discharges,  wide  for  the 
1.8  MA  discharge. 

A  survey  spectrum  of  the  emitted  line  radiation  for  the  1 .8  MA  discharge, 
obtained  by  the  TFTR  X-ray  Pulse  Height  Analysis  (PHA)  system^,  is  shown  in  Fig. 
2a.  The  spectrum  shows  the  Ka-line  from  Kr  and  various  metal  impurities  Cr,  Fe 
and  Ni,  superimposed  on  a  bremsstrahlung  continuum.  Since  the  energy  resolution 
of  the  PHA  system  is  230  eV,  the  peak  at  13  keV  includes  radiation  from  all  the 
krypton  charge  states  present. 

Figure  2b  shows  the  time  history  of  the  Kr  peak  at  13  keV  for  the  three 
discharge  currents.  The  intensity  of  the  Ka  peak  represents  a  chord  integrated 
measurement  along  the  horizontal  sight-line  of  the  PHA  system  through  the  plasma 
center;  it  is  expected  to  increase  with  the  electron  temperature  profile.  Further 
analysis  shows  that  at  the  peak  in  Fig.  2b  the  krypton  density  is  0.1%  of  the  total 
electron  density. 

The  total  radiated  power  (see  Fig.  3)  as  measured  by  the  TFTR  bolometer 
array^  increased  from  40%  to  100%  of  the  ohmic  power  input,  for  the  1 .4  and  1 .6 
MA  discharges  and  from  40  %  to  70%  of  the  ohmic  power  input  for  the  1.8  MA 
discharge.  The  fact  that  at  1 .8  MA  the  intensity  of  the  Ka  radiation  increased  while 
the  radiated  power  decreased,  indicates  that  at  broader  electron  temperature  profiles 
the  ionization  equilibrium  was  shifted  to  higher  charge  states  which  emit 
significantly  less  line  radiation  than  lower  charge  states. 

This  conclusion  is  confirmed  by  Fig.  4,  which  shows  high  resolution  spectra 
of  KrXXXV  obtained  for  plasma  currents  of  1.4  MA  (Fig.  4a)  and  1.8  MA  (Fig. 

4b).  The  data  were  accumulated  for  6  and  18  discharges  respectively.  For  a  peaked 
electron  temperature  profile.  Fig.  4a,  only  berylliumlike  and  boronlike  features  are 
present,  whereas  for  the  wide  electron  temperature  profile,  all  the  features,  including 
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heliumlike  lines  w,  x,  y,  and  z  and  their  associated  lithiumlike,  berylliumlike  and 
boronlike  satellites  (Fig  4b).  The  features  have  been  identilSed  using  the 
instrumental  dispersion  and  theoretical  wavelengths  from  Vainshtein  and 
Safronova^;  they  are  listed  in  Table  1 . 

A  noticeable  difference  between  the  KrXXXV  spectrum  and  high  resolution 
sateUite  spectra  of  TiXXI,  CrXXIII ,  FeXXVand  NiXXVII  is  that  the  heliumlike 
features  w,  x,  y,  and  z,  are  well  separated  from  each  other,  and  that  the  lithiumlike 
and  beryllium  features  are  now  much  closer  to  the  resonance  line  w.  An  important 
aspect  of  the  results  is  that  the  forbidden  line  z  seems  to  be  isolated  from  lithiumlike 
and  boronlike  features,  its  intensity  is  comparable  to  that  of  the  line  w,  but  unlike  w, 
z  is  uncontaminated  by  n  >  3  satellites.  It  would  therefore  appear  that  the  z  line  is  a 
better  candidate  for  Doppler  broadening  measurements,  and  preliminary  analysis 
gives  Tf  =  2.5  ±  0.4  keV,  in  agreement  with  Ti(0)  =  2.3  keV  from  FeXXV  which 
was  simultaneously  observed  on  another  arm  of  the  TFTR  vertical  spectrometer.  In 
contrast  the  measurement  from  the  KrXXXV  resonance  line  yielded  the  higher  value 
Ti(0)  =  6.3  ±  0.5  keV. 

CONCLUSIONS 

We  have  injected  krypton  into  ohmically  heated  TFTR  plasmas  in  order  to 
study  the  emitted  line  x-ray  radiation,  both  for  the  purpose  of  measurement  of  the 
central  ion  temperature  and  the  controlled  release  of  plasma  energy  for  future 
tokamaks  such  as  ITER.  We  have  observed  that  the  ionization  equilibrium  and  the 
loss  of  power  via  radiation  depend  on  the  electron  temperature  profile.  High 
resolution  spectra  have  been  obtained  showing  heliumlike  KrXXXV  lines  and  the 
associated  lithiumlike,  berylliumlike  and  boronlike  satellites.  Preliminary 
measurements  of  Tj  appear  to  be  in  agreement  with  established  methods.  Further 
studies,  both  theoretical  and  experimental  atomic  physics,  and  with  large  tokamaks 
are  needed  in  order  to  better  understand  the  effects  that  krypton  injection  will  have 
on  ITER. 
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TABLE  I:  Theoretical  wavelengths  for  the  transitions  in  heliumlike  krypton,  KrXXXV, 
and  the  main  n  =  2  satellites  from  Vainshtein  and  Safronova  (Ref.  14) 


Key 

Transition 

Wavelength 

(A) 

w 

ls2  ISo  -  ls2p  IPi 

0.94538 

X 

ls2  ISo  -  ls2p  3P2 

0.94708 

s 

1s22s  2$  1/2  -  Is2p2s(3p)2p2y2 

0.94746 

t 

1s22s  2$  1/2  -  Is2p2s(3p)2pj^ 

0.94804 

q 

1s22s2Si/2  -  ls2p2s(lp)2p3/2 

0.94961 

k 

ls22p  2pi/2  -  Is2p2  2D3/2 

0.94995 

j 

ls22p  2p3/2  -  Is2p2  2D5/2 

0.95137 

y 

ls2  ISo  -  ls2p  3pi 

0.95156 

r 

1s22s  2Si/2  -  Is2p2s  (Ip)  2py2 

0.95288 

P 

1s22s2  ^Sq  -  Is2p2s2  Ipi 

0.9529* 

z 

ls2  ISo  -  ls2s  3Si 

0.95525 

e 

ls22p  2p3/2  -  Is2p2  ^P5/2 

0.95615 

u 

1s22s  2$  1/2  -  ls2s2p(3p)4p3/2 

0.95652 

*from  present  calculations 
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Figure  Captions 


Fig.  1  Radial  profiles  of  the  electron  temperature  of  ohmically  heated  TFTR  plasmas  with 
different  plasma  currents:  (a)  Ip  =  1.4  MA,  (b)  Ip  =  1.6  MA  and  (c)  Ip  =  1.8  MA  .  The 
data  were  obtained  from  Laser  Thomson  scattering  measurements. 

Fig.  2  Results  from  the  TFTR  X-ray  Pulse  Height  Analysis  system:  (a)  X-ray  spectrum 
consisting  of  the  bremsstrahlung  continuum  and  Ka  line  radiation  from  Cr,  Fe,  Ni  and 
Kr,  (b)  time  evolution  of  the  krypton  Ka  peak,  where  t  =  0  corresponds  to  the  starting 
time  of  the  krypton.injection. 

Fig.  3  Total  radiated  power  as  measured  by  a  bolometer  array  for  (a)  Ip  =  1.4  MA,  (b)  Ip 
=  1.6  MA  amd  (c)  Ip  =  1.8  MA.  The  dotted  lines  represent  data  from  comparison 
discharges  without  krypton  injection. 

Fig.  4  Krypton  Ka  spectra  as  measured  from  one  channel  of  the  TFTR  Vertical  Crystal 
Spectrometer  with  a  vertical  sightline  through  the  center  of  the  vacuum  vessel  at  Rq  = 
2.65  m.  The  data  were  accumulated  (a)  for  6  nearly  identical  ohmic  discharges  with  Ip 
=  1.4  MA,  and  (b)  for  18  nearly  identical  ohmic  discharges  with  Ip  =  1.8  MA. 
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Progress  in  Multicenter  Molecular  Integrals  Over  Slater-Type  Orbitals 

Herbert  W.  Jones 
Department  of  Physics 
Florida  A&M  University 
Tallahasee,  FL  32307 


Introduction 

Quantum  Mechanics  was  first  applied  to  the  hydrogen  molecule  in  1927  by 
Heitler  and  London  [1].  They  used  a  trial  solution  consisting  of  a  linear 
combination  of  atomic  orbitals  (LCAO)  of  the  hydrogen  atom  type  which  were 
solutions  to  the  Schrodinger  wave  equation  and  applied  the  variational  principle. 
This  method  evolved  and  was  generalized  to  the  use  of  combinations  of  Slater-type 
orbitals  as  basis  sets  and  applications  to  more  complicated  molecules.  At  most  we 
must  deal  with  four-center  and  two-electron  integrals.  During  the  1950  s  and 
1960’s,  the  University  of  Chicago  under  the  leadership  of  Mulliken  and  Roothaan 
systematized  the  molecular  orbital  method  and  developed  formulas  [2]  for  molecular 
integrals  over  Slater-type  orbitals  (STOs).  Difficulties  with  the  formulas  were 
encounted  .when  they  were  programmed  on  computers  with  finite  word  lengths. 
More  formidable  difficulties  developed  when  calculating  integrals  needed  for  three 
and  four-atom  molecules.  These  overwhelming  difficulties  gave  STOs  the 
reputation  of  being  "intractable". 

This  "bottleneck"  in  molecular  integrals  was  broken  by  Boys  [3]  with  the 
introduction  of  Gaussian-type  orbitals  (GTOs),  characterized  by  the  functional  form 
exp(-x2)  rather  than  exp(-x)  as  characterized  by  STOs.  Mathematically,  any 
function  can  be  represented  by  a  complete  set  of  STOs  or  GTOs.  But  in  the  case  of 
molecules,  generally  speaking,  STOs  are  three  times  more  efficient  than  GTOs 
because  they  represent  the  basic  physical  characteristics  of  electron  clouds  near 
nuclei  and  at  large  distances  from  nuclei.  And  in  specialized  situations,  such  as 
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NMR,  Mossbauer  effect,  and  collisions  of  particles,  STOs  are  very  much  more 
efficient  than  GTOs, 

In  1976  Sharma  [4]  showed  that  Lowdin  a-ffinctions  cotdd  be  represented  by 
a  finite  number  of  terms  of  exponentials  and  powers.  At  Florida  A&M  University, 
Jones  and  Weatherford  [5]  used  this  knowledge  to  reformulate  this  development 
using  computer  algebra  to  generate  a  "C  matrix"  with  integer  elements  for 
expressing  the  a-function.  The  basic  ideal  of  expanding  a  displaced  STO  in 
spherical  harmonics  with  functional  coefficients  originated  with  Coolidge  [6],  but 
was  formalized  and  used  extensively  by  Lowdin  [7].  We  were  able  to  demonstrate 
the  power  of  a  computer  algebra  approach  to  molecular  integrals  in  a  series  of 
papers  that  generated  formulas  for  two-center  overlap.  Coulomb,  hybrid,  exchange, 
and  three-center  and  four-center  [8]  integrals.  The  cancellation  problem  was  solved 
by  using  computer  algebra  that  could  produce  numbers  of  arbitrary  precision.  Our 
success  and  the  support  of  Dr.  William  A.  Lester,  Jr.,  made  it  possible  to  organized 
the  first  international  conference  on  molecular  integrals  in  Tallahassee  in  August, 
1981.  The  publication  [9]  of  the  proceedings  gave  a  boost  to  the  field. 

Progress  continues  to  be  made  on  the  multicenter  molecular  integral  problem. 
(Our  latest  innovation  is  in  the  use  of  Mathematica).  Computer  programs  continue 
to  be  developed  for  speed  in  Fortran.  Much  work  remains  to  be  done  and  the  issue 
of  a  competitive  program  with  GTOs  remains  in  doubt.  But  the  challenge  has  been 
taken  up  by  groups  in  Germany  [10],  Spain  [11],  Russia  [12],  France  [13],  and  of 
course,  at  FAMU. 

Lowdin  Alpha  -Function 

A  review  of  the  a-function  method  together  with  computer  algebra  generation 
of  formulas  has  been  given  before  [14].  In  essence,  a  displaced  STO  in  its  coordinate 
system  (R,0,  (}))  can  be  written  as  an  infinite  series  of  spherical  harmonics  with  a- 
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NMR,  Massbauer  effect,  and  collisions  of  particles,  STOs  are  very  much  more 
efficient  than  GTOs. 

In  1976  Sharma  [4]  showed  that  Lowdin  a-functions  could  be  represented  by 
a  finite  number  of  terms  of  exponentials  and  powers.  At  Florida  A&M  University, 
Jones  and  Weatherford  [5]  used  this  knowledge  to  reformulate  this  development 
using  computer  algebra  to  generate  a  "C  matrix'  with  integer  elements  for 
expressing  the  a-function.  The  basic  ideal  of  expanding  a  displaced  STO  in 
spherical  harmonics  with  functional  coefficients  originated  with  Coolidge  [6],  but 
was  formalized  and  used  extensively  by  Lowdin  [7].  We  were  able  to  demonstrate 
the  power  of  a  computer  algebra  approach  to  molecular  integrals  in  a  series  of 
papers  that  generated  formulas  for  two-center  overlap,  Coulomb,  hybrid,  exchange, 
and  three-center  and  four-center  [8]  integrals.  The  cancellation  problem  was  solved 
by  using  computer  algebra  that  could  produce  numbers  of  arbitrary  precision.  Our 
success  and  the  support  of  Dr.  William  A.  Lester,  Jr.,  made  it  possible  to  organized 
the  first  international  conference  on  molecular  integrals  in  Tallahassee  in  August, 
1981.  The  publication  [9]  of  the  proceedings  gave  a  boost  to  the  field. 

Progress  continues  to  be  made  on  the  multicenter  molecular  integral  problem. 
(Our  latest  innovation  is  in  the  use  of  Mathematical  Computer  programs  continue 
to  be  developed  for  speed  in  Fortran.  Much  work  remains  to  be  done  and  the  issue 
of  a  competitive  program  with  GTOs  remains  in  doubt.  But  the  challenge  has  been 
taken  up  by  groups  in  Germany  [10],  Spain  [11],  Russia  [12],  France  [13],  and  of 
course,  at  FAMU. 

Lowdin  Alpha  -Function 

A  review  of  the  a-function  method  together  with  computer  algebra  generation 
of  formulas  has  been  given  before  [14].  In  essence,  a  displaced  STO  in  its  coordinate 
system  (R,©,  <]))  can  be  written  as  an  infinite  series  of  spherical  harmonics  with  a- 
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function  coefficients  when  viewed  from  a  distance  a  in  the  coordinate  system  (r,  9, 
()))  displaced  in  the  z  -  direction: 


;t  =  AR^-'e-^®Y“(0,0) 


A 

■(2L-f-l)(L-i-M)!‘ 

K  - 

Y 

4;r(^-t-M)! 

^N-l 

4;r(L-M)! 

£-M 

_(2^-i-l)(^-M)!_ 

X  af^(^a,Cr)Yf 


where 


af^(Ca.Cr)  = 


X 


(2/+l)(^-M)!  nlm 

2(^  +  M)!  ife  U  ‘ 


(i.j) 


and 


e  e^'^  e'^''j,  r  <  a 

e’^'‘[(-l)'  e^®  -  e'^*].  r>  a 


A  (2C)N+  1/2  [(2N)!]-i/2  is  normalization  factor;  N,  L,  and  M  are  the  quantum 
numbers  of  the  orbital,  and  ^  is  the  screening  constant  or  orbital  exponent. 

It  leads  to  greater  efficiency  to  first  simplify  the  a-function  by  summing  over 
the  constant  values  to  obtain  T  and  X  one-dimensional  matrices: 


i=0 


x,(j)  - r'-'  T  cr''(i,j)[(-l)'e«-- e  '-](4:a)'-'''-'. 


r  >  a. 
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Hence, 


£T,(j)[(-l)-'e^"  -e-^'‘]r'-'-\  r  <  a 

j=0 

N+^ 
j=0 

Overlap  Integral 

As  an  example  of  our  methods,  the  overlap  integral  will  be  worked  out.  The 
simplest  non-trivial  example  is  the  two-center  overlap  integral,  i.e.,  the  integral  of 
the  product  of  two  orbitals  located  at  (0,0,0)  and  (0,0, a).  Thus 

S= 

If  we  take  Is  orbitals  (N=l,  L=0,  M=0,Ca=Cb=Q  the  C  matrix  is 

nloo  _  /I  i\ 

^0  -  \i  o) 

;(,=  Ar'‘''e«'Y“  («,♦)=  Y‘{e,<p). 

X^=  AR''-'e"‘Y"(e,«.)=  X  ci,y;(9.4,). 

S  =  A^^Jr^dr  e'^'’ jd0  sin  0  d</>  Y” Y°(0,^). 

The  orthogonal  conditions  for  spherical  harmonics  are 
jde  sin  9d<>  Y”(9,d)  Yp-  (0,d)=  s„  d„„. 

Hence,  only  one  term  (  ^=o)  survives  in  the  infinite  summation. 

S  =  A^  Jr^  dr  e'^''  a»,  with  A  =  2. 
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“•=4  i  i  Cr(io)  «■'■[(- ly  e''  -  e''ka)‘-‘(4:r)J-‘,r<  a 

“•  =  ?  i  i  CJ“(i,j)  «'■  -  e''‘teaf'(fry-‘,r>  a. 

^  1=0  j»0  ^ 

By  use  of  computer  algebra  one  may  arrive  at  the  well  known  formula 

S  =  e-^*[l  +  Ca  +  (Ca)V3]- 

/ 

On  the  other  hand,  the  integral  may  be  done  by  numerical  integration.  This 
approach  is  tinder  current  investigation. 

One  disturbing  aspect  of  the  C  matrix  is  the  explosive  growth  of  the  values  of 
the  C  matrix  elements  with  high  quantum  numbers.  But  this  is  just  a  reflection  of 
the  ever  present  danger  of  cancellation  errors,  i.e.,  the  subraction  of  nearly  identical 
numbers  with  a  resulting  loss  of  significant  figures.  We  have  succeeded  in  writing  a 
computer  algebra  program  to  expand  the  a-fimction  for  small  r  values  in  terms  of  an 
E  and  F  matrix  [15],  which  avoids  large  numbers.  A  definitive  and  simple  program 
for  C,E,  and  F  matrices  together  with  an  overlap  program  has  been  written  in  terms 
of  Mathematica  [16].  This  program  has  been  used  to  find  numerous  errors  in  Tai’s 
recent  article  on  overlap  integrals  [17].  Although  Mathematica  always  gives  correct 
answers,  it  is  far  too  slow  to  compete  against  FORTRAN  in  a  production  mode. 

Current  Research 

FAMU  has  gained  recognition  as  one  of  the  main  world  centers  for  the  study 
of  STO  molecular  integrals.  The  first  center  that  produces  a  complete  computer 
package  for  evaluating  STO  integrals  will  be  in  a  position  to  do  work  of 
unprecedented  accuracy  over  the  entire  range  of  molecular  phenomena.  Some 
groups  are  working  with  transforms  into  momentum  space,  but  we  continue  to 
make  progress  in  coordinate  space.  Jones  and  Etemadi  [18]  have  just  obtained  the 
most  accuracy  calcinations  to  date  for  the  ground  state  of  H2'‘’  using  LCAO  methods 
with  orbitals  up  to  principal  quantum  number  6  and  angular  momentum  number  5. 
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For  our  LCAO  we  used  a  very  systematic  sequence  of  orbitals: 

¥=t  I  Cm.  r"- e  '- Yf(0.0) 

N=1  L=0 

Note  that  only  one  screening  constant  was  used.  Surpisingly,  only  one  screening 
constant  was  needed  for  HeH2+.  Our  program  was  written  in  FORTRAN  and  run  on 
a  CRAY-YMP.  However,  it  was  checked  using  Mathematica.  The  strategy  of  using 
Mathematica  to  check  FORTRAN  is  the  key  to  our  current  research,.  In  order  to 
insure  accuracy,  we  are  installing  an  interface  program  so  that  a  program  running 
in  FORTRAN  may  occassional  obtain  crucial  evaluations  in  Mathematica.  This  is 
what  we  must  be  prepared  to  do  in  order  to  break  through,  if  we  meet  a  road  block. 

Two  of  our  students,  Andrew  Jackson  and  Jay  Jackson,  are  assisting  in 
implementing  a  semianalytic  approach  [20]  to  molecular  integrals.  All  -  analytic 
methods  can  be  slow  and  may  require  arbitrary  precision.  For  the  important  two- 
electron  integrals,  we  plan  to  evaluate  a-functions  along  a  grid  and  do  the  second 
integral  for  energy  by  using  Gauss-Legendre  and  Gauss-Laguerre  quadratures. 

Our  numerical  integrals  should  be  trouble-free,  since  the  ot-functions  and 
potential  do  not  have  singularities.  (It  is  singularities  and  function  pathologies  that 
plague  investigators  using  transform  methods). 

Our  mood  remains  optimistic;  we  believe  that  the  strategies  developed  here 
for  dealing  with  problems  that  "do  not  compute"  will  be  generally  useful. 
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A  JPL  demonstration  of  Deep-space  Optical  communications 
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Pasadena,  California  91109 


ABSTRACT 

The  successful  completion  of  the  Galileo  Optical  Experiment  (GOP^), 
represented  the  accomplishment  of  a  significant  milestone  in  JPL's  optical 
communication  plan.  The  experiment  demonstrated  the  first  transmission  of  a 
narrow  laser  beams  to  a  deep-space  vehicle.  Laser  pulses  were  beamed  to  the 
Galileo  spacecraft  by  Earth-based  transmitters  at  the  Table  Mountain  Facility 
(TMF),  California,  and  Star  fire  Optical  Range  (SOR),  New  Mexico.  The 
experiment  took  place  over  an  eight-day  period  (December  9  through 
December  16,  1992)  as  Galileo  receded  from  Earth  on  its  way  to  Jupiter,  and 
covered  ranges  from  1-6  million  km.  At  6  million  km  (15  times  the  Earth- 
Moon  distance),  the  laser  uplink  from  TMF  covered  the  longest  known  range 
for  laser  beam  transmission  and  detection.  This  demonstration  is  the  latest  in 
a  series  of  accomplishments  by  JPL  in  the  development  of  deep-space  optical 
communications  technology. 


T.  INTRODUCTION 

JPL’s  Deep  Space  Network  is  singular  in  its  ability  to  track  and  communicate  with  deep- 
space  probes.  As  we  continue  to  explore  and  better  understand  our  solar  system,  there  will  be 
demands  to  return  increasing  volumes  of  data  from  deep-space  probes  to  Earth.  Existing  70 
meter  DSN  systems  will  be  unable  to  support  this  increased  demand,  and  communication 
systems  at  optical  and  microwave  (32  GHz,  Ka  band)  frequencies  are  considered  the  most  viable 
emerging  technologies  for  the  missions  of  the  twenty-first  century.  With  NASA's  directive  of 
faster,  better  and  cheaper,  the  mini  and  micro  spacecraft  are  expected  to  play  an  ever 
increasing  role  in  the  future.  Weight  and  size  are  the  primary  driving  considerations  for 
subsystems  on  mini  and  micro  spacecraft  and  small  light  weight  optical  telescopes  coupled  to 
miniature  laser  transmitters  are  being  favorably  considered  as  an  attractive  alternative  to  the 
large  microwave  antennae. 

Although  there  has  been  significant  progress  in  the  component  technologies  for  deep- 
space  optical  communications,  it  was  becoming  apparent  that  a  systems-level  demonstration  to 
show  the  viability  of  optical  communications  was  needed.  The  Galileo  spacecraft's  second  flyby 
of  Earth,  part  of  the  Venus-Earth-Earth  Gravity  Assist  (VEEGA)  trajectory  [1],  afforded  a  unique 
opportunity  to  perform  a  deep-space  optical  uplink  with  the  spacecraft  as  it  receded  from  Earth 
on  its  way  to  Jupiter.  The  Galileo  Optical  Experiment  (GOPEX)  was  conducted  over  the  period 
December  9  through  December  16  from  transmitter  sites  at  Table  Mountain  Facility  (TMF), 
California,  see  figure  1,  and  at  the  Starfire  Optical  Range  (SOR),  New  Mexico;  see  figure  2.  The 
spacecraft's  Solid  State  Imaging  (SSI)  camera  was  used  as  the  optical  communications  uplink 
receiver.  The  experiment  had  three  principal  objectives,  namely; 

•  Demonstrate  laser  beam  transmission  to  a  spacecraft  at  deep-space  distances 


63 


•  Verify  laser-beam  pointing  strategies  applicable  to  an  optical  uplink  based  solely  on 
spacecraft  ephemeris  predicts 

•  Validate  the  models  developed  to  predict  the  performance  of  the  optical  link 

Galileo's  phase  angle  after  its  second  Earth  flyby  was  approximately  90°.  Thus  as  the 
spacecraft  receded  from  Earth,  it  looked  back  at  a  half-illuminated  Earth  image.  This  geometry 
allowed  laser  beam  transmission  against  a  dark-Earth  background,  and  the  experiment  was 
conducted  between  3:00  a.m.  and  6:00  a.m.  Pacific  Standard  Time.  The  nighttime  transmission 
had  Pvo  distinct  advantages: 

•  It  allowed  the  uplink  to  be  performed  at  the  frequency-doubled  Nd:YAG  laser 
wavelength  of  532  nm,  where  the  responsivity  of  the  solid-state  imaging  (SSI)  camera  is 
high. 

•  Long-exposure  camera  frames  could  be  taken.  This  facilitated  the  identification  of  the 
detected  laser  transmissions.  Analysis  of  the  stray-light  intensity  in  the  focal  plane  of 
the  camera  showed  that  the  camera  shutter  could  remain  open  for  up  to  800 
milliseconds  before  the  scattered  light  from  the  bright  Earth  saturated  the  pixels  that 
detected  the  laser  uplink. 

The  camera  was  scanned  across  the  Earth,  parallel  to  the  Earth's  terminator,  during  each 
exposure  to  facilitate  the  identification  of  the  laser  uplink  from  spurious  noise  counts  in  the 
camera  frame.  Using  this  strategy,  the  laser  uplink  appeared  as  a  series  of  evenly  spaced 
bright  dots  within  the  camera  frame,  and  was  quite  distinct  from  other  features  in  the  frame, 
^ser  uplink  data  were  received  on  each  of  the  seven  days  of  the  experiment  with  detections  on 
50  of  the  159  GOPEX  frames  taken.  The  demonstration  covered  a  period  of  eight  days,  but  other 
spacecraft  activities  precluded  laser  transmission  on  Day  5.  Because  of  an  unanticipated  bias  in 
the  scan  platform  pointing  ,  pulses  were  detected  on  only  two  of  the  frames  with  exposure 
times  less  than  400  milliseconds.  Inclement  weather,  aborted  transmissions,  and  restrictions 
imposed  by  regulatory  agencies  and  by  the  Project  Galileo  team  accounted  for  the  loss  of  data 
on  the  remaining  frames. 

The  two  GOPEX  laser/ telescope  transmitters  and  SSI  camera  receiver  are  discussed  in 
Section  II  and  III  of  this  paper.  Section  IV  describes  the  telescope  pointing  strategy  which  used 
bright  stars  in  the  vicinity  of  Galileo  as  references,  and  the  GOPEX  results  are  presented  in 
Section  V.  Conclusions  are  discussed  in  Section  VI. 


II.  GOPEX  Laser  Transmitters 

The  laser  transmitters  at  both  sites  consisted  of  a  frequency-doubled  Nd:YAG  laser 
operating  at  532  nm  coupled  to  an  optical  telescope  through  a  coude  mount  arrangement.  The 
transmitter  characteristics  are  given  in  Table  1. 

A  0.6-meter  equatorial-mount  astronomical  telescope  was  used  for  the  TMF  transmitter. 
This  is  the  same  telescope  that  was  used  in  1968  to  perform  the  laser  transmission  to  the 
Surveyor  7  spacecraft  on  the  Moon.  The  telescope  is  f/36  at  the  coude  focus,  and  the 
appropriate  beam-forming  lens  set  was  inserted  into  the  optical  train,  figure  3,  to  achieve  the 
required  laser  beam  divergence.  See  Table  1.  The  optical  train  was  designed  so  that  the 
transmitted  laser  beam  illuminated  a  15  cm  sub-aperture  of  the  primary  mirror.  The  principal 
benefit  of  this  technique  over  full  aperture  illumination  was  that  it  elirninated  the  large  loss  in 
transmitted  energy  that  would  have  been  caused  by  occultation  from  the  0.2-meter  secondary. 
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Table  1.  GOPEX  laser  transmitter  characteristics. 


Characteristic 


Table  Mountain 
Facility 


Starfire  Optical 
Range 


Wavelength,  nm 
Pulse  energy,  mj 
Repetition  rate,  Hz 
Pulse  width,  ns 
Beam  divergence, 
tirad 
Days  1-4 
Days  6-8 
Telescope  mirror 
diameter 
Primary,  m 
Secondary,  m 
Optical  train 
transmission 


532 

532 

250 

350 

15-30 

10 

12 

15 

no 

80 

60 

40 

0.6 

1.5 

0.2 

0.1 

60% 

43% 

s-ra  sa.  ,.r  ... 

0.6-meter  aperture. 

SOR’s  telescope  was  the  1.5-meter  system  that  is  used  aS^k'^^kd 

this  facility.  A  thin-film-plate  polarizer  served  stars  to’ be  observed  by 

microradian  beam  divergence,  respectively. 

The  SOR  laser  output  was  transmitted  through  the  full  1.5-meter  aperture  of  the 
occultation  by  the  secondary. 


TTT.  GOPEX  Receiver 

The  Galileo  SSI  camera  was  used  to  detect  the  GOPEX  uplink.  The  camera  is  mounted  on  the 
spacecraft  scan  platform  located  on  '“'f  focaT  length  f/8.5 

Sgmin°rsc*o?2  fS.^“heTgm^^^^ 

feeZs^‘ir"lsiTa"s^s?air*l'^^ 
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days,  and  the  more  sensitive  gain  state  3  with  160  detected  photoelectrons  per  dn  was  used  on 
subsequent  days. 

Field  correction  elements,  an  eight  position  filter  wheel,  and  a  two  blade  shutter  are 
positioned  along  the  optical  train  between  the  telescope  primary  and  the  focal  plane  of  the  SSI 
camera.  The  latter  two  elements  were  inherited  from  the  Voyager  program.  The  filter  wheel 
contained  one  clear  filter,  one  infrared  transmitting  filter  and  six  20  nm  bandpass  color- 
filters.  These  were  rotated  into  the  optical  train,  as  required,  to  enable  the  color  reconstruction 
of  an  imaged  scene.  The  green  color-filter  with  50  %  transmission  at  532  nm  (the  peak 
transmission  was  90%  af  560  nm)  was  used  for  GOPEX. 

The  SSI  shutter  is  operable  in  any  one  of  28  exposure  times  ranging  from  4.16 
milliseconds  to  51.2  seconds.  Exposure  time  selection  was  based  on  the  estimate  of  the  best 
balance  between  the  conflicting  requirements  of  short  duration  to  reduce  stray  light  effects, 
and  long  duration  to  ensure  that  enough  pulses  were  detected  to  confirm  the  laser  uplink. 
Improved  estimates  of  the  scattered  light  intensity  were  made  using  data  taken  at  the  Gaspra 
encounter.  Final  estimates  of  the  scattered  light  rates  due  to  Earth-shine  ranged  from  a  high 
of  110  eVmsec  on  the  end  of  the  first  day  for  the  SOR  location  to  32  e'/msec  on  day  8  for  both 
sites.  On  the  basis  of  these  estimates,  the  GOPEX  imaging  sequence  was  designed  with  shutter 
times  ranging  from  133  msec  to  800  msec.  Actual  scattered  light  rates  measured  during  GOPEX 
ranged  from  8  to  10  eVmsec.  The  low  actual  scatter  levels  would  have  allowed  longer  camera 
exposure  times  and  the  accumulation  of  more  data. 

IV.  GOPEX  telescope-Dointing  strategy 

Telescope-pointing  files  for  TMF  and  SOR  were  generated  from  updates  of  the  spacecraft 
ephemeris  files  that  were  provided  to  the  GOPEX  team  on  December  8  and  December,  11.  The 
strategy  was  to  off-point  the  telescope  from  reference  stars  located  within  0.5  of  the 
spacecraft  position.  Over  the  eight-day  period,  six  guide  stars  of  magnitudes  6  to  10  were  used  to 
point  the  TMF  telescope  at  Galileo. 

Transmission  to  Galileo  was  accomplished  by  using  a  “point  and  shoot”  approach.  Here, 
the  telescope  was  set  to  track  the  reference  star  in  the  intervals  between  the  three-second 
bursts  of  laser  transmissions.  Two  and  one-half  minutes  before  laser  transmission,  the 
reference  star  was  positioned  in  the  center  of  the  field  of  view  of  the  focal  plane  aperture  at 
coude  and  the  telescope  was  calibrated.  Ten  seconds  prior  to  transmission,  the  telescope  was 
pointed  to  Galileo’s  predicted  location  and  set  to  track  the  spacecraft  for  the  next  thirteen 
seconds.  This  procedure  was  repeated  during  the  three-minute  to  six-minute  intervals  between 
the  laser  transmissions.  Because  the  telescope  calibration  was  performed  just  before 
transmission,  the  pointing  errors  introduced  by  mount  sag  were  reduced  significantly.  In 
addition,  the  high  elevation  of  the  spacecraft  during  the  uplink — the  experiment  was  conducted 
when  the  spacecraft's  elevation  from  TMF  was  greater  than  30  — and  the  proximity  of  the 
reference  stars  to  Galileo's  position  obviated  the  need  to  implement  atmospheric  refraction 
compensation  techniques  while  pointing  to  the  spacecraft. 

The  beam  width  of  the  rf  beam  to  Galileo  was  approximately  forty  times  that  of  the  laser 
beam,  2  millirad  versus  60  microrad.  To  verify  the  accuracy  of  the  pointing  predicts  SOR 
dithered  the  laser  beam  in  a  85  microrad  radius  circle  about  the  spacecraft's  predicted  position 
while  the  TMF  transmitter  pointed  directly  to  Galileo’s  predicted  position.  This  strategy  was 
implemented  for  several  of  the  long-duration  frames  (frames  with  exposure  times  greater  than 
400  milliseconds)  on  the  first  day.  The  results  from  one  of  the  dithered  uplinks  are  shown  in 
Figure  4.  Nine  pulses  can  be  clearly  discerned  in  the  figure;  seven  are  from  the  15-hertz  TMF 
transmitter,  and  two  are  from  the  10-hertz  SOR  transmitter.  Without  beam  scanning,  a  total  of 
four  pulses  would  have  been  detected  from  the  SOR  transmitter.  The  presence  of  only  two 
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pulses  from  SOR  and  of  seven  from  TMF  clearly  demonstrates  that  the  error  in  the  telescope 
pointing  predicts  was  significantly  less  than  85  microradians.  This  was  further  confirmed  by 
the  successful  use  of  a  60-microradian  beam  from  TMF  for  laser  transmissions  on  the  last  three 
days  of  GOPEX. 


V.  GOPHX  RESULTS 

Table  2  gives  a  summary  of  the  detected  GOPEX  laser  transmissions  over  the  duration  of 
the  experiment. 


Table  2.  Summary  of  detected  laser  signals. 


Day 

Shutter  speed 
(ms) 

Frames  received 

Frames  with 
detections 

1 

133 

9  of  10 

0 

200 

24  of  25 

0 

400 

19  of  20 

6 

800 

5  of  5 

4 

2 

200 

5  of  5 

0 

267 

15  of  15 

0 

533 

15  of  15 

11 

800 

5  of  5 

5 

3 

200 

5  of  5 

0 

267 

10  of  10 

0 

533 

5  of  5 

5 

4 

200 

3  of  3 

lb 

267 

4  of  4 

0 

533 

3  of  3 

2b 

5 

No  activity 

planned 

6a 

133 

3  of  3 

0 

267 

6  of  6 

0 

533 

3  of  3 

3 

7a 

200 

3  of  3 

1 

400 

3  of  4 

3 

800 

3  of  3 

3 

8a 

267 

2  of  2 

0 

533 

4  of  4 

4 

800 

2  of  2 

2 

^  Adverse  weather  at  Starfire  Optical  Range  precluded  laser 

transmission  on  this  day. 

b  Adverse  weather  at  Table  Mountain  Observatory  precluded 

laser  transmission  on  this  day,  and  it  was  cloudy  at  Starfire  Optical 
Range. 
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Over  the  eight-day  period,  transmissions  to  the  spacecraft  were  made  over  a  range 
beginning  at  600,000  kilometers  on  the  morning  of  December  9  and  ending  at  6,000,000 
kilometers  on  the  morning  of  December  16.  Signals  were  successfully  detected  on  each  of  the 
experiment  days,  although  not  on  all  frames  within  a  given  day.  Unfavorable  weather  (which 
caused  outages),  regulatory  agency  restrictions  on  transmissions,  temporary  signal-to-noise 
anomalies  on  the  downlink,  and  an  unexpected  camera-pointing  bias  error  resulted  in  the  lack 
of  detection  on  several  frames.  Final  results  show  that  the  laser  uplink  was  successfully 
detected  on  50  camera  images  during  the  experiment  window.  Figures  4  and  5  show  two 
representative  images  of  the  detected  laser  pulses. 

Adverse  weather  at  the  sites  and  not  telescope  pointing  was  the  most  severe  impediment 
to  successful  detection  of  the  laser  transmissions.  Winter  storms  at  TMF  and  SOR  brought  snow, 
heavy  clouds,  and  ground  fog  to  these  facilities.  Transmission  from  TMF  was  most  affected  on 
the  first  and  fourth  days  of  the  experiment.  The  last  seven  frames  obtained  on  the  first  day 
were  taken  with  TMF  completely  overcast  and  SOR  in  daylight.  On  the  fourth  day,  falling  snow 
at  TMF  precluded  transmission  from  this  facility;  also  on  that  day,  during  only  one  of  the  ten 
transmissions  was  there  clear  sky  between  the  SOR  transmitter  and  the  spacecraft.  Falling 
snow  and  heavy  cloud  cover  prevented  transmission  from  SOR  on  the  last  three  days. 

Restrictions  from  regulatory  agencies  also  caused  data  outages.  Transmission  of  the 
GOPEX  laser  beam  into  space  required  the  concurrence  of  the  US.  Space  Defense  Operations 
Center  (SPADOC).  On  the  first  day,  SPADOC  restrictions  prevented  TMF  from  transmitting  during 
four  frames.  An  additional  frame  was  lost  because  the  ground  receiving  station  (at  Goldstone, 
California)  momentarily  lost  lock  on  the  Galileo  spacecraft  downlink  signal.  Owing  to  the  loss 
of  downlink  signal,  the  orientation  of  the  spacecraft  could  not  be  confirmed,  and  since  one  of 
the  GOPEX  concurrence  conditions  was  that  laser  uplink  would  proceed  only  if  the  spacecraft 
orientation  was  known,  no  laser  transmissions  were  sent  during  this  data  outage. 

During  the  first  two  days  of  GOPEX,  the  spacecraft  orientation  resulted  in  the  low-gain 
antenna  being  pointed  away  from  Earth.  This  resulted  in  a  low  signal-to-noise  ratio  (SNR)  of 
the  spacecraft  downlink  and  was  evidenced  by  the  numerous  burst  errors  in  the  data  files.  The 
GOPEX  images  for  these  days  showed  numerous  streaks  across  the  frames  and  made  it  difficult  to 
discern  successful  laser  transmissions  on  the  images.  On  the  second  day,  just  after  the  GOPEX 
uplink,  a  planned  spacecraft  maneuver  that  increased  the  SNR  of  the  radio  frequency 
downlink  was  executed.  This  resulted  in  clearer  GOPEX  images  for  the  remainder  of  the 
demonstration. 

The  GOPEX  demonstration  required  that  the  SSI  camera  be  operated  in  a  mode  for  which 
it  was  not  designed  (that  is,  slewing  the  camera  during  imaging).  To  get  the  GOP^  transmitter 
sites  in  the  field  of  view  during  the  slew,  the  camera  was  initially  pointed  to  a  position  above  or 
below  the  targeted  direction  and  the  shutter  was  opened  at  a  prescribed  time  after  the  start  of 
the  slew.  Uncertainties  in  the  stray-light  intensity  in  the  focal  plane  of  the  SSI  camera  dictated 
the  shutter  times  used  for  GOPEX.  The  times  chosen  ranged  from  133  milliseconds  to  800 
milliseconds,  and  these  were  loaded  into  the  spacecraft  sequence  of  events  prior  to  GOPEX.  As 
the  experiment  progressed,  it  was  observed  that  laser  transmissions  were  consistently  detected 
only  on  frames  with  greater  than  400-millisecond  exposure  times;  laser  pulses  were  detected  in 
one  200  msec  exposure  frame  taken  on  day  4  and  again  on  day  7.  The  consistent  absence  of 
detections  on  the  shorter  duration  frames  was  traced  to  a  pointing  error  caused  by  the  scan 
platform  acceleration  being  slower  than  predicted.  As  a  result,  no  clear  evidence  of  laser 
transmission  was  observed  on  88  of  the  90  frames  taken  with  exposure  times  less  than  400 
milliseconds. 
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After  determining  the  reasons  for  the  missing  pulse  detections  on  selected  frames,  the 
remaining  frames  were  analyzed  to  calculate  the  detected  pulse  energy  statistics  and  to 
compare  those  statistics  with  theoretical  predictions.  Figure  6  shows  a  typical  comparison.  The 
histogram  of  the  detected  pulse  energies  was  constructed  from  a  total  of  51  pulses  received  on 
this  day,  and  it  is  plotted  along  with  the  theoretical  distribution  using  atmospheric  turbulence 
data  and  the  appropriate  turbulence  model.  The  error  bars  in  the  histogram  heights  are  large 
because  of  the  small  data  set.  They  show  the  expected  range  of  the  relative  frequency  of 
detections  for  the  GOPEX  demonstration  performed  again  under  the  similar  atmospheric 
conditions.  The  poor  fit  near  the  lower  laser  energies  is  due  to  the  difficulty  in  estimating  low 
detected  laser  energies.  The  log-normal  intensity  distribution,  based  on  statistics  of  the 
measured  data,  is  also  shown.  The  data  show  good  agreement  between  the  measured  and 
theoretically  predicted  distributions. 


VI.  CONCLUSION 

GOPEX  represents  the  achievement  of  a  significant  milestone  in  deep-space  optical 
communications.  The  laser  transmission  was  performed  from  transmitters  located  at  TMF, 
California,  and  SOR,  New  Mexico.  The  laser  uplink  was  detected  on  every  day  of  the  experiment — 
out  to  a  range  of  6,000,000  kilometers  for  the  TMF  transmission.  The  camera  images  returned 
from  Galileo  and  the  analysis  of  the  data  show  that  all  the  experiment  objectives  were  achieved. 
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Figure  1. 


Laser  transmission  from  0.6-meter  telescope  at  TMF,  Wrightwood,  CA. 


Figure  2. 


Laser  Transmission  from  1.5-meter  telescope  at  SOR,  Albuquerque,  NM. 


PROFILE  MONITOR 


Figure  3.  GOPEX  optical  train,  at  TMF.  Relay  mirrors  (labeled  "A"  in  the  figure)  are 
appropriately  inserted  into  the  optical  train  to  obtain  the  required  beam 
divergence. 
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Figure  4.  Camera  image  showing  detected  pulses  from  TMF  and  SOR.  To  the  right  of  the 
figure  is  the  blurred  image  of  the  portion  of  the  Earth  illuminated  by 
sunlight.  Detections  of  the  TMF  laser  pulses  are  shown  on  the  left  and  of  SOR 
pulses  closer  to  the  right.  The  missing  pulses  are  due  to  the  spatial  scanning 
of  the  SOR  beam. 


Figure  5.  Camera  image  showing  laser  uplink  detections  from  TMF  (15-hertz  repetition 
rate)  and  SOR  (10-hertz  repetition  rate)  on  Day  2  of  GOPEX.  SOR  scanning  was 
off. 
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”111  International 


CarwiI  James 

A  speech  prepared  for  the  National  Society  of 
Black  Physicists  Annual  Meeting, 

Tallahassee,  Florida 

The  XXIV  International  Physics  Olympiad  will  he  hosted  this 
College  of  William  and  Mary  in  Williamsburg,  Virginia.  TheU.S.  is  hosting  this 
evlntontZheXoflts  suclessful  competition  oi  the  XXIII  Olympiad  in 
Finland  That  comvetition  was  both  a  proud  achievement  for  the  American 
^^konVf  nyTiT,  Teacher,  and  the  many  dedicated  individuals  involved 
with  the  United  States  Physics  Team,  and  a  rewarding,  though  challenging, 
experience  for  the  five  U.S.  participants. 

Last  Tulv,  I  had  the  privilege  of  representing,  along  with  four  other 
high  school  students,  the  United  States  of  America  at  the  XXIII 
Physics  Olympiad  in  Helsinki,  Finland.  The  contest,  which  took  place  from 
Tulv  5  to  July  13,  involved  178  participants  and  99  team  leaders  and  observers 
from  39  countries,  and  included  ten  hours  of  competition,  was  the  largest 
ever  International  Physics  Olympiad.  For  all  of  those  people  it  marked  the 
end  of  a  very  long  process  of  competition,  training,  and  other  preparations. 

The  International  Physics  Olympiad  is  so  named  bec^se  t  e 
International  Olympic  Committee  claims  exclusive  use  of  the  term  Olympics 
and  the  Olympic  Rings.  As  an  alternative,  the  U.S.  Physics  Olympiad  uses  a 
torch,  but  the  International  Organizers  don't  bother.  The  competition  itself  u 
a  test  of  individual  physics  skill.  While  teams  are  created  on  the  basis  o 
nationality,  the  teams  themselves  do  not  compete,  only  individuals  do.  O 
course,  it  isn't  hard  to  add  up  the  five  scores  of  each  team,  but  that  s  not  what 
the  competition  is  about.  Rather,  the  competition  is  the  meeting  o  physics 
students  from  around  the  world  in  both  a  competitive  and  a  social  setting. 
After  my  experiences  in  Helsinki,  I  truly  wished  for  the  opportunity  to  visi 
an  Olympic  village  and  see  if  the  athletes  developed  the  same  international 
friendships.  The  competition  part  of  the  experience  is  divided  into  a 
theoretical  competition  with  three  problems  and  an  experimenta 
competition  with  two.  Each  of  the  five  problems  carries  equal  weight. 

The  theoretical  competition  consists  of  three  questions  (almost  always 
multipart)  in  five  hours.  Very  few  people  finish  all  the  parts  and  a  perfect 
score  only  occurs  once  every  few  years.  Our  quesbons  dealt  with;  ^ 
ferris  wheel-type  satellite  in  low  earth  orbit,  satellites  as  blackbody  radiators 
interacting  with  the  sun,  and  the  nature  of  straight  chain  moleci^es  as  masses 
connected  by  ideal  springs  with  application  to  the  radiation  of  CO2  molecules. 
The  experimental  section  provides  much  less  information  to  the  participants. 
As  the  key  to  the  problem  is  devising  an  appropriate  procedure  tor 
measurement,  only  the  quantity  to  be  measured  and  detailed 
the  equipment  are  provided.  The  first  experiment  we  did  was  determinin,^ 
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the  proportionality  of  spark  gap  potential  to  distance  in  air  by  the  use  of 
piezoelectric  material  as  a  voltage  source.  The  second  experiment  was 
determining  a  variety  of  properties  of  some  filters  (including  estimates  of 
their  frequency  ranges)  and  the  width  of  CD  grooves,  using  a  penlight,  Duplo 
bricks,  and  a  protractor. 

The  test  is  composed  by  the  hosting  nation  and  translated  into  both 
English  and  Russian.  Team  leaders  then  conduct  debate  in  both  languages 
about  the  merit  of  the  questions,  seeking  to  gain  as  much  correspondence  as 
possible  between  their  team's  knowledge  and  the  topics  on  in  the  questions. 
After  this  process,  each  team  leader  translates  from  either  English  or  Russian 
to  his  native  language.  Once,  this  process  was  put  to  democratic  scrutiny, 
leading  to  hours  of  debate  on  translation.  This  has  been  replaced  by  the  honor 
system,  resulting  in  a  good  night's  sleep  for  both  the  contestants  and  the 
coaches. 

The  questions  in  the  Olympiad  can  deal  with  practically  anything.  In 
theory,  the  questions  should  involve  problems  from  a  specified  list  of  topics. 
Also,  Calculus  should  not  be  necessary  for  the  problems.  In  reality  the 
problems  usual  exceed  these  limits.  As  I  was  going  through  the  backlog  of 
problem  sets,  I  found,  as  warned,  that  they  became  more  and  more  difficult 
through  the  years.  Every  time  students  begin  to  approach  the  upper  limit  of 
scoring,  the  test  is  made  more  difficult  the  next  year.  After  23  years,  the 
difference  is  extraordinary. 

The  selection  of  U.S.  participants  in  the  International  Physics  Olympiad 
is  a  three  tiered  process,  each  tier  more  involved  than  the  last.  During  winter, 
the  American  Association  of  Physics  Teachers  (AAPT)  sends  out  letter  to 
numerous  students  and  teachers,  soliciting  teachers  to  recommend  their  top 
physics  students  to  take  the  preliminary  exam  and  the  students  to  try  and  take 
it.  I  received  such  a  letter  both  years  I  was  in  high  school,  undoubtedly  due 
my  participation  in  a  number  of  math  contests  in  prior  years.  Whatever  their 
source,  I  responded  by  taking  the  preliminary  examination  both  years  with 
the  blessing  of  my  physics  teacher,  Mr.  Robert  Shurtz.  My  first  year,  I  was  a 
student  in  AP  Physics  B,  clearly  unprepared  to  do  very  well  in  high  level 
competition.  Needless  to  say,  I  did  not  advance  beyond  this  first  phase  of 
competition,  which  was  a  combination  of  a  multiple  choice  test  and  a  series  of 
open  response  questions. 

In  my  junior  year,  I  was  one  of  four  students  to  take  the  exam  at  my 
school.  I  was  joined  by  my  friends  and  classmates,  Omar  Antar,  Mike 
Antonelli,  and  Charlie  Saulino.  Omar,  Mike  and  I  all  advanced  to  the  second 
stage  of  competition,  the  Semi-Final  exam.  While  Mike  opted  not  to  take  the 
second  exam,  Omar  and  I  went  on  to  the  Semi-Final  exam  with  the  aid  of 
crash  courses  from  Mr.  Shurtz.  This  exam  included  two  sections.  The  first  had 
four  challenging,  but  doable,  open  response  questions  in  one  hour.  The 
second  had  two  problems  which  we  couldn't  truly  expect  to  fully  complete 
within  the  allotted  two  hours.  Omar  and  I  both  took  the  test,  using  our  15- 
minute  intermission  between  the  sections  to  see  what  we  could  of  our  high 
school  talent  show,  and  waited.  Several  weeks  later,  1  found  out  I  had  been 
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selected  along  with  nineteen  other  students  to  be  on  the  1992  United  States 
Physics  Team  and  was  asked  to  attend  training  camp  during  my  last  week  of 
classes  in  high  school.  I,  of  course,  was  more  than  happy  to  oblige. 

When  I  arrived  at  Baltimore-Washington  International  Airport,  there 
were  no  representatives  of  the  United  States  Physics  Team  at  the  gate.  Finding 
this  surprising  and  unusual,  I  looked  at  my  table  of  arrivals  and  went  to  the 
gate  of  the  next  person  arriving.  There  still  were  no  USPT  representatives,  but 
when  the  flight  arrived  I  found  someone  else  on  the  team,  Jose  Lorenzo  who 
is  from  Miami.  It  was  not  even  five  minutes  before  we  ran  into  Chris,  who 
was  to  be  one  of  our  lab  instructors  at  training  camp.  The  people  meeting  us 
at  the  airport  had  been  waiting  for  us  at  the  entrance  to  the  terminal. 

From  the  time  we  checked  into  our  rooms  at  the  Adult  Education 
Center  at  the  University  of  Maryland,  College  Park  campus,  we  had  physics 
problems  to  work  on.  A  typical  day  include  a  two  hour  theoretical  test  with 
two  problems,  a  three  hour  lab,  a  three  hour  morning  lecture,  and  some 
discussion  of  problems  in  the  evening.  Both  the  laboratory  experiments  and 
the  theoretical  problems  got  more  and  more  difficult  as  the  week  went  on, 
serving  as  both  practice  for  us  and  selection  criteria  for  our  coaches.  We  did, 
once  in  a  while,  have  designated  free  time  at  training  camp,  but  mostly  we 
used  whatever  time  we  could  find.  We  were  taught  and  encouraged  and  led 
by  our  two  academic  directors.  Professor  Larry  Kirkpatrick  of  Montana  State 
University  and  Dr.  Avi  Hauser  of  Bell  Labs  as  well  as  Mr.  Ted  Vittitoe  of 
Libertyville  High  School  in  Illinois  and  our  academic  advisor,  Arthur 
Eisenkraft  who  was  at  the  time  working  primarily  on  the  1993  Olympiad  in 
Williamsburg.  Everyday,  a  researcher  came  to  speak  to  us  about  their  work. 
Their  speeches  discussed  things  ranging  from  scanning  tunneling  microscopy 
to  lasers  to  cosmic  rays.  While  not  directly  applicable  to  the  competition,  these 
speeches  did  help  broaden  our  knowledge  of  physics  and  prepared  us  for 
more  diversity  in  Olympiad  questions. 

The  U.S.  Team  was  quite  an  interesting  group  of  people.  We  had  a 
veteran  of  the  1991  Olympiad,  Eric  Miller,  who  was  always  looking  for  links 
to  things  he  had  done  before.  There  was  my  extremely  inquisitive  roommate, 
Dean  Jens,  the  very  verbose  and  confident  Mikhail  Sunitsky,  the  modest  (in 
attitude,  not  ability)  Victoria  Yung,  the  insightful  if  a  bit  soft-spoken  Szymon 
Rusinkiewicz,  the  most  sociable  Michael  Schultz,  and  the  extremely 
enthusiastic  Mary  Pat  Campbell.  I  don't  need  to  point  out  that  all  these 
people,  and  the  other  dozen  as  well,  were  very  competent  in  the  area  of 
physics  and  each  brought  their  own  perspective  to  the  problems  we  faced. 

Aside  from  our  daily  academic  routine,  we  took  the  time  to  visit 
Washington,  DC  for  half  a  day  during  training  camp.  Our  first  stop  was  the 
Old  Executive  Office  building,  where  we  met  with  former-President  Bush's 
Advisor  for  Science  and  Technology.  The  visit  included  a  speech  by  the 
Advisor  and  a  short  photo-op.  We  then  walked  over  to  the  National  Science 
Foundation,  where  things  began  to  get  interesting.  While  the  head  of  the 
Foundation  was  away  at  a  meeting,  we  did  get  to  talk  to  three  officials  who 
encouraged  us  to  discuss  the  status  of  science  education  in  the  United  States. 
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We  were,  they  said,  what  was  going  right  in  science  education,  and  hopefully 
we  could  tell  them  what  needed  to  be  done  in  science  education  today.  The 
most  striking  thing  about  this  conversation  was  that  there  was  no  program  or 
policy  that  we  could  point  to  as  our  inspiration.  We  all  noted  the  importance 
of  contests  such  as  the  one  we  were  representing  to  encourage  students  to  try 
harder.  Mostly  though,  we  considered  ourselves  the  result  of  the  dedication 
of  our  families  and  of  our  teachers,  as  individuals  not  as  part  of  some  broader 
program. 

At  the  Department  of  Education  which,  by  the  way,  shares  its  building 
with  NASA  as  a  prominent  symbol  of  its  importance  to  the  government,  the 
discussion  continued.  Once  again  the  aides  to  the  Education  Secretary,  Lamar 
Alexander,  were  sent  in  to  hear  our  comments.  Mr.  Alexander  emerged  only 
for  the  photo  and  quickly  retreated.  The  workers  in  the  Department  are 
proud,  I  remember  one  woman  there  who  was  happy  to  inform  me  of  the 
symbolism  of  the  Department's  seal.  Her  knowledge  of  that  piece  of 
information  for  some  reason  impressed  me.  During  our  discussion,  we  talked 
a  lot  about  what  could  be  improved  in  American  education.  Standard 
textbooks  were  debated,  as  was  how  much  the  Federal  Government  can  help. 
On  the  vast  majority  of  issues,  though,  we  agreed:  the  fundamentals  of 
science  must  be  taught  earlier,  the  middle  school  science  curriculum  should 
serve  as  a  bridge  from  the  experimenting  and  tinkering  of  elementary  science 
class  to  the  rule  based  systems  of  high  school,  elementary  physics  should  be 
taught  earlier,  a  bad  book  can  teach  little  or  disenchant,  and  above  all,  we 
have  to  change  the  system.  One  government  initiative  drew  particular  fire; 
the  America  2000  plan.  Not  one  of  us  expects  the  United  States  to  be  first  in 
the  world  in  science  and  math  by  the  end  of  the  millennium.  While  they 
argued  we  couldn't  aim  for  second  or  third,  we  made  it  clear  that  without 
more  reasonable  goals,  nothing  would  get  done. 

The  last  full  day,  the  team  members  were  announced  in  random  order. 
They  were  Szymon  Rusinkiewicz,  Michael  Schultz,  Dean  Jens,  Eric  Miller, 
and  myself.  The  alternate,  who  had  to  study  for  a  month  like  us,  but  without 
the  promise  of  a  trip,  was  Michael  Agney.  The  six  of  us  were,  naturally, 
ecstatic  and  received  the  support  of  the  other  students  immediately.  I 
personally  had  expected  to  be  in  the  second  group  of  five  people  at  the  camp 
and  so  I  was  pleasantly  surprised.  Of  course,  none  of  us  could  have  been  truly 
crushed  had  we  not  been  selected.  It  really  is  an  honor  just  to  be  on  the  team 
whether  or  not  you  represent  them  in  the  Olympiad. 

As  team  members,  we  were  all  benefactors  of  the  support  of  a  wide 
variety  of  corporations,  foundations  and  other  agencies.  The  U.S.  effort  was 
primarily  supported  by  AT&T  and  IBM  (Bell  Labs'  Dr.  Hauser  insisted  on 
alphabetical  order),  but  some  of  our  more  direct  support  was  in  the  form  of 
in-kind  donations.  We  received  a  wide  variety  of  textbooks  both  at  training 
camp  and  afterwards.  Even  last  December,  I  received  a  new  graphing 
calculator  form  Texas  Instruments. 

From  the  last  night  at  training  camp,  I  was  vividly  aware  that  our 
interaction  with  others  at  IPhO  would  be  in  some  way  representative  of  the 
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United  States.  We  had  instantly  acquired  the  need  for  diplomatic  courtesy  in 
all  our  dealings  as  participants.  That  night  I  made  a  mental  note  not  to  pack 
my  Tienanmen  Square  T-shirt  when  I  went  to  Finland.  There  were  very  few 
occasions  when  I  felt  constrained  by  this  diplomatic  courtesy^  but  they  existed. 
The  night  before  one  of  the  competitions,  the  Romanian  team  went  around 
offering  samples  of  a  folk  drink  from  a  wooden  flask.  We  were  warned  by  Eric 
and  Mike  to  expect  them.  When  they  did  visit  us,  we  could  not  refuse, 
though  rumor  had  it  that  the  proof  of  the  concoction  had  three  digits. 
Nonetheless,  we  all  had  about  half  a  mouthful.  Our  diplomacy  cost  merely  a 
head  rush  thirty  seconds  later. 

Between  training  camp  and  the  actual  competition,  the  real  training 
took  place.  To  train  for  a  physics  Olympiad  is  to  enter  into  a  degree  of 
seclusion.  Everyday  included  a  five  hour  exam,  checking  your  work,  and 
about  twice  a  week  there  was  a  lab  experiment.  I  took  the  time  to  enjoy  the 
festivities  surrounding  my  semi-graduation,  so  I  had  ail  of  the  experience  of 
training  in  less  time.  We  were  each  given  every  problem  from  the  previous 
twenty-two  olympiads,  each  harder  than  the  previous  one.  Each  day,  we  were 
asked  to  visualize  ourselves  in  the  competition,  going  from  problem  to 
problem  and  succeeding.  As  verification  we  sent  Arthur  Eisenkraft  a  three 
sentence  mantra  on  a  postcard  everyday.  The  third  part  of  this  training  was 
completion  of  the  freshman  physics  laboratory  course  at  a  local  university.  I 
performed  this  part  of  my  training  at  Case  Western  Reserve  University, 
where  my  high  school  physics  teacher,  Mr.  Shurtz  had  enough  authority  to 
make  it  easy  for  me  to  get  access  to  the  laboratory  program.  This  lab 
experiment  included  some  truly  interesting  experiments  including 
measuring  the  charge  to  mass  ratio  of  an  electron,  something  I  think  every 
physics  student  should  see  and  do  for  him  or  herself. 

The  night  before  I  left,  I  did  not  sleep  at  all  in  an  effort  to  adjust  to  the 
time  change.  If  it  had  gone  as  planned  I  would  have  slept  on  the  flight  and 
awoke  at  about  10  AM  Finnish  time.  It  didn’t,  and  I  ended  up  getting  the  first 
part  of  that  sleep  in  the  airport  in  New  York,  a  few  more  hours  on  the  plane 
and  the  remainder  when  I  fell  over  on  the  bed  in  my  hotel  room. 

We  arrived  on  the  third  of  July,  two  full  days  before  the  Olympiad 
convened  in  order  to  adjust  to  the  time  change  and  to  experience  a  little  bit  of 
Helsinki  for  ourselves.  After  giving  us  time  to  settle  into  our  rooms  (or  fall 
into  them  as  the  case  may  be),  we  were  taken  on  a  walk  around  Helsinki 
including  a  visit  to  the  embassy  district  and  the  front  of  the  U.S.  embassy.  We 
found  out  the  location  of  the  Independence  Day  festivities,  which  just 
happened  to  be  on  July  3rd,  and  left.  We  also  visited  a  small  snacks  restaurant 
and  saw  a  lot  of  the  Baltic  Sea  before  we  headed  back  to  our  hotel,  the  Ramada 
Presidentii. 

Our  only  celebration  of  Independence  Day  actually  on  July  4th  was  the 
lighting  of  a  pyramid  of  matches  in  symbolic  fireworks.  Eric,  perhaps 
remembering  what  we  had  just  seen,  insisted  on  having  a  full  glass  of  water 
ready  for  the  occasion.  What  we  had  just  seen  was  one  of  the  most  unusual 
cartoons  ever  produced  outside  of  Japan.  The  protagonist,  who  is  in  a  dream 
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sequence,  got  cold  in  his  house.  So,  naturally  he  built  a  fire  on  his  coffee  table. 
As  his  room  fills  with  smoke,  he  goes  outside  to  cut  a  hole  in  the  roof.  Now, 
unsatisfied  by  the  heat  generated,  he  proceeds  to  put  his  television  in  the 
flames  and  shortly  thereafter,  it  explodes.  I  won't  get  into  the  ending,  but  it 
nearly  as  confusing  as  the  Japanese  film  Akira. 

We  took  the  Helsinki  street  car  to  the  Independence  Day  Celebration, 
meeting  a  group  of  Canadians  living  for  two  years  as  missionaries  in 
Helsinki.  We  agreed  to  remain  silent  as  to  their  nation  of  origin.  We  arrived 
late  to  the  Seursaari  Show  Area,  site  of  the  festivities  and  missed  the  opening 
by  the  U.S.  ambassador.  On  the  other  hand,  we  missed  very  little  of  the  actual 
party.  Surrounded  by  other  Americans,  we  thoroughly  enjoyed  ourselves.  We 
took  part  in  the  volleyball  and  other  sports  and  took  home  a  sizeable  amount 
from  the  raffle. 

The  next  day,  we  took  the  officially  recommended  walking  tour  of 
downtown  Helsinki  visiting  sites  such  as  the  beautiful  Upensky  cathedral 
which  along  with  another  cathedral  and  the  parliament  building  dominates 
the  skyline.  There  are  no  sky  scrapers  in  Helsinki,  apparently  under  a  six  to 
eight  story  building  limit  which  makes  the  vast  majority  of  buildings  of 
uniform  height.  In  many  ways,  the  design  of  the  buildings  and  the  streets 
resembles  an  eight-story  Legoland  town  set  up.  If  the  streets  were  in  an 
American  city,  we  would  find  them  impassibiy  narrow,  but  the  Finns  don't 
drive,  they  walk  and  take  public  transport  and  they  have  a  great  system  with 
which  to  do  that. 

The  contest  itself  took  place  in  Espoo,  rather  than  Helsinki  proper. 
Espoo  is  Finland's  second  largest  city,  with  a  population  of  175,000,  and  is  just 
to  the  west  of  Helsinki.  We  stayed  in  the  easternmost  administrative  district 
of  Espoo,  Otaniemi,  home  of  the  Helsinki  University  of  Technology.  So,  the 
fifth  we  boarded  a  bus  and  moved  to  Espoo  where  we  were  to 
stay  in  the  Summer  Hotel  Dipoli  and  our  leaders  in  the  actual  Hotel  Dipoli, 
Unfortunately,  the  Finnish  government  had  other  plans. 

At  the  same  time  as  the  Olympiad,  another  gathering  of  many  nations 
was  taking  place  in  the  Helsinki  metropolitan  area,  the  Council  on  Securitv 
and  Cooperation  in  Europe  (CSCE).  While  I’m  sure  much  was  accomplished 
at  that  meeting  (including  the  expulsion  of  Yugoslavia  from  the 
organization),  the  only  thing  we  could  see  of  the  conference  were  the 
incredible  number  of  police  lines  strung  up  all  across  the  city  which  we  passed 
at  least  daily.  More  directly  affecting  us,  the  CSCE  preempted  us  from  our 
original  hotel.  As  a  result  the  team  leaders,  who  saw  the  questions  in 
advance,  were  staying  two  floors  above  the  team  members,  who  hadn't. 
While  not  truly  presenting  a  problem  for  anyone,  it  did  create  a  certain  degree 
of  anxiety  at  times  when  we  saw  our  coaches.  For  us  to  say  anything  more 
than  hi  would  have  potentially  aroused  suspicion.  Protocol  was  thus  kept  to 
a  minimum. 

In  keeping  with  international  tradition  we  brought  gifts  from  our 
home  cities  and  states  and  exchanged  them  with  other  teams.  We  each  built 
up  a  sizeable  collection  of  pins,  boomerangs,  maps  and  postcards  from  around 
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the  world.  To  be  honest,  the  boomerangs  were  only  from  Australia,  but  nearly 
everyone  brought  something.  It  all  made  the  occasion  very  festive  and  added 
to  experience  of  meeting  people  from  other  cultures. 

The  participants  at  Helsinki  were  clearly  representatives  of  the  post¬ 
cold  war  world.  In  the  opening  ceremonies,  we  sat  next  to  Ukraine.  Five  other 
participating  countries  had  been  independent  for  less  than  three  years.  The 
members  of  these  new  teams  were  very  open,  warm  and  friendly  people.  By 
contrast,  we  were  told  of  the  Soviet  teams  from  earlier  years.  They  remained 
austere  and  extremely  serious.  Like  us,  they  considered  the  competition 
important  and  strove  to  do  as  well  as  possible,  but  these  cold  war  teams 
thought  of,  or  were  told  to  think  of,  more  than  that.  They,  like  the  Soviet  bloc 
Olympic  athletes  contemporary  to  them,  were  representatives  of  a  political 
and  economic  system  on  the  world  stage.  Their  performance  was  a  triumph 
for  their  nation,  an  advance  in  the  cold  war.  Now  that  the  Soviet  bloc  has 
disappeared,  the  sole  austere  team  was  the  People's  Republic  of  China. 

For  the  first  few  days,  up  until  the  competition  was  over,  the  Chinese 
team  wore  white  dress  shirts  and  black  trousers  to  every  meal.  They  were  the 
picture  of  unity,  and,  at  least  to  me,  unapproachable.  This  is  no  impeachment 
of  their  character:  afterwards  they  were  as  friendly  as  any  other  team  there, 
joining  freely  in  the  gift  giving  which  so  often  marked  our  contact  with  those 
for  whom  the  language  barrier  was  too  great.  And  notably,  they  gave  us  a  slip 
of  paper  with  their  addresses,  ages  and  other  information. 

Of  the  178  participants,  only  four  were  female.  For  this  ratio  to  exist  in 
my  major  continually  disappoints  me.  We  must  work  to  interest  just  as  many 
young  women  as  men  in  physics.  Likewise,  of  the  thirty-nations  in 
attendance:  twenty-four  were  European;  three.  North  American;  a  pair.  South 
American;  and  seven,  Asian.  Even  granted  the  location  was  more  convenient 
for  the  European  teams,  especially  the  eastern  bloc  nations  for  whom  science 
is  rapidly  becoming  an  unaffordable  luxury,  the  geographic  distribution  of 
teams  is  in  need  of  improvement.  I  hope  that  the  Olympiad  in  Williamsburg 
will  include  nations  from  six  continents  rather  than  five,  all  of  the  most 
populated  nations,  and  will  inspire  students  of  all  nations  to  educate 
themselves  as  a  symbol  of  pride  in  their  homelands. 

The  United  States  had  the  single  largest  contingent  at  Helsinki, 
participating  both  as  representatives  of  a  nation  of  250  million  and  as  the  next 
hosts  of  the  competition.  Besides  the  five  of  us,  there  were  our  two  coaches, 
Avi  Hauser  and  Larry  Kirkpatrick,  textbook- writer  Anthony  P.  French,  A  APT 
representative  Delores  Mason,  the  coordinator  of  the  XXIV  Olympiad,  Arthur 
Eisenkraft,  a  high  school  teacher,  and  two  members  of  the  press  paid  by  the 
American  Institute  of  Physics. 

Helsinki  is  a  city  very  familiar  with  foreigners,  even  somewhat  shaped 
by  their  presence.  The  two  official  languages  of  Finland,  Finnish  and  Swedish, 
are  accompanied  by  English  and  German  nearly  everywhere  in  Helsinki, 
even  on  the  public  transit.  Downtown  Helsinki,  which  we  took  a  walking 
tour  of  the  first  two  days,  is  filled  with  governmental  centers,  commercial 
sites,  and  cultural  landmarks  integral  to  the  heart  of  Finland.  Nonetheless,  I 
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didn’t  feel  that  I  had  seen  real  Finnish  life  until  four  of  us  visited  our  guide's 
home  in  northern  Helsinki.  Although,  many  residential  districts  sprawled 
out  before  us  on  the  perspective  map  I  purchased  for  seven  dollars  from  the 
Finnish  national  tourism  organization,  the  life  of  the  Finns  would  have 
remained  unknown  to  me  had  I  not  made  this  visit. 

The  experimental  competition  was  conducted  in  two  flights.  The 
evening  before  the  experimental  competition,  we  learned  that  we  would  be  in 
the  A  shift  of  the  competition,  scheduled  for  8:00  AM  to  1:00  PM.  My 
roommates.  Dean  Jens  and  Szymon  Rusinkiewicz  took  this  news  with  som^ 
alarm,  and  as  I  showered  that  evening  made  preparation  to  go  to  sleep  by  8:45. 
They  would  awake,  they  assured  me,  at  6:15  with  plenty  of  tinae  to  shower 
and  get  ready  for  breakfast  at  7. 1  remained  awake  that  night  until  nearly  ten, 
entering  the  Dreaming  to  the  sound  of  Billy  Joel's  Vienna.  The  next  morning 
at  6:45,  I  awoke  and  glanced  at  my  watch.  When  I  told  my  sleeping 
roommates  what  time  it  was,  they  were  not  pleased.  Once  the  experimental 
session  started  (at  least  30  minutes  late),  we  were  divided  into  two  groups. 
Each  group  worked  on  one  experiment  at  first  and  then  we  switched.  This 
procedure  was  simple  enough  in  principle  but  in  practice  was  confused  in  at 
least  two  major  ways.  First,  the  descriptions  of  both  were  in  our  lab  books 
allowing  us  extra  time  to  prepare  for  the  second  experiment.  Second,  they 
decided,  or  were  forced  by  procedural  difficulties,  to  give  us  time  to  eat 
between  experiments.  Besides  putting  us  further  behind  schedule,  this  gave 
us  the  opportunity  to  prepare  for  the  2nd  experiment  and  for  some  teams  to 
discuss  that  experiment. 

Just  as  I  felt  less  nervous  about  the  experimental  competition  than  my 
teammates  did,  I  felt  calm  before  theoretical  competition.  Once  I  had  gathered 
the  materials  I  needed  for  it  (ruler,  calculator,  and  writing  utensils),  I 
removed  my  yearbook  from  the  Merimekko  bag  they  had  given  each  of  us 
and  turned  up  my  Walkman,  as  I  read  through  what  my  friends  at  school  and 
from  the  U.S.  Physics  Team  had  written  inside  it.  My  thoughts  of  the  conflict 
between  my  personal  feelings  and  the  chaos  of  getting  ready  for  the 
competition  which  surrounded  me  are  only  in  hindsight.  Then,  I  thought  of 
nothing  but  the  goodwill  scattered  in  the  formerly  blank  pages  and  margins  of 
my  yearbook. 

The  relations  between  participants  should  not,  in  theory,  reflect  the 
political  situation  of  their  nations.  After  all,  the  purpose  of  international 
interaction  at  this  level  is  to  foster  friendship  and  community.  Nonetheless, 
one  case  where  such  harmony  was  not  so  easy  quickly  comes  to  mind. 
Among  the  gifts  brought  by  the  Cypriot  team  were  pamphlets  from  their 
government,  both  advertising  the  island  nation's  tourism  potential  and 
condemning  the  1974  Turkish  invasion.  While  I  recognize  the  imperative  to 
make  such  statements  as  members  of  international  society,  it  seemed  that  the 
presence  of  the  Turkish  team  must  have  resulted  in  a  little  bit  of  trouble  on 
both  sides.  I  hoped  they  learned  not  to  condemn  each  other  for  their 
governments. 
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As  the  week  went  along,  I  went  to  sleep  progressively  later  and  later, 
excepting  the  night  before  the  experimental  competition.  One  of  the  things 
that  kept  me  awake  so  late  was  the  card  game,  Bartog  which  is  basically  Uno 
played  with  real  cards  and  with  one  twist.  After  each  round,  the  winner  gets 
to  make  up  one  rule.  In  a  good  game  there  is  a  fair  mix  between  verbal  rules 
and  procedural  rules  and  everyone  is  confused  beyond  belief,  but  someone 
always  remembers  a  rule  when  you  don’t.  We  had  one  or  two  really  good 
games  with  about  five  verbal  rules  which  made  each  play  of  the  cards  into  a 
speech. 

One  of  the  least  interesting  mandatory  activities  which  we  participated 
in  was  the  visit  to  the  Porvoo  works  of  Neste  Technology.  No  matter  how 
science  oriented  a  group  of  teenagers  is,  I  doubt  I  could  find  a  group  of  them 
who  would  revel  in  the  opportunity  to  visit  and  industrial  plant.  When  we 
first  arrived,  we  saw  a  film  designed  to  give  an  overview  of  Neste 
Corporation.  Three  economics  students  decide,  in  order  to  complete  a  project 
for  class,  to  make  an  energy  company  based  in  Scandinavia.  As  they  do  more 
and  more  research  they  refine  the  idea;  skimping  on  their  social  lives  to 
design  the  perfect  corporation.  What  they  decided  on,  of  course,  was  in  fact 
what  Neste  sees  itself  as,  a  growing  energy  supplier  in  a  strong  market, 
branching  out  into  other  fields.  As  a  corporation  introducing  itself,  the  film 
had  merit,  as  entertainment  it  needed  serious  help.  Bus  and  walking  tours, 
and  a  very  restful  presentation  on  solar  energy  completed  this  excursion.  The 
afternoon  was  to  redeem  our  day. 

That  afternoon  we  visited  Heureka,  the  Finnish  Science  Center  located 
in  Vantaa,  a  city  just  north  of  Helsinki  and  the  third  largest  city  in  Finland. 
There  was  the  typical  conglomeration  of  hands-on  science  activities. 
Commodore  Amiga  and  Macintosh  terminals  with  educational  programs, 
and  short  video  presentations.  In  true  form,  nearly  everything  was  translated 
into  both  Swedish  and  English.  While  the  featured  exhibit  was  a  celebration 
of  Finland  s  75  years  of  independence,  the  things  I  found  most  interesting 
were  elsewhere.  There  was  an  especially  large  area  devoted  to  language. 
Finnish  is  a  Finno-Ugric  language  similar  to  Hungarian  and  a  distant  relation 
of  Korean.  There  was  a  huge  globe  with  earphone  jacks  corresponding  to 
various  languages.  While  the  U.S.  speaker  says  "I  speak  English,"  the  one 
from  the  United  Kingdom  says,  "I  speak  British  English. "We  spent  3  hours 
looking  at  computer  simulations,  examining  World  War  I  artifacts,  having 
our  names  printed  in  Cuneiform,  and  searching  for  obscure  Finnish  coins. 

During  the  walk  on  the  first  day  I  noticed  the  ships  in  the  harbor  which 
went  back  an  forth  to  the  various  Baltic  nations,  serving  a  market  that  would 
not  have  existed  a  few  years  before.  The  Estonia  New  Line  boat  service 
advertised  a  one  day  visa-free  trip  to  Talinn,  Estonia  which  looked  very 
promising.  By  the  end  of  the  week  there  was  a  troop  of  four  of  us  willing  to 
go,  if  we  could.  Eric  Miller,  myself  and  two  members  of  another  team  all 
woke  up  at  some  insane  hour  and  proceeded  downtown.  That  morning,  as 
the  four  of  us  walked  to  the  offices  of  the  Estonia  New  Line,  I  had  a  definite 
feeling  that  our  trip  to  Talinn  was  not  to  be.  Twenty  minutes  later,  we  found 
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out  that  no  tickets  renaained  for  that  day.  Before  heading  back  to  Espoo,  I 
purchased  a  breakfast  of  strawberries  from  one  of  the  merchants  at  the 
Helsinki  central  market  whose  stand  was  less  than  fifty  feet  from  police  lines 
put  up  for  a  CSCE  motorcade.  The  most  enjoyable  part  of  the  experience  was 
the  contact  with  people  from  around  the  world,  something  no  day  excursion 
could  match,  although  Estonia  might  have  come  close. 

When  we  went  downtown  again,  I  spent  all  day  walking  around  with 
Michael  Schultz.  The  first  place  we  went  to  was  the  Forum  Mall.  Finally 
having  a  chance  to  see  the  European  Legos,  we  were  duly  impressed.  As  there 
was  a  conveniently  available  supply  of  technically  minded  people  available  at 
IPhO,  we  decided  it  would  be  a  good  idea  to  get  a  group  of  people  and  buy  one 
of  the  larger  sets  of  Technic  Legos  (formerly  known  as  the  Expert  Builder  sets), 
and  assemble  them  as  a  group  and  raffle  them  off.  Some  of  the  more 
advanced  models  even  included  a  programmable  interface.  Keeping  this  in 
mind  we  explored  the  city,  hoping  to  run  into  people  who  would  team  up  to 
buy  with  us  and  accept  a  lottery  among  the  builders  to  decide  who  should 
keep  it.  While  we  did  see  a  lot  of  Helsinki  that  day  including,  possibly,  the 
Chinese  restaurant  with  the  worst  quality  to  cost  ratio  on  the  planet.  We  only 
found  two  fellow  participants.  They  were  willing  to  buy  into  the  idea,  but  not 
to  pay  a  fourth  of  what  could  be  up  to  $180.  Taking  a  large  risk,  Mike  bought  a 
$150  set  anyway,  with  his  parents’  credit  card. 

Starting  that  night  at  dinner,  we  lobbied  for  builders  and  for  people 
willing  to  contribute  10  markka  for  a  shot  at  winning  the  completed  model. 
We  easily  found  a  group  of  about  15  builders  willing  to  work  on  it.  The 
number  of  people  buying  into  the  raffle  lagged  behind  the  number  needed  to 
break  even,  but  we  pressed  forwards  anyway.  With  so  many  people 
constructing  the  model,  which  was  a  working  pneumatic  bulldozer  and 
backhoe,  we  always  had  about  twice  as  many  people  as  could  successfully 
work  together  on  something  of  that  scale.  But  keeping  about  eight  people 
working  at  a  time  wasn't  a  problem  with  lack  of  interest,  overcrowding  and 
the  need  for  sleep.  We  completed  the  device  that  very  night  between  2:00  and 
3:00  AM. 

Regrettably,  as  a  direct  consequence,  our  team  was  less  than  fully 
prepared  for  the  closing  ceremonies  the  next  morning.  While  we  had  known 
our  awards  from  talking  to  our  leaders  the  previous  day,  the  experience  of 
going  down  to  the  stage  and  receiving  the  awards  was  just  as  impressive.  Eric 
and  Szymon  both  received  gold  medals;  Mike,  a  silver;  and  Dean  and  I  each 
got  honorable  mention.  The  United  States  had  done  the  best  it  had  ever  done 
in  the  International  Physics  Olympiad.  We  placed  fifth  among  the  nations 
participating  based  on  both  our  average  performance  and  our  medal  count. 
First  was  China  whose  team  garnered  five  gold  medals.  Of  course  these  team 
scores  are  unofficial,  teams  don't  receive  awards,  only  individuals  do. 
Through  the  day,  we  got  enough  participants  in  the  raffle  to  break  even,  at 
which  point  we  stopped  selling.  (Dur  last  buyer  was  Arthur  Eisenkraft  who 
gave  his  ticket  to  our  guide,  Amy  Wassenich  upon  realizing  the  potential 
problems  if  a  U.S.  team  leader  won.  We  scheduled  the  raffle  and  departed  for 
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downtown  Helsinki.  That  night,  our  coaches,  and  by  extension  all  the 
organizations  supporting  us,  treated  the  entire  United  States  entourage  to  a 
very  expensive  dinner.  While,  I  made  sure  not  to  convert  back  to  dollars 
until  I  went  home,  I  knew  it  was  a  lot.  The  meal,  which  averaged  about  $40 
per  person,  was  the  best  one  we  had  in  Helsinki  and  a  great  end  to  a  great 
week. 

We  somehow  managed  to  be  only  about  45  minutes  late  to  the  raffle 
and  collected  an  audience.  We  gave  a  hotel  employee  the  honor  of  selecting 
the  winner.  After  looking  up  the  name  next  to  the  number  on  the  list,  we 
found  that  the  winner  was  in  fact,  Amy  Wassenich.  I'm  sure  Amy  and  her 
fiancee,  who  joined  us  at  the  closing  ceremonies  have  already  had  hours  of 
fun  with  their  pneumatically  powered  bulldozer /backhoe  Lego  model.  That 
night,  I  got  as  many  people  as  I  could  to  sign  my  yearbook  and  didn't  get 
much  sleep.  My  return  trip  was  uneventful  excepting  the  last  twenty  minutes 
on  approach  to  New  York.  As  the  overhead  monitor  indicated  the  exact  flight 
path  of  the  plane,  we  all  watched  one  of  the  best  definitions  of  zig-zag  take 
shape  before  our  eyes.  And  then,  triumphant,  we  were  back  in  the  United 
States. 

Even  as  I  finish  talking  about  my  experiences  of  last  year,  many  of  the 
dedicated  people  who  made  that  U.S.  a  reality  are  working  on  hosting  this 
year's  Olympiad  in  Williamsburg,  Virginia.  With  many  other  countries  eager 
to  host  the  contest,  this  may  be  the  last  chance  to  see  and  help  out  with  such 
an  event  in  America  for  a  long  time.  I  encourage  all  of  you  to  get  involved 
and  to  contribute  to  yet  another  rewarding  experience  for  students  from 
around  the  world. 
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difficult  part  of  planning  is  accepting  that  it  must  be  done  in  spite  of  great 
and  ignorance.  The  business  plan  is  part  of  a  ''prediction/correction  state 
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Licensee  Questions 


Do  you  see  an  iinnsediate  need  in  the  market  place  or  is  market 
development  required? 

What  do  you  estimate  to  be  the  size  of  the  market? 

Estimate  the  market  duration. 

What  competition  do  you  see? 

Will  you  pursue  foreign  markets?  If  so,  what  is  your  experience  in 
foreign  markets?  What  experience  do  you  have  with  export  control 
regulations  and  for  what  countries? 

What  will  be  your  marketing  strategy? 

What  marketing  expertise  do  you  have  for  this  product? 

Will  the  product  be  distributed  separately  or  bundled  with  another 
product? 

When  do  you  anticipate  providing  the  product  for  sale? 

What  is  your  estimated  pricing  structure? 

What  maintenance  agreements  will  you  offer  to  your  customers? 

What  kind  of  user  support  will  you  provide,  such  as  training, 
newsletters  user  group  meetings  or  hot  lines? 

What  kind  of  license  do  you  seek?  (Exclusive,  nonexclusive,  area, 
use  or  hybrid. ) 

What  royalty  base  do  you  have  in  mind? 

Describe  your  current  and  future  projects  in  this  technology. 
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Information  on  Business  Funding 
(from  the  Wall  Street  Journal,  February  19,  1993) 

American  Bankers  Association 
800  872  7747 

Business  Consortium  Fund 
212  872  5590 

National  Association  of  Community  Development  Loan  Funds 
215  923  4754 

National  Association  of  Investment  Companies 
202  289  4336 

/ 

National  Association  of  Minority  Contractors 
202  347  8259 

National  Bankers  Association 
202  588  5432 

National  Business  Incubator  Association 
614  593  4331 

National  Council  of  State  Legislators 
303  830  2200 

National  Federation  of  Independent  Business 
800  552  6342 

National  Minority  Supplier  Development  Council 
212  944  2430 

National  Small  Business  United 
202  293  8830 

Overseas  Private  Investment  Corporation 
202  336  8400 

Small  Business  Legislative  Council 
202  639  8500 

U.S.  Office  of  Small  and  Disadvantaged  Business/ AID 
703  875  1551 

U.S.  Small  Business  Administration 
800  827  5722 
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So  what  is  the  problem? 
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-  Your  market  is  that  portion  of  the  population  or  of  businesses  that  you  are  offering 
product  to.  The  function  of  marketing  is  to  uncover  the  true  purchaser  of  your  good 
make  them  your  market. 
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-  A  small  business  needs  sales  persons  far  more  that  marketing  employees,  but  it  is  critical 
that  thought  is  given  to  where  your  sales  effort  should  be  directed. 


Promotion  of  your 
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media  designed  for  impulse  buying 
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-  Listen  to  how  people  approach  you  over  the  telephone  when  they  are  selling  to  you.  You 
will  find  there  are  a  wide  range  of  talents  out  there. 


Selling  Continued 
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Telephone  Survey 
STATUS : 


Date: 


Time: 
Time: 

Not  a  Lead 


Immediate  Follow  Up  _ 

Call  Back 


Appt  Date: 
on  Date : 


Prospect's  Ncime 
Phone  Number 


Irlc-lt**********  **********************5lr********Hr********************** 

Hello.  This  is  _ of  BLES  Scientific.  We  are  conducting  a 

survey  of  manufacturing  firms  in  this  area  and,  if  I  could,  I'd 
like  to  ask  you  a  couple  of  questions.  Do  you  have  two  minutes  to 
answer  these  questions? 

:if  no;  When  would  be  a  good  time  to  call  back? 

1.  What  type  of  products  do  you  manufacture  at  your  facility? 

Answer : _ 

2.  Are  you  currently  using  computers  to  assist  you  in  your 
manufacturing  efforts? 

:  if  no;  Let  us  send  you  a  brochure  that  describes  the  benefits 
you  can  gain  from  our  production  support  software.  The  information 
may  help  you  when  you  configure  your  new  system. 

:  if  yes;  What  kind  of  hardware?  _ 

What  kind  of  software?  _ _ 

Is  your  system  meeting  your  needs?  _ : _ 


3.  Are  you  thinking  of  getting/upgrading  your  software  capability 

in:  0-3  months  _ 

3-6  months  _ 

6-12  months  _ 

4.  Thank  you  very  much  for  your  time.  Good  bye. 

*************  *!Hf******:fr****:*r***********:«r**************************** 


Immediate  Impressions: 


Hot 

Possible 

Cold 
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-  Set  bold  goals  but  derive  conservative  plans  to  reach  them.  This  is  ray  resolution  to  all 
the  words  I  have  heard  and  schemes  I  have  tried. 
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Properties  of  the  Solutions  to  the 
Time— Independent,  Nonlinear  Cubic,  Schrodinger  Equation 

Ronald  E.  Mickens 
Clark  Atlanta  University 
Department  of  Physics 
Atlanta,  Georgia  30314,  USA 

Abstract 

We  study  the  mathematical  properties  of  the  time-independent,  nonlinear  cubic, 
Schrodinger  equation.  This  ordinary  differential  equation  has  complex-valued  solutions 
and  corresponds  to  two  coupled,  second-order  differential  equations.  Using  a  polar  co¬ 
ordinate  representation,  we  show  that  the  radial  equation  reduces  to  a  one-dimensional 
conservative  classical  problem  for  which  the  potential  has  edl  bound  states.  A  knowledge  of 
the  radial  function  then  allows  the  direct  calculation  of  the  phase  fxmction.  Approximations 
to  the  solutions  axe  obtained  by  means  of  a  perturbation  procedure. 

I.  Introduction 

The  time-dependent,  nonlinear  cubic  Schrodinger  equation  (TDNSE) 

/,\  d^u  ,  ,0  ,  .s 

(1)  =  +  “  = 

is  ubiquitous  in  the  sciences  and  engineering^ In  particular,  this  equation  provides  a 
mathematical  model  for  shallow  water  waves^,  the  transport  of  vibrational  energy  along 
a-helix  proteins  through  the  formation  of  solitons^,  and  nonlinear  wave  propagation  in 
optical  fibers.^  Our  purpose  is  to  investigate  the  properties  of  the  solutions  to  the  time- 
independent,  nonlinear  cubic,  Schrodinger  equation 

(Pu 

(2)  =  0,  u  =  u{x). 

For  futiure  references,  Eq.  (2)  will  be  denoted  as  TINSE:  the  time-independent,  nonlinear 
Schrodinger  equation.  While  a  great  deal  of  information  exists  on  the  various  solutions  of 
Eq.  (1),  in  particular  its  soliton  solutions^’^,  little  is  known  concerning  the  properties  of 
the  solutions  to  the  TINSE. 
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Certain  preliminary  results  on  Eq.  (2)  have  been  obtained  by  D.  Lockett  and  presented 
in  his  Master-of-Science  thesis.®  There,  he  used  the  method  of  harmonic  balance®  to 
calculate  an  analytic  approximation  to  the  general  solution  of  Eq.  (2).  In  this  paper,  we 
extend  the  investigations  of  Lockett  and  study  the  general  mathematical  properties  of  u(x). 

The  next  section  provides  certain  background  information  needed  for  our  analysis. 
Section  III  gives  the  polar  representation  of  the  solution  aind  shows  how  exact  first -integrals 
can  be  calculated.  In  Section  IV,  we  give  two  special  exact  solutions  to  the  TINSE.  They 
follow  from,  respectively,  requiring  the  amplitude  and  phase  functions  to  be  constant. 
Section  V  is  devoted  to  a  detailed  study  of  the  mathematical  properties  of  the  solutions  to 
Eq.  (2).  We  find  that  much  can  be  learned  about  these  solutions  in  spite  of  the  fact  that  an 
explicit  representation  for  them  does  not  exist.  In  Section  VI,  we  show  that  approximate 
solutions  for  the  TINSE  can  be  calculated  using  a  perturbation  procdure. 

II.  Preliminaries 

Examination  of  Eq.  (2)  shows  that,  in  general,  u(x)  is  a  complex-valued  fvmction  of 
the  real  variable  x.  The  solutions  to  Eq.  (2)  can  be  written  in  either  the  cartesian  or  polar 
representations: 

(3)  u{x)  =  X{x)  +  iY{x), 


(4) 


«(x)  = 


where  X,  Y,  R  and  9  axe  real  functions  of  x.  Substitution  of  these  results  into  Eq.  (2)  and 
setting  the  real  and  imaginary  parts  equal  to  zero  gives 


X"  +  {X^  +  Y^)X  =  0, 

+  +  =  0, 


and 


(6) 


1  Re"  +  2R'9'  =  0, 
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where  the  prime  (^)  denotes  tahing  the  derivative  with  respect  to  x.  These  representations 
indicate  that  the  TINSE  is  equivalent  to  two  second'-order,  coupled  differential  equations. 
For  the  remainder  of  this  paper  only  the  polar  forms  given  by  Eqs.  (4)  and  (6)  will  be 
considered.  Also,  with  no  loss  of  generality,  the  following  initial  conditions  will  be  selected. 

(7)  i2(0)^0,  it:'(0)  =  0;  6(0)  ^0,  ^'(0)^0. 

IIL  Details  of  the  Polar  Representation 

Equations  (6)  can  be  integrated  once  to  obtain  first-integrals.  First,  note  that  the 
second  of  Eqs.  (6)  can  be  integrated  to  give 

(8)  R^d'  =  A  =  [il(O)]20'(O), 

where  the  integration  constant  A  is  expressible  in  terms  of  the  initial  conditions.  Second, 
solving  Eqs.  (8)  for  O'  and  substituting  it  into  the  first  of  Eqs.  (6)  gives 

(9) 

This  equation  has  the  following  first-integral 

(10)  ^^l^  +  U{R)  =  E=  constant  >0, 
where 

(11)  =  + 

Note  that  R{x)  can  be  considered  to  be  the  dependent  variable  describing  a  one- 
dimensional  dyneunic  system  having  the  convex  potential  energy  function  (/{RY’^;  conse¬ 
quently,  Eq.  (10)  is  the  energy  integral.  Since  both  terms  on  the  left-side  of  Eq.  (10)  are 
non-negative,  it  follows  that  for  i?  >  0,  all  the  integral  curves  in  the  phase-plane  (R,  R') 
are  closed.^  Hence,  R{x)  is  a  periodic  function  of  x. 

The  following  comments  are  in  order  at  this  point:  (i)  Examination  of  Eqs.  (9), 
(10)  and  (11)  show  that  the  radial  equation  is  equivalent  to  a  one-dimensional  classical 
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nonlinear  oscillator  problem,  (ii)  From  the  definition  of  the  radial  function  i2(r),  it  follows 
that  R{x)  >  0.  For  the  potential  function  U{R)y  this  condition  is  seen  to  hold  for  any 
finite  energy  E  >  0,  (iii)  The  potential  function  U{R)  has  a  single  minimum  located  at 

(12) 

for  which  it  has  the  value 

1^1“/^  >  0. 

The  phase  function  d{x)  can  be  determined  once  is  known.  From  Eq.  (8),  we 

have 

""  [iZ(x)]2‘ 

The  solution  of  this  equation  is 

In  spite  of  the  fact  that  R{x)  is  periodic,  i.e.,  the  integral  of  a  periodic  function  does  not, 
in  general,  give  a  periodic  function. 

IV.  Exact  Solutions 

Equations  (6)  have  two  special  exact  solutions.  The  first  is  obtained  by  requiring 

(16)  R{^)  =  constant. 

The  solution  for  this  case  is 

(17)  u(x)  =  .R=lA|'/^ 

and  correspnds  to  the  “motion”  of  the  one— dimensional  system  being  confined  to  the 
bottom  of  the  potential  well  U{R)  at  R  =  R. 

The  second  solution  is  gotten  by  requiring 

(18)  ^(x)  =  9o  =  constant. 
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This  condition  gives 


(19)  u(x)  =  e'^°i?(0)c„[i?(0)x;  l/>/2], 

where  c„['  •  •]  is  the  Jacobi  elliptic  cosine  function®. 

Both  of  these  special  solutions  are  periodic.  In  each  case,  the  period  depends  on  the 
initial  conditions.®  It  is  unlikely  that  other  special  exact  solutions  exist  to  the  TINSE.  In 
any  case,  the  general  solution  cannot  be  expressed  in  terms  of  a  finite  combination  of  the 
elementeiry  functions.  In  Section  VI,  we  show  how  to  determine  approximations  to  the 
exact  solutions  of  the  TINSE  for  “motions”  near  the  bottom  of  the  potential  well. 

V.  Mathematical  Properties  of  u{x) 

Theorem  1.  R(x)  >  0. 

Proof.  This  follows  directly  from  Eqs.  (11)  and  (13). 

Theorem  2.  R{x)  is  periodic  with  period  x. 

Proof.  See  the  arguments  immediately  following  Eq.  (11).  Note  that  while  the  value  of 
z  is  not  known,  we  can  be  sure  that  it  does  exist. ^’®  (Section  VI  will  provide  a  means  to 
calculate  its  value  using  a  pertmrbation  procedure.) 

Theorem  3.  R(x)  can  be  represented  as  a  “cosine”  Fourier  series.  This  means  that  R(x) 
is  an  even  function,  i.e., 

(20)  R{-x)  =  R(x). 

Proof.  This  result  is  a  direct  consequence  of  the  selection  of  initial  conditions  for  R(x) 
as  given  in  Eq.  (7). 

Theorem  4.  6{x)  is  an  odd  function,  i.e., 

(21)  0(-x)  =  -e(x). 
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Proof.  From  Eq.  (15),  we  have  that 


(22) 

Now  let  z  =  —t.  Therefore, 


dt 

IW' 


(23) 


e{-x)  =  A  f 

Jo 


i-dz) 


-e{x). 


Theorem  5.  u{—x)  =  [■u(x)]*. 

/ 

Proof.  The  (*)  denotes  complex-conjugation.  From  Eq.  (4),  it  follows  that 

(24)  u(-x)  =  R(-x)e'^^-‘\ 

However,  Theorems  3  and  4  state  that  i?(x)  is  even  and  ^(x)  is  odd.  Consequently,  using 
these  fcicts  in  Eq.  (24)  gives 

(25)  u(-x)  =  i2(x)e-‘^(*)  =  [u(x)]*. 


Theorem  6.  There  exists  a  constant  cq  such  that 

(26)  u(x  +  x)  =  e‘"»*u(x). 

Proof.  As  previously  defined,  x  is  the  period  of  the  periodic  function  R(x),  i.e.,  i2(x+x)  = 
R(x).  Let 

(27)  u  =  —. 

X 

Then  from  Theorem  3,  it  follows  that  i?(x)  hais  the  Fourier  representation 

oo 

(28)  >  R(x)  =  4-  cos(kux). 

ib=i 

Since  i?(x)  >  0  (Theorem  1),  there  exists  coefficients  {c/a:}  such  that 

(29)  =  J  +  E  *  cos(ku,x). 


105 


Hence,  from  Eqs.  (15),  it  follows  that 


(30) 

where 

(31) 

(32) 

and 

(33) 

Consequently, 

(34) 


0(x)  =  A 


^ + E  (^) 


=  CqX  +  ^l(x), 


Co  = 


Ado 


<50  /  J  \ 

di{x)  =  A  ^  f  ]  sin(A:wx), 

fc=i  ' 


9\{x  +  x)  =  ^i(x). 


u{x)  =  i?(x)e‘^(*)  =  i2(x)e‘t'=«==+®‘(^)J, 


and 


(35)  u{x  +  x)  =  R{x  +  =  J?(x)e'‘'(^)e‘‘^°^  =  e'"°^u(x). 


Theorem  7.  u(x)  has  the  representation 

(36)  tz(x)  =  $(x)e"^®^ 
where 

(37)  <&(x  +  x)  =  $(x). 

Proof.  This  result  follows  from  Theorem  6  where  Eqs.  (26)  is  treated  as  a  finite  difference 
equation.®’^® 

VI.  An  Approximate  Solution 

While  Eq.  (9)  cannot,  in  general,  be  solved  exactly,  we  can  calculate  good  approximate 
solutions  for  motions  near  the  bottom  of  the  potential  well  at  R  =  R. 
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Define  the  parameter  a  to  be 

(38)  a  =  R{0)-R, 
and  the  new  variable  y(x)  as 

(39)  y{x)  =  R{x)  -  R. 
We  assume  that 

(40)  \a\  <  R. 
From  this  it  follows  also  that 


(41)  ly(x)l  <  R. 

Substitution  of  Eq.  (39)  into  Eq.  (9)  gives,  after  a  large  amount  of  algebraic  manipulation, 
the  result 


(42)  y"  +  [6|A|2/3]y  _  [3|A|i/3]y2  +  Hy^  =  0, 

where  only  terms  through  order  have  been  retained.  The  initial  conditions  for  y(x)  are 

(43)  y(0)  =  a,  y'(0)  =  0. 

Equations  (42)  and  (43)  are  an  example  of  a  mixed-parity  perturbation  problem^  ^ 
and  unless  one  is  very  careful  difficulties  can  arise  in  the  construction  of  a  uniformly  valid 
solution.  In  a  recent  paper,  Mickens  ha^  shown  how  to  resolve  the  various  issues  related  to 
having  the  perturbation  solution  satisfy  exactly  the  initial  conditions  of  Eq.  (43).  Using 
these  results  of  Mickens^ ^  the  following  perturbation  solution  is  obtained  for  Eqs.  (42)  and 
(43): 


R{x)  =  R  +  a  cos(u;x) 
^2  1 


126il 


(44) 


+ 


[—3  +  2  cos(a;x)  +  cos(3a;x)] 
'43,311 


[I323(i?)2j  L 


-1  - 


576 


cos(wx)  +  (  X  )  cos(2u;x) 


+  r^|^)cos(3wi)|  +0(<.<), 
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where  the  parameter  a  plays  the  role  of  the  perturbation  expansion  parameter.  The  angular 
frequency  is  given  by  the  expression 


(45) 


LJ 


,  /14,533\  a2  ,  4, 

(20,168)  (^)2 


The  phase  function  d{x)  can  be  calculated  by  direct  substitution  of  Eq.  (44)  into  the 
right-side  of  Eq.  (15).  We  do  not  give  tlois  calculation,  but,  leave  it  as  an  exercise  in  Taylor 
series  expansion  and  elementary  integration  for  the  reader  to  complete. 


VII.  Summary 

We  have  shown  that  many  of  the  important  mathematical  properties  of  the  TINSE 
C2in  be  determined  by  use  of  a  polar  representation  for  the  solution.  The  theorems  of 
Section  V  summarize  the  results  found  to  date. 

Note  that  u(x),  in  general,  is  not  periodic.  This  follows  from  the  relation  given  by 
Eqs.  (36)  £ind  (37).  There  are  two  periodic  behaviors  associated  with  u(x).  The  first  is 
related  to  the  fact  that  the  radial  function  R{x)  can  be  considered  as  the  solution  to  a 
one-dimensional  classical  nonlinear  oscillator  problem;  this  period  Ti  =  x  and  its  value 
depends  on  the  initial  conditions.  The  second  period  is  determined  by 

(46) 

Co 

If  T1/T2  is  a  rational  number,  then  u{x)  is  periodic,  otherwise,  u{x)  is  queisi-periodic  or 
almost  periodic.  In  general,  for  arbitrary  initial  conditions,  we  expect  the  latter  situation 
to  hold. 

It  is  of  interest  to  see  what,  if  any,  relationships  exist  between  the  solutions  of  the 
TDNSE  (1)  and  the  TINSE  (2).  We  are  currently  investigating  this  issue. 
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Inteeration  of  a  large  number  of  components  on  a  single  chip  requires 
sophisticaSd  interconnections  to  minimize  signal  delays  and 
ont^ze  the  packing  density.  Aluminum  has  been  most  widely  used  for  contacts 
SdTnteroXSs  in  both  bipolar  and  MOS  integrated  ‘eclmo lo^^  Low 

temperature  interditfusion  of  aluminum  and  siliconlil  j 

nnci^vation  or  packaging  of  the  device  can  result  in  gain  (kgration  and  mcreasea 

function  leakage.  The  use  of  A1  films  containing  1%  to  2%  Si  has  been 

minimize  contact  pitting  BI  but  precipitation  and  f 

doped  Si  in  the  contact  window  du^g  the  cooling  f 

increases  the  contact  resistance  of  small-area  contacts  drastically,  uev 

reliability  can  be  improved  by  interposing  a  barrier  layer  between  the  A1  and  the  Si 

which  reduces  the  mass  transport  in  the  contact  structure  during  processing  o 

operation  of  the  device. 

With  the  evolution  of  VLSI  technology  into  the  submcron  dimenmon  therms 
an  increase  in  the  proximity  of  layers  of  different  materials  and  arguments  the  role 
of  sTmcturhl  difects  inherent  in  these  layers.  This  increases  the  possibility  of 
l“«actts  in  multilayer  thin  film  assemblies.  Therefore 
number  of  studies  are  devoted  to  the  development  of  smtaWe  barrier  layers 
reduce  or  prevent  such  materials  reactions.  In  spite  of  these  efforts,  rules  r^^fi^st 
Teen  foldto  date  which  would  permit  to  select  a  barraier  fy®*- 
SiplTs  However,  general  trends  in  the  barrier  properties  of  materials  used  in 
Lntac't  structures  to  silicon  semiconductor  devices  ‘>ocome  apparent  Smilar 
trends  can  be  attributed  to  materials  belonging  to  the  same  materials  c 

REQUIRMENTS  FOR  BARRIER  LAYERS 

The  general  approach  to  form  an  electrical  contact  to  a  silicon  semiconductor 
device  is  a  multilayer^ thin  film  structure.  It  consists  of  a/onabination  of  different 
mlJerials  which  Ive  bound  from  each  other  by  interfaces.  In  addition,  the 
SdltTdual  aims  contain  structural  defects,  such  as  gram  be-rndams  ^d 
dislocations.  Unfortunately,  such  an  assembly  'S by 

eauilibrium  Therefore,  atom  migration  will  take  place  to  establish  7 

lowering  the  total  free  energy  of  the  assembly.  The  driving  force  for  ^e  atom  c 
motion  fs  the  gradient  of  the  chemical  potential  of  the  atoms  m  the  assembly. 

BARRIER  FOR  CHEMICAL  INTERDIFFUSION 

The  net  flow  of  atoms  is  commonly  described  as  diffusion  and  charactenzed 
by  a  diffusion  coefficient  D.  In  a  multi-layer  thin  film  system  the  combined  effect 
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several  diffusing  species  is  expressed  by  a  chemical  interdiffusion  coefficient 


D  s  X  xi  Di 

« 

1 

where  xi  and  Di  are  the  atomic  fraction  and  diffusitivity,  respectively,  of  the  species 
i.  It  is  D  which  is  determined  in  an  ordinary  diffusion  experiment.  The 
temperature  dependence  of  D  is  generally  written  in  the  form  D  =  Do  exp  (-Ea/kT), 
which  is  known  as  the  Arrhenius  equation.  Ea  is  the  experimental  activation 
energy  and  Dq  is  the  frequency  factor.  In  the  case  of  self-diffusion  and  diffusion  in 
dilute  substitutional  solid  solutions  a  specific  meaning  in  terms  of  atomic  processes 
can  be  attributed  to  the  quantities  Ea  and  Dq.  However,  in  the  general  case  of 
chemical  interdiffusion  the  nature  of  the  activation  energy  becomes  complex. 


Structural  defects  have  a  profound  effect  on  the  diffusion  rate  of  atoms  in 
thin  films.  From  a  large  number  of  experimental  measurements  it  is  generally 
concluded  that  diffusion  at  free  surfaces  is  the  quickest.  Diffusion  along  grain 
boimdaries  and  dislocations  is  the  intermediate,  and  diffusion  in  interior  of  crystals 
is  the  slowest  process  in  thin  film  diffusion  mechanisms.  Each  of  these  process  have 
a  different  activation  energy  which  results  in  different  temperature  behavior.  In  a 
polycrystalline,  thin  film  lattice  diffusion  is  dominant  at  high  temperatures.  The 
reason  is  that  the  total  cross-sectional  area  of  the  grain  boundaries  is  small 
compared  to  that  of  the  gfrain  boundaries,  resulting  in  a  small  contributions  of  the 
grain  boundaries  to  the  overall  diffusion  mechanism.  At  low  temperatures, 
however,  processes  of  low  activation  energy  become  progressively  important.  Thus, 
grain  boxindary  diffusion  dominates  at  low  temperature.  In  addition,  the  smaller 
the  grain  size,  the  greater  the  importance  of  boundaries  in  the  diffusion  process. 
The  crossover  of  the  dominant  mechanism  from  lattice  to  grain  boundary  diffusion 
in  metals  occurs  at  about  171  to  3/4  of  the  melting  point  temperature. 

Naturally,  one  would  prefer  a  single  crystal  material  as  barrier  layer  against 
chemcal  interdiffusion,  because  lattice  diffusion  has  the  lowest  diffiisivity.  But  for 
device  processing,  single  crystal  or  epitaxial  layers  are  impractical,  because  they 
can  rarely  bef  formed  in  a  contact  structure.  Therefore,  the  next  best  choice  for  a 
diffusion  barrier  layer  is  a  large-grained  polycrystalline  material.  In  addition, 
solutes  ^  segregated  at  the  grain  bondaiies  may  decrease  the  grain  boundary 
diffiisivity.  However,  this  is  not  a  general  rule,  since  it  has  been  observed  that  the 
converse  can  also  occur,  For  this  reason,  the  concept  of  segregating  impurities  or 
second  phase  particles  at  the  grain  boundaries  of  a  polycrystalline  barrier  matterial, 
does  not  necessarily  yield  an  efficient  diffusion  barrier. 

Another  class  of  possible  diffusion  barrier  materials  are  amorphous  thin- 
films.  The  diffiisivity  in  these  materials  is  intermediate  between  lattice  and  grain 
boundary  diffusion.  19-101  However,  due  to  the  very  limited  data  available,  the 
diffusion  mechanism  in  amorphous  materials  is  poorly  understood.  More  studies 
are  needed  to  evaluate  the  usefulness  of  amorphous  materials  as  diffusion  barriers. 
It  should  be  noted  that  diffusion  barriers  do  not  eliminate  the  driving  force  for 
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diffusion,  also  known  as,  the  concentration  gradient.  They  merely  reduce  mobility 
for  atomic  migfration.  The  efficiency  of  a  diffusion  barrier  material  is  thus 
determined  by  how  well  it  diminishes  mass  transport  through  its  diffusion  paths 
present. 

BARRIER  FOR  CHEMICAL  INTERACTION 

The  total  free  energy  barrier  of  a  thin-film  assembly  can  also  be  lowered  by  a 
chemical  reaction  resulting  in  the  formation  of  a  new  phase.  The  interfacial 
reaction  between  two  dissimilar  materials  generally  requires  mass  transport  across 
the  intermatter  of  kinetics.  The  formation  of  a  new  phase  is  a  matter  of  nucleation 
and  growth.  For  a  new  phase  to  nucleate  the  free  energy  of  an  atom  in  the  new 
phase  must  be  lower  than  that  in  the  inital  phase. 

A  multilayer  thin-film  system  is  in  stable  chemical  eqinlibrium  as  long  as  any 
change  in  temperature,  concentration,  or  pressure  does  not  decrease  the  total  free 
energy  of  the  system.  If  the  total  free  energy  is  reduced,  a  spontaneous  chemical 
reaction  may  occur.  A  thermodynamic  criterion  for  a  chemical  reaction  to  take  place 
at  a  given  temperature  is  a  negative  Gibbs  free  energy  change  for  the  reaction  G  = 
AH  -  T  AS,  where  H  is  the  change  in  enthalpy  and  S  is  the  change  in  enthropy.  In 
many  cases  the  contribution  form  S  is  small  and  the  reaction  is  mainly  determined 
by  H.  Thus,  the  main  driving  force  for  an  interfacial  reaction  is  a  difference  in 
chemical  free  energy. 

Barrier  layers  can  be  used  to  reduce  the  total  free  energy  of  a  thin-film 
system.  Two  thin-films  A  and  B  in  contact  with  each  other  have  a  tendency  to  react 
and  form  a  compound  AB  if  the  fee  energy  of  the  system  compared  to  the  A/B 
configuration.  By  the  addition  of  the  layer  C,  the  system  becomes  more  stable.  The 
layer  C  has  therefore  been  named  reaction  barrier. 

In  general,  the  barrier  layer  C  does  not  establish  chemical  equilibrium  in  the 
system  A/C/B.  An  interfacial  reaction  will  therefore  still  occur,  but  at  a  reduced 
rate.  This  reaction  consumes  the  layer  C  prevent  reactions  between  the  layers  A 
and  B.  This  situation  is  illustrated  in  Figure  2.  Contact  structure  in  Figure  2 
consists  of  a  silicide  layer,  a  transition  metal  layer  as  reaction  barrier,  and  the 
aluminum  overlayer.  An  increase  in  temperature  during  processing  or  operation  of 
the  device  will  initiate  interfacial  reactions  between  the  aluminum  and  the 
transition  metal  layer,  and  the  transition  metal  and  silicide  layer.  The  reaction 
products  are  an  aluminum  intermetallic  compoxmd  and  a  transition  metal  silicide. 
In  many  cases,  the  reaction  rates  can  be  measured  experimentally.  Provided  that 
the  reaction  is  laterally  uniform,  the  point  in  time  when  they  transition  metal 
barrier  layer  is  consumed  by  the  interfacial  reaction  is  then  pre^ctable  for  a  given 
temperature. 

Thus,  a  minimum  barrier  layer  thickness  can  be  found  for  any  given  time- 
temperature  cycle.  Upon  consumption  of  the  barrier  layer,  rapid  interdiffusion 
generally  causes  a  fast  deterioration  of  the  contact  structure.  This  concept  of  a 
consumable  barrier  layer,  which  has  been  termed  "sacrificial  barrier"  by  Nicolet, 

[11, 12]  can  only  be  applied  if  atomic  diffusion  through  the  barrier  layer  is  negligible. 
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RESISTIVITY  AND  ADHESION  REQUIREMENTS 

Barrier  layers  have  to  fulfill  two  other  requirements  for  device  application, 
namely,  low  contact  resistance  and  good  adhesion.  It  is  clear  that  the  hairier  layer 
should  have  metallic  properties.  However,  a  low  resistivity  is  generally  not  of 

primary  concern.  A  typical  thickness  for  barrier  layers  is  lOOOA.  For  a  IpD  cm2 
contact  area,  a  barrier  layer  material  of  lOOfiO  cm2  resistivity  would  only  contribute 

0.1Ocm2  to  the  total  contact  resistance.  More  important  is  a  low  contact  resistivity 
when  the  barrier  layer  is  used  as  an  Ohmic  contact  to  silicon.  For  shallow  junction 
devices,  contact  resistivities  in  the  10"®  to  low  10"'^Dcm2  range  are  required. 

Barrier  layer  material  should  be  compatible  with  device  processes,  such  as 
patterning,  passivation,  multilevel  metallization,  and  packaging.  Good  adhesion  to 
insulators  and  corrosion  resistance  in  various  gas  ambients  may  also  be  required. 

ALUMINIUM-ALLOY  METALIZATIONS: 

Simple  single  layer  aluminium  metallizations  are  unstable  even  at  moderate 
temperatures.  Grain  degradation  or  Junction  shortening  is  a  result  of  low 
temperature  •  interdiffusion  of  Aluminium  and  Silicon  during  contact  intering, 
passivation  or  packaging.  The  aluminum  tends  to  mound  up  into  hillocks  during 
post  disposition  anneal  cycles  degrading  the  texture  and  reflectively  of  the  films 
causing  patterning  problems.  Grain  growth  at  400  degrees  C  is  a  result  of  the  low 
melting  point  of  aluminum  and  high  solubility  of  silicon. 

Reduction  of  interdiffusion  of  in  the  Contact  System  can  be  reduced  by 
interposing  a  barrier  layer  in  between  the  aluminum  and  the  silicon. 

An  example  of  a  Barrier  Metal  "  is  titanium.  The  structures  consisting  of 
aliuninum  as  an  interconnect  refractory  metal,  titanium  as  the  diffusion  barrier, 
and  Silicon  as  the  main  contact  material  imporves  the  contact  stability  but  only  up 
to  400OC.  At  400^0,  aluminum  and  titanium  start  to  react  with  each  other.  They 
form  a  new  compound  TiA13.  As  the  temperature  exceeds  400  degrees  C,  the 
mununuin  will  pBnotrate  into  tho  Silicon  substrate.  The  end  result  of  this  is  a 
binary  and  ternary  compound.  Figure  I  illustrates  the  effect  of  the  phase 
transformations  in  the  vicinity  of  the  contact  hole.  First,  the  TiA13.  As  the 
temperature  exceeds  400  degrees  C,  the  aluminum  will  penetrate  into  the  Silicon 
substrate.  The  end  result  of  this  is  a  binary  and  ternary  compound.  Figure  I 
Ulustrates  the  effect  of  the  phase  transformations  in  the  vicinity  of  the  contact  hole. 
First,  the  TiA13  formation  be^ns  to  across  the  entire  Ti  layer.  When  all  the  Ti  is 
consumed,  Si  begins  to  react  with  the  TiA13  layer  to  produce  the  ternary  phase.  Si 
then  migrates  from  the  contact  and  diffuses  into  the  Al  layer  while  A1  moves  into 
the  area  of  the  removed  silicon.  The  ternary  layer  then  appears  to  extend  beyond 
the  contact  window  as  the  Si  saturated  Al  provides  Si  to  convert  the  Tial3. 
Titanium  is  a  good  diffusion  barrier  for  Silicon  as  long  as  the  titanium  layer  in  not 
totally  consumed,  it  is  an  effective  bamer  to  AlSi  interdiflfusion.  It  has  been  found 
that  stuffing  the  refractory  metal  titanium  with  oxygen  or  nitrogen  is  seen  to 
greatly  improve  the  barrier  properties.  Titanium  nitride  has  been  used  as  a 
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salicided  CMOS  technology.  The  A1  /  TiN  /  TiSi2  /  Si  and  A1  /  TiN  /  Si  systems  are 
both  stable  up  to  annealing  temperatures  of  SSO^C  for  30  min  anneal.  Significant 
interdiffusion  occurs  at  BOO^C.  Therefore,  interdiffusion  of  aluminum  has  become 
an  intrinsic  defect  with  the  barrier  and  contact  material  in  the  salicided  structures 
since  processing  temperatures  range  from  700^  -  900®C. 

As  the  device  dimensions  are  scaled  to  the  submicron  regime,  the  reduction  in 
the  dimensions  of  metallization  film,  juction  depth,  and  contact  size  affecting  the 
electro  and  stress  migrations  leading  junction  leakage  currents  and  increased 
contact  resistance.  Even  at  moderate  current  densities,  the  bulk  self  dififusivity  of 
aluminum  is  sufficient  to  cause  voids  and  open  circuits.  In  an  effort  to  improve  the 
electro  and  stress  migrations,  several  met^s  such  as  Cu,  Si,  Ti,  Pd,  etc  have  been 
added  as  impurites  to  pure  aluminum.  Allowing  of  aluminum  decreases  the  grain 
boundary  self  dififusivity.  The  main  conducting  layer  improves  electromigration. 
Addition  of  cooper  primarily  reduces  the  number  of  hillocks,  improves  reflectively, 
increases  hardening  etc.  It  improves  electromigration  resistance  by  blocking  grain 
boundary  diffusion  paths  while  Si  addition  eliminates  A1  penetration  into  the  Si 
substrate  at  the  contact  areas  allowing  to  decrease  contact  resistance.  At  the 
present  time,  A1  alloys  containing  lwt%  Si  and  .1  -  .5wt%  Cu  are  widely  used  in 
submicron  CMOS  technologies.  The  .5wt%  of  Cu  addition  improves  the  uniformity 
of  the  grain  structure  and  further  increase  of  Cu  impurity  content  degrades 
electromigration  strength.  So  it  is  important  to  xmderstand  the  interfacial  reactions 
of  these  layered  structures  while  studying  the  effect  of  implanted  aluminum  alloy 
layers  with  underlying  barrier  and  contact  layers.  The  dopants  As,  P,  or  BF2 
distribute  themselves  in  the  shallow  jimctions  during  rapid  thermal  processing  or 
rapid  thermal  nitridation  and  also  believed  to  be  interacting  with  the  new 
intermetallic  phases  in  the  salicided  structures  effect  the  sheet  resistance, 
reflectively,  and  microstructures  of  the  film.  So  we  present  our  results  of  our 
experiments  on  Al-Si-  Cu  /  TiN  /  TiSi2  /  Si  and  Al-Cu  /  Tin  /  TiSi2  /  Si.  Unlike  the 
earlier  studies,  this  work  aims  to  understand  how  multilayered  metallizations  in 
salicides  behave  during  heat  treatments,  and  their  impact  in  the  submicron 
technology. 

EXPERIMENTAL 

Salicide  films  are  formed  on  Silicon  by  ion  implantation  and  subsequent  annealing. 
Titanium  nitride  films  were  sputter  deposited  in  the  TiSi2.  The  thickness  of  the 
TiSi2  and  TiN  layer  were  BOO^C  each.  Al-0.5wt%Cu  -  lwt%Si  film  was 
subsequently  sputter  deposited  on  the  top  of  the  TiN  layer.  Since  conventional 
techniques  do  not  satisfy  via  filling  requirements  for  submicom  devices.  MX214IOE 
bias  sputtering  metallization  techniques  are  used  for  filling  and  planarization  of 
aluminxam  by  taking  advantage  of  low  melting  point  of  aluminum  high  target  power, 
substrate  bias  and  substrate  temperature  to  enhance  aluminum  mobility,  in  order 
to  achieve  uniform  grain  size  distribution.  The  samples  were  subsequently 
annealed  in  pure  nitrogen  for  30  minutes  at  temperatures  of  400°  to  550®C.  The 
sheet  resistance  of  the  samples  were  measured  with  a  PROMETRIX  four  point 
probe  apparatus.  SIMS  depth  profiles  were  obtained  using  a  camera  IMA-35  ion 
microanalyser.  X-ray  diffraction  data  were  recorded  using  a  Rigaku  glancing  angle 
X-ray  diffractometer  with  a  monochromatic  Cu  radiation.  Cross  sectional  TEM 
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analysis  were  done  using  Phillips  120  keV  transmission  electromicrograph. 
Specimens  for  cross  sectional  TEM  analysis  were  obtained  by  means  of  mechanical 
polishing  and  ion  milling.  The  surface  orphology  was  obtained  using  SEM  and 
hillocks  were  cotmted  per  picture.  In  order  to  see  the  reflectively  changes  with 
respect  to  hillock  formation,  conventional  optical  microscopes  were  used.  AES 
analysis  was  used  to  characterize  the  A1  alloy  layer.  Most  of  the  experimental  data 
were  obtained  from  AT&T  Bell  Laboratories  ,  Allentown,  PA  18103, 

RESULTS 

We  begin  this  section  with  our  sheet  resistance  and  reflectively 
measurements.  Figures  2  and  3  show  sheet  resistance  of  the  n-Al-Ca  and  n-Al-Si 
-Cu  samples  as  a  ftmction  of  annealing  temperatures.  The  sheet  resistance  of  the 
samples  annealed  below  450^0  is  about  80%  of  that  A1  deposited  sample.  The 
sheet  resistance  increases  in  between  400^0  -  SSO^C.  The  increase  of  sheet 
resistance  is  significant  at  525^0.  Also,  both  n-Al-Cu  and  Al-Si-Cu  show  more 
increase  in  sheet  resistance  compared  to  P-Al-Cu  and  P-  Al-Si-Cu.  The  reflectance 
is  umform  up  to  400®C  anneal  and  decreases  drastically.  Also,  the  formations  of 
hillocks  has  increased  in  temperatures  between  400^0  -  550^C  which  in  turn 
affected  the  reflectance.  Figure  3  shows  the  reflectance  versus  temperature  plot  of 
the  samples.  As  the  hillocks  increased  reflectively  decreased  in  the  temperature 
range  of  400^  -  SSO^C,  Table  1  shows  the  number  of  hillocks  formed  during 
different  heat  treatments  for  the  samples  studied.  Figure  4  shows  the  SIMS  depth 
profiles  of  samples  annealed  at  400®C  -  550®C.  SIMS  intensity  of  each  element  is 
correlated  with  the  total  ion  density.  At  550®C,  the  temperature  at  which  the 
increase  in  sheet  resistance  is  observed  Ti  and  N  move  significantly  in  the  A1  layer 
towards  the  surface  A1  and  Si  move  in  the  TiN  layer.  In  addition  to  the  TiN 
thickness,  A1  moved  in  to  Si  substrate  at  the  temperature  SSO^C.  The  TiA13  phase 
is  not  identified  in  the  temperature  400®C  since  copper  addition  reduced  the 

interdiffusion.  However,  at  the  temperature  550®C,  weak  X-ray  Diffraction  signal 
have  been  identified. 

CONCLUSIONS 

The  interdiffusion  of  the  A1  -  .5wt%  Cu  -  lwt%Si/  TiN  /  TiSi2  /  Si  due  to  the 
different  heat  treatments  at  temperatures  400®C  and  550®C  in  nitrogen  was 
investigated  using  XRD,  SIMS,  AES,  RBS  and  TEM,  The  inter-diflusion  resulted  in 
increase  of  sheet  resistance,  decrease  of  reflectively  and  increase  of  hillocks. 
Titanium  started  to  diffuse  into  the  A1  alloy  layer  at  400®C.  However,  the  Ti 
reaction  product  could  not  be  detected  by  the  5QID  and  TEM  at  400^0.  Titanium  on 
the  otherhand,  n  diffused  and  reacted  with  A1  at  550®C.  Copper  segregation  is 
more  at  the  film  surface  or  film  substrate  interface  in  n  type  material  than  the  p+ 
type  material.  Also  the  increase  of  sheet  resistance  is  more  in  n  type  samples  than 
the  p+  samples. 
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Classical  and  Quantum  Statistical  Physics:  An  Exact  Approach 


U.  F.  EdgaJ  and  D.  L.  Huber 
Department  of  Physics 
University  of  Wisconsin-Madison 
Madison,  Wisconsin  53706 


A  scheme  analogous  to  that  developed  for  ’efficiently’  and  ’exactly’  determining 
the  free  energy  of  classical  fluids  with  realistic  many-body  interaction  potentials  at 
arbitrary  densities  and  temperatures  has  been  evolved  for  quantum  systems.  We 
present  in  this  paper,  a  brief  outline  of  the  scheme.  The  case  for  clsissical  systems 
which  WM  recently  published  is  reaidily  shown  to  be  a  limiting  csae  obtainable  from 
the  more  general  quemtum  treatment. 


I.  INTRODUCTION 

The  free  energy  of  equilibriTim  systems  is  said  to  be  characterized  by  the  di¬ 
mensionless  parameter  e.  By  employing  an  exact  integral  recursion  relation  for  the 
partition  function,  a  highly  nonlinear  ’differential  algebraic  equation’  [1]  governing 
e  is  formulated.  An  important  parameter  governing  the  nonlinear  equation  involves 
the  function  termed  the  ’general  point  process  neeirest  neighbor  probability  density 
function’(NNPDF)[2].  NNPDF’s  are  used  to  replace  radial  distribution  functions 
(and  more  generally,  n-body  correlation  functions)  for  describing  the  microstructure 
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of  equilibrium  systems.  They  may  also  be  exactly  formulated  in  terms  of  the  e  param¬ 
eter.  In  quantum  systems,  the  analogue  of  the  NNPDF  is  the  QNNPDF  (Q  standing 
for  Quantum),  e  is  determined  from  the  governing  differential  algebraic  equation  via 
an  iterative  scheme.  The  case  for  classical  systems  presented  in  detail  elsewhere[2]  is 
shown  to  resTilt  from  the  present  quantum  treatment  as  a  limiting  case. 


n.  SUMMARY  OF  RESULTS 

The  summary  given  in  this  section  is  elaborately  presented  in  a  paper[3]  presently 
under  preparation  for  submission  (for  publication).  The  canoniced  partition  function 
of  a  system  of  N  identical  classical  or  quantum  (fermi  or  bose)  particles  in  a  volume 
V  may  be  written  without  loss  of  generality  as: 

Z{N,V)  =  {enV)^IN\  (1) 


where  ClV  may  be  seen  as  a  ’one  particle  phase  space’,  fl  is  connected  with  the 
volume  in  ’interned  coordinate’  space  of  the  particles.  For  particles  with  no  internal 
structure,  f2  =  A“^  where  A  is  the  thermal  wavelength  (/3/i^/27rm)^/^;/3  =  l/AiflT,  T 
is  temperature,  kg  is  Boltzmann  constant,  h  is  planck’s  constant  and  m  is  mass  of  a 
particle.  The  following  relation  is  found  to  hold  for  Z(N,V): 

Z{N,  V)  =  r  Z{N  -  1,  Vi)PndVi  (2) 

Jo 


where  for  quantum  systems 


(-^1  >  X\ ,  ... ,  Xn)  O  /  Vr  tV  \77n  J  v 

fxr  (V  v\  An)lli_iaAi 

(n  ^  10) 


(3) 


Wn+\{X[,  Xi^ Xr^  is  the  ’slater  sum’[4]  for  n-|-l  particles  located  at  X(,.Yi,...,Xn. 
The  coordinate  X[  is  situated  on  the  boundary  of  volume  V\  and  is  taken  as  the 
origin  (The  boundary  of  voltime  Vi  may  usually  be  eissumed  to  be  locally  flat  ’almost 
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always  ).  Xi  is  the  coordinate  (relative  to  the  origin)  of  the  i-th  nearest  neighbor 
to  the  origin.  is  the  general  point  process  QNNPDF  describing 

the  probability  that  the  first  neaxest  neighbor  to  the  origin  is  located  (relative  to 
the  origin)  at  Xi ,  the  second  nearest  neighbor  to  the  origin  is  located  (relative  to  the 
origin)  at  X2  etc.  The  term  general  point  process  implying  a  paxticle  is  not  necessarily 
situated  at  the  origin.  We  may  write: 

Xn)  =  Kexp{~i,rll3P')Wn{Xx, je„) 

(n  ^  10)  (4) 

where  r„  is  the  radial  part  of  and  P  is  the  system’s  pressure,  is  a  normalization 
constant.  Eqn8.(l)-(4)  are  exact  relations  which  have  couterparts  for  classical  systems. 
In  the  classical  limit  where  the  average  distance  between  paxticles  is  »  A  (usually 
applicable  under  weakly  degenerate  conditions [4]),  ...,  where 

is  the  potential  energy  of  the  n  particles  at  With  tHs,  Eqns.(3)  and  (4) 

are  easily  shown  to  reduce  to  their  classical  counterparts [2].  The  validity  of  Eqns.(2) 
to  (4)  as  ’exact’  equations  for  the  ’exact’  determination  of  e  is  made  possible  by  the 
fact  that  the  slater  sum  possesses  a  ’product  property’  which  is  a  generalization  of 
that  usually  stated  in  the  literature[4];  ie:  If  we  may  write  W„  =  e*"",  then  it  is  true 
that: 


Wat  = 

where  |  wn,  |,|  wn,  >j  /  |;Wi,iV2  »  1  and  iVi  +  iVz  =  iV  (It  is  assumed 
we  have  two  ’clusters’  of  coordinates).  We  may  employ  (1)  and  (2)  (after 

differentiation)  to  arrive  at  the  differential  algebraic  equation  (assuming  T  constant): 

(5) 

p  is  the  particle  number  density  N/V.  The  equation  of  state  of  the  system  is  easily 
obtained  as: 
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where  4>  =  0P' •  In  the  full  version  of  this  paper[3],  we  briefly  describe  how  (5)  and 
(6)  may  be  solved  iteratively.  The  efficiency  of  computation  depends  mainly  on  how 
accurately  the  slater  sum  may  be  evaluated.  The  slater  sum  (the  diagonal  matrix 
element  of  the  Boltzmann  operator  in  coordinate  representation)  may  be  written 
(after  employing  a  complete  set  of  states)  as: 

WJ^X)  =  ^  ^  <  X  1  m  ><  m'  1  X  ><  m  I  |  m  >  (7) 

m' 

X  stands  for  the  coordinates  ;  m  and  m'  stand  for  occupation  numbers  of 

a  set  of  basis  states.  In  second  quemtized  formalism,  Hn  may  be  written  as  sums  or 
integrals  (with  infinite  limits)  over  terms  involving  creation  and  annihilation  opera¬ 
tors.  We  may  readily  show  that  the  Boltzmann  operator  may  be  accurately  written 
as  Y!i=o{-l3Hn)^ li\  (for  arbitrary  temperatures),  where  k  is  some  finite  number  which 
may  always  be  chosen  large  enough.  Also,  it  may  be  shown[3]  that  the  limits  of  the 
sums  in  (7)  (including  the  limits  of  the  sums  or  integrals  introduced  by  expressing 
Hn  in  second  quantized  formalism)  may  be  chosen  to  be  fimte  withbut  limiting  the 
accuracy  in  the  evaluation  of  (7).  Hence  we  may  expect  that  given  a  sufficiently  pow¬ 
erful  computing  environment,  (7)  may  always  be  accurately  evaluated  for  a  variety  of 
configurations  X  at  reasonable  ’computing  cost’. 

m.  REMARKS 

We  observe  that  our  above  scheme  first  reduces  a  many  problem  to  a  ’few  body’ 
problem  involving  computations  with  n  ~  10  particles;  the  method  of  second  quan¬ 
tization  (which  may  be  said  to  be  weU  adapted  for  accurately  handling  few  body 
problems)  is  then  employed  in  the  computation  of  Wn  (which  provides  probably  the 
most  computationally  intensive  aspect  of  our  scheme).  In  which  case,  the  ’’power”  of 


second  quantization  as  a  computational  tool[4]  (book-keeping  device)  is  fully  realized 
when  used  in  conjuction  with  our  scheme.  It  is  interesting  to  note  that  NNPDF’s  and 
QNNPDF’s  provide  fundamental  ways  alternative  to  ’n-body  correlation  functions’ 
for  describing  microstructure.  Presently  the  above  scheme  is  being  applied  to  the 
ideal  fermi  and  bose  systems  at  all  densities  and  temperatures.  Also,  multicompo¬ 
nent  systems  are  under  development  and  these  are  found  to  have  the  same  general 
mathematical  structure  as  that  for  single  component  systems.  Additionally,  the  sim¬ 
ulation  method  required  to  realize  our  scheme  is  under  development.  We  hope  by 
the  above  that  various  properties  of  material  systems  (including  liquid  He*,  high  Tc 
superconductors  etc.)  may  be  investigated  with  renewed  vigor. 


122 


references 


1  ZwiUinger  D.,  Handbook  Of  Differential  Equations  (2nd  edition)  Academic  Press: 
Boston  (1992) 

2  Edgal  U.  F.,  J.  Chem.  Phys.,  24,8179(1991);  ibid,  94,  8191(1991) 

3  Edgal  U.  F.  and  Huber  D.  L.,  ’Quantum  Statistical  Physics:  An  Exact  Approach’ 
( 1 993)-U  npublished 

4  Hu«ng  K.,  (2nd  edition)  John  Wiley;  Ne»  York  (1987); 

Negele  J.  W.  and  Orland  H.,  QuantnmM anvParticleSystms,  Addison- Wesley; 
Redwood  City,  California  (1988) 


123 


Characterization  of  short  pulse  laser-produced  plasmas  at  the  Lawrence 
Livermore  National  Laboratory  Ultra  short-Pulse  Laser'*' 

Ronnie  Shephard,  Dwight  Price,  Bill  White,  Susana  Gordan'*"*',  Albert 
Osterheld,  Rosemary  Walling,  William  Goldstein,  and  Richard  Stewart 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550. 

Ultra-short  pulse  (USP)  laser-produced  plasmas  have  opened  a  new, 
exciting  regime  in  plasma  physics  [1-3].  The  USP  experiments  at  LLNL  have 
focused  on  the  study  of  high  energy-density  matter.  The  K-shell  emission 
from  porous  aluminum  targets  is  used  to  infer  the  density  and  temperature 
of  plasmas  created  with  800  nm  and  400  nm,  140  fs  laser  light.  The  laser  beam 
is  focused  to  a  minimum  spot  size  of  5  pm  with  800  nm  light  and  3  pm  with 
400  nm  light,  producing  a  normal  incidence  peak  intensity  of  =  10^  ® 
Watts/cm2,  A  new  800  fs  x-ray  streak  camera  is  used  to  study  the  broadband 
x-ray  emission.  The  time  resolved  and  time  integrated  x-ray  emission 
implies  substantial  differences  between  the  porous  target  and  the  flat  target 
temperature. 

* 

1.  J.A  Cobble,  et  al.,  Phys.  Rev.  A.  39, 454  (1989). 

2.  M.M.  Murname,  H.  Kapteyn,  and  R.  Falcone,  Phys.  Rev.  Lett.  62, 155  (1989). 

3.  H.M.  Milchberg,  R.R.  Freeman,  S.C.  Davey,  and  R.M.  More,  Phys.  Rev.  Lett. 

61,2364  (1988). 

■•■Lawrence  Livermore  National  Laboratory  under  contract  number  W-7405- 
ENG-  48. 

■•‘■'■physics  Department,  University  of  California,  Berkeley. 
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Why  the  Speed  of  Light  is  Reduced  In  a  Linear  Dielectric  Medium 


Mary  B.  James 

Department  of  Physics 
Reed  College 
Portland,  Oregon 

Maxwell's  electro  dynamic  theory  predicts  that  the  speed  of  a  plane 
wave  in  a  linear  dielectric  medium  is  reduced  from  its  speed  in  vacuum  by  a 
factor  of  n,  the  index  of  refraction  of  the  medium,  which  is  closely  related  to 
the  material's  electric  susceptibility.  While  Maxwell's  description  of  this 
phenomenon  is  straightforward,  it  does  little  to  illuminate  the  mechanism 
responsible  for  the  reduced  speed  of  propagation.  I  describe  a  mechanistic 
model  in  which  the  main  molecular  dipoles  induced  in  the  material  by  a 
plane  wave  travelling  at  speed  c  "conspire"  to  produce  a  single  plane  wave 
travelling  through  the  medium  at  the  reduced  speed,  c/n. 
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Recent  Results  of  Research  on  1.3  GHz  Annular  Electron  Beam  Powered 
Multi-Gigawatt  Microwave  Amplifier 

Walter  R.  Fayne,  Kyle  J.  Hendricks,  Michael  D.  Haworth"'’,  Robert  C.  Platt"*", 
James  Wells"*",  Thomas  A.  Sperncer,  Lester  A.  Bowers,  Moe  J.  Arman^' 
Raymond  W.  Lemke*,  Michael  Mazarakis*,  M.  Collins  Clark* 

Electromagnetic  Sources  Division  (PL/WSR) 

Phillips  Laboratory 
Kirkland  AFB,  NM  87117-6008 

Recent  experimental  and  computer  simulation  results  of  our  annular 
electron  beam  amplifier  will  be  presented.  Our  research  includes  the 
generation  of  a  annular  electron  beam  using  the  IMP  pulser  (500  KV,  5  Q,  300 
ns),  the  modulation  of  this  beam  by  a  two  cavity  relativistic  Klystron 
amplifier,  and  the  radiation  of  the  microwave  power.  Our  work  also  includes 
the  measurements  of  power  used  to  modulate  the  electron  beam,  the 
resulting  modulating  voltage,  the  resulting  electron  beam  modulation,  and 
the  excitation  of  the  idler  (second)  cavity.  Included  in  our  work  is  the  design 
and  construction  of  a  1.5  Tesla  pulsed  magnet  to  constrict  the  beam,  and  the 
recently  revised  design  of  a  more  efficient  cathode.  We  will  also  present 
comparisons  of  our  experimental  results  with  the  computer  simulated  results 
from  the  2  1/2-D  code  MAGIC  and  the  1-D  code  RKA. 

■•■SAIC,  Albuquerque,  New  Mexico 

^National  Research  Council  Senior  Fellow 

*Sandia  National  Laboratory,  Albuquerque,  New  Mexico 
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Soft  X-Ray  Lithography  and  Nanostructures  Program 
Keith  H.  Jackson 

Center  for  X-ray  Optics 
Lawrence  Berkeley  Laboratory 
Berkeley,  CA 

As  the  semiconductor  industry  begins  to  consider  the  challenge  of 
lithography  and  pattern  transfer  at  critical  dimensions  below  0.25  pm,  it  is 
becoming  increasingly  clear  that  the  development  of  ultra  high  precision  x- 
ray  optical  elements  will  be  essential.  The  fabrication  issue  for  diffraction 
limited  imaging  optics  at  a  wavelength  of  130  A  can  be  generalized  with  the 
following  statement,  "If  the  optical  figure  can  be  measured  the  optic  can  be 
fabricated."  This  requires  a  technique  that  can  measure  the  surface  figure  of 
complex  aspheric  optics  with  a  precision  better  than  a  tenth  of  a  wavelength, 
i.e.,  less  than  one  nanometer  (10  A).  The  most  powerful  optical  testing  tool  of 
surface  figure  is  interferometry.  The  scale  of  precision  in  interferometry  is 
determined  by  the  wavelength  of  the  source.  The  need  to  perform  these 
measurements  at  soft  x-ray  wavelengths  is  clear.  The  Advanced  Light  Source 
(ALS)  finishing  completion  at  LBL  is  the  ideal  source  for  interferometric 
measurements.  The  low  emittance,  relativistic  electron  beam  in  the  ALS 
storage  ring  will  provide  partially  coherent  radiation  at  soft  x-ray 
wavelengths.  Undulators,  periodic  magnetic  structures,  will  increase  the 
intensity  of  the  radiation  and  narrow  the  bandwidth,  providing  soft  x-ray 
source  of  extraordinary  brilliance. 

The  Center  for  X-ray  Optics  at  the  Lawrence  Berkeley  Laboratory  is 
currently  constructing  two  x-ray  beamlines  at  the  ALS  that  are  dedicated  to 
optical  metrology,  and  a  Nanostructures  fabrication  laboratory.  The 
processing  lab  will  include  a  direct  write  e-beam  lithography  system,  a 
Nanowriter,  that  will  allow  the  fabrication  of  reflective  masks,  zone  plates, 
and  device  structures  with  feature  sizes  as  small  as  0.05  pm  over  an  8"  wafer. 
This  development  of  soft  x-ray  optical  metrology,  advanced  e-beam  writing, 
and  processing  technology  will  allow  soft  x-ray  projection  lithography,  for 
pattern  transfer  in  the  0.1  pm  range  and  beyond. 
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Semianalytic  Methods  for  Slater-Type  Orbitals 

Andrew  Jackson,  Jay  Jackson,  Babak  Etemadi  and  Herbert  Jones 


Department  of  Physics 
Florida  A&M  University 
Tallahassee,  FL  32307 

It  would  be  desirable  to  use  Slater-type  orbitals  (STOs),  characterized  by 
exp(-x),  rather  than  Gaussian-type  orbitals  (GTOs),  characterized  by  exp(-x2), 
basis  sets  in  all  calculations  involving  molecules  and  their  interactions.  This  is 
because  STOs  represent  electron  clouds  surrovmding  nuclei  with  greater  fidelity. 

The  difficulty  with  STOs  are  computational:  when  molecular  integrals  are 
derived  as  formulas,  they  are  subject  to  enormous  cancellation  errors;  when  they 
are  set  up  with  Fourier  transforms  for  numerical  integration,  pathologies  and 
singularities  often  appear. 

Our  approach  uses  computer  algebra  to  generate  a  C  matrix  for  use  in  the 
Lowdin  alpha-function  method,  and  then  performs  a  Gauss-Legendre  numerical 
quadrature.  This  semi-analytic  method  is  being  initiated  by  finding  overlap 
integrals  over  s  orbitals  with  an  application  to  the  hydrogen  molecular  ion,  H2''' 
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Diamond/Diamond'Like  Thin  Film  Growth  and  the  Alkanes: 
Methane,  Ethane,  and  Butane 

■  Elvira  Williams,  Johnnie  Richardson  Jr.,  Donald  Anderson,  and  Kristen  M.  Starkey 

Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 
Greensboro,  North  Carolina  27411 


Abstract 

Determination  of  synthesis  conditions  for  diamond  (and  diamond-like)  thin  films 
from  various  hydrocarbon  plasmas  is  an  important  aspect  of  Plasma  Enhanced 
Chemical  Vapor  Deposition  (PECVD)  research.  In  the  current  work,  PECVD  synthe¬ 
sis  of  diamond  was  investigated  in  radio  frequency  (Rf)  generated  plasmas  of  alkane 
and  atomic  hydrogen  mixtures.  The  alkanes,  methane  (CH^),  ethane  (C2H0)  and  bu¬ 
tane  (C^j^q)  were  used. 

•  Some  experimental  resvilts  related  to  the  study  of  diamond/diamond-like  thin 
films  on  imheated,  unetched  p-tjrpe  silicon  (100)  substrates  produced  in  different  al¬ 
kane  plasmas  are  reported.  The  relationships  between  film  index  of  refraction,  Rf 
power,  density,  molecular  mass,  hydrocarbon  structure,  hydrocarbon  flow  rate  and 
hydrogen  flow  rate  were  determined.  Index  of  refraction,  measured  with  an  elhp- 
someter  was  used  to  characterize  the  films 
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Diamond  Thin  Films  Synthesis  in  a  Non-Hydrogenated  Methane  Plasma  on 
Unheated,  Unetched  Si(lOO)  N-type  Substrates 

Elvira  Williams,  Johnnie  Richardson  Jr.,  Donald  Anderson  and  Bryan  Brown 

Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 
Greensboro,  North  Carolina  27411 

Abstract 

The  growth  of  diamond  thin  films  on  various  substrates  by  a  variety  of  methods 
is  an  important  part  of  diamond  thin  film  research.  Both  scientific  and  practical 
application  are  almost  limitless. 

Successful  diamond  growth  experiments  typically  employ  a  hydrogen-hydrocar¬ 
bon  gas  mixture.It  is  widely  accepted  by  researchers  in  the  diamond  thin  film  field 
that  such  films  can  be  grown  in  hydrogenated  methane  plasmas.  The  objective  of 
this  work,  however,  is  to  explore  the  possibility  of  growing  diamond  films  using  the 
hydrocarbon,  methane,  in  the  absence  of  atomic  hydrogen. 

Diamond  thin  films  were  synthesized  in  a  non-hydrogenated  methane  plasma 
using  the  technique  of  Plasma  Enhanced  Chemical  Vapor  Deposition  (PECVD).  A 
series  of  experiments  were  performed  to  characterize  the  film.  Measurements  of 
index  of  refraction,  composition,  thickness  and  mass  were  made  for  these  films. 
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A  Study  of  Butane-Hydrogen  Plasma  Enlianced  Chemical  Vapor  Deposition 

of 

DiamonThinFilms  on  Unheated,  Unetched  Silicon  Substrates 

Elvira  Williams,  Johnnie  S.  Richardson,  Jr.,  Donald  Anderson, 

Vernon  Simmons 

Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 
Greensboro.  North  Carolina  2741 1 


ABSTRACT 


Plasma  Enhanced  Chemical  Vapor  Deposition  (PECVD)  diamond-growing 
techniques  have  been  used  to  grow  very  thin  diamond  films  on  materials  since 
the  early  1950's.  However,  the  matter  gained  little  attention  for  the  next  thirty 
years.  The  primary  reason  was  that  It  was  thought  to  have  little  commercial 
fxjtentlal  due  to  the  extremely  slow  growth  rate.  In  the  past  few  years. 
Increasing  attentionn  has  been  drawn  to  the  special  form  of  carbon.  In  1990, 
diamond  was  named  "Molecule  of  the  Year"  by  Science  magazine.  The  thin 
films  have  been  named  i-carbon,  diamond-like,  hard  carbon  or  a-C.  Some  of 
the  methods  used  to  produce  the  film  are  evaporation,  sputtering,  ion  beam 
and  plasma  deposition.  Ion  beams  and  plasma  deposition  produce  films  whose 
properties  resemble  those  of  diamond.  If  made  under  apprprlate  deposition 
conditions,  the  films  are  extremely  hard  and  electrically  insulating. 

Very  important  aspects  of  diamond  deposition  research  involves 
investigating  diamond  films  characterization  and  deposition  on  various 
substrates.  The  tremendous  excitement  of  growing  diamond  is  largely 
generated  by  the  wide  variety  of  applications.  A  thin  coating  could  be  placed 
on  a  materials  which  would  enhance  their  longevity;  an3d:hing  from  tools  and 
lenses  to  turbine  blades  could  be  coated.  In  these  experiments,  the  Si  (100) 
substrate  was  unheated  and  unetched. 

Not  heating  the  substrate  allows  the  coating  of  all  materials  which  remain 
stable  at  room  temperature.  The  growth  process  is  more  convenient  when  the 
surface  of  the  substrate  is  unetched.  Experiments  under  these  conditions  were 
used  to  Investigate  the  growth  of  diamond  in  a  very  low  frequency  using  butane 
(C4H  jq)  as  the  reactant  gas  in  the  presence  of  hydrogen. 

A  series  of  experiments  were  performed  in  which  film  index  of  refraction, 
film  thickness,  and  mass  were  measured  cis  functions  of  radio  frequency  (Rf) 
power  and  butane  flow  rate.  The  hydrogen  flow  rate  remained  constant 
throughout  the  series  of  experiments.  Index  of  refraction,  measured  with  an 
eUlpsometer,  was  used  to  characterize  the  films. 


131 


A  Study  of  Plasma  Enhanced  Chemical  Vapor  Deposition  of  Diamond  Thin 
Films  on  Unetched,  Unheated  Aluminum  Substrates 

Elvira  Williams.  Johnnie  S.  Richardson.  Jr..  Donald  Anderson. 

John  M.  Brown.  Jr. 


Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 
Greensboro,  North  Carolina  27411 


ABSTRACT 

Diamond  thin  films  offer  a  wide  variety  of  unique  thermal,  chemical, 
electrical,  and  optical  properties.  These  diamond  thin  films  have  properties 
similar  to  bulk  diamond:  high  thermal  conductivity,  high  electrical  resistivity, 
extreme  hardness,  radiation  resistance,  large  index  of  refraction,  and 
ultraviolet  to  infrared  transparency.  In  the  past  few  years,  interest  in  diamond 
research  has  gained  much  momentum  due  to  the  potential  applications  of 
diamond  films. 

There  are  many  methods  and  techniques  to  produce  diamond  thin  films. 
Some  of  the  methods  are  evaporation,  sputtering,  ion  beams,  and  plasma 
deposition.  Of  these,  the  plasma  and  ion  beam  methods  produce  diamond  thin 
films  with  properties  most  resembling  those  of  bulk  diamond. 

The  primary  objective  of  the  current  investigation  is  to  deposit  from  a 
hydrogen -methane  plasma,  uniform  layers  of  diamond  on  unetched,  unheated 
aluminum  (metal)  substrates  using  PECVD  techniques. 

A  series  of  experiments  were  conducted  in  which  the  film  thickness,  film 
mass,  and  film  index  of  refraction  were  measured  as  a  function  of  the  Rf  power. 
The  methane  and  hydrogen  flow  rates,  substrate  chamber  position,  ^d  the 
deposition  times  were  kept  constant  throughout  this  series  of  experiments. 
Index  of  refraction  was  used  to  characterize  the  films. 
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The  Effects  of  Variation  in  Methane  Flow  Rates  on  the  Index  Of  Refraction 
of  Diamond  Thin  Films  on  Unetched  and  Unheated  Silicon  Substrates 

Dr.  Elvira  Williams 
Dr.  Johnnie  S.  Richardson  Jr. 

Donald  Anderson 
Stacey  Barrett 

L’Tonya  Jefferson 

Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 

Greensboro,  NC 

'  ABSTRACT 

Diamond  and  diamond-like  films  were  produced  from  methane  and  hydrogen 
mixture  by  a  capacidvely  coupled  radio  frequency  (R.f.)  reactive  plasma  deposition 
process  imder  particular  deposition  conditions.  The  index  of  refraction,  mass,  and 
thickness  were  determined.  Variations  in  the  index  of  refraction  could  be  correlated 
with  the  methane  flow  rate.  Characterization  of  the  diamond  films  were  done  by  the 
technique  of  ellipsometry.  The  deposition  of  diamond  and  diamond-like  thin  films  on 
variety  of  substrates  can  be  studied  by  varying  the  processing  parameters  and  using 
different  reaction  gases.  Thin  films  with  material  properties  similar  to  diamond  or 
graphite  can  be  deposited  using  the  technique  of  Plasma  Enhanced  Chemical  Vapor 
Deposition  (PECVD).  Diamond  is  most  noted  for  it's  extreme  hardness  and  high 
thermal  conductivity.  In  addition,  diamond  has  excellent  electrical  properties  such  as 
it’s  large  dielectric  constant  and  its  high  resistivity  which  are  important  in  some  o 
today's  emerging  technologies.  The  unique  structure  of  diamond  has  the  potential  of 
providing  improved  arc  resistance  and  dielctric  breakdown  strength,  which  are 
critical  properties.  Furthermore,  chemical  vapor  deposition  techniques  offer  the 
possibility  of  developing  conductive  n  and  p-type  diamond  films  which  would  be 
enabling  materials  for  these  new  technologies.  The  chemical  vapor  deposition 
technique  uses  hydrocarbon  gases  as  the  carbon  source.  The  general  reaction  is  that 
these  hydrocarbon  gases  react  via  the  supplied  Rf  energy  to  produce  solid  reactants 
products  and  gaseous  by-products.  Methane  and  hydrogen  were  heated  with  Rf 
energy  to  produce  diamond  and  some  gaseous  by-product.  Rf  plasma  deposition  is  a 
technique  for  the  controlled  hetergeneous  polymerization  of  thin  films  from  active 
species  produced  in  an  Rf  generated  glow  discharge.  The  index  of  refraction  and  the 
amount  of  deposited  mass  decreased  rapidly  while  the  film  thickness  increased, 
rapidly  as  the  flow  rate  of  methane  gas  was  increased  from  22.1  to  22.9  SCCM. 
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Deposition  of  Diamond/Diamond-Like  Carbon  Films  on  Unetched. 

Unheated  Silicon  (100)  Substrates  produced  in  an  Ethane  Plasma 

Dr.  Elvira  Williams 
Dr.  Johnnie  S.  Richardson  Jr. 

Donald  Anderson 
Katina  Wilson 

Department  of  Physics 

North  Carolina  Agricultural  and  Technical  State  University 

Greensboro,  NC 

ABSTRACT 

Diamond  and  diamond-like  films  were  produced  from  ethane 
mixture  by  a  capacitively  coupled  radio  frequency  (R.f.)  ^ 

process.  The  films  were  deposited  at  room  temperature  on  unetched  p  typ 
flOO)  substrates.  The  films  were  characterized  by  ellipsometer  . 

index  of  refraction.  The  Department  of  Physics  at 

Technical  (A&T)  State  University  is  currently  engaged  m  research  to  study  th 
dSSn  of  diamond  and  dimnond-like  thin  films  on  variety  of  substrates.  By 
varying  the  processing  parameters,  thin  films  with  materi^ 
diamond  or  graphite  can  be  deposited  using  Plasma  Enhanced  Che^c^  p 
Deposition  (PECVD).  Diamond  is  most  noted  for  it's  extreme  hardness  and  g 
Lrmal  conductivity.  In  addition,  diamond  has  excellent  electrical  ^ 

it's  large  dielectric  constant  and  its  high  resistivity  which  are  ipiport^t  ^  some  of 
today’s  emerging  technologies.  The  unique  structure  of  ^  a?e 

providing  improved  arc  resistance  and  dielctnc  breakdown  strength, 
critical  properties.  Furthermore,  chemical  vapor  deposition  techniques 
possibility  of  developing  conductive  n  and  p-type  diamond  films  which  would  be 
Labling  materials  for  these  new  technologies.  The  chemical  vapor  deposinon 
technique  uses  hydrocarbon  gases  as  the  carbon  source.  The  general 
these  hydrocarbon  gases  react  via  the  supplied  Rf  energy  to  produce  sohd  reactan 
products  and  gaseL  by-products.  Specifically  in  this  expe^ent,  ethane  ^d 
hydrogen  are  heated  with  Rf  energy  to  produce  diamond  and 
product.  Rf  plasma  deposition  is  a  technique  for  the  controlled  hetergeneous 
polymerization  of  thin  films  from  active  species  produced  in  an  Rf  generated  glow 
discharge.  The  index  of  refraction  decreased,  the  film  thickness  mcreased,  and 
amount  of  deposited  mass  increased  as  the  flow  rate  of  ethane  gas  was  mcrease 

from  19  to  23  SCCM. 
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"International  Physics  Olympiad" 
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Department  of  Physics 
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6:15  pm 
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Department  of  Physics 
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